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DIFFERENTIAL EQUATIONS 

BETWEEN TWO OR MORE VARIA' LES. 



SECTION I. 

DIFFERENTIAL EQUATIONS OF THE FIRST ORDER AND 

DEGREE. 

1. In that part of the Integral Calculus which relates 
to the integration of explicit functions of one variable, we 
have to determine the relation between y and or from the 
equation 

in the present portion, we have to determine it from the 
equation 

or to assign the relation between Wy y^ z (where zi& bl function 
of the independent variables m and y)^ or between a greater 
number of variables and their functions^ from the equation 

f{d„%y dyXy Zj Of J y) = 0, 

or from other equations in which a greater number of variablea 
and differential coefficients of higher orders are involved* 

2. A differential equation is said to be of the nf^ order,, 
when the differential coefficient of the highest order which it 
involves is the n^\ 



.' '*■ 



A differential equation of any order is said, moreover, to 
be of the first, second, &c., degree, when the differential co- 
efficient which marks its order, is raised to the first, second, &c., 
power. 

To integrate a difi^rential equation of any order, is to 
pass to the primitive equation between the variables and the 
constants, from which the proposed may have been derived 
by the process of differentiation. 

3. We shall begin with the simplest case^ viz. that of 
differential equations of the first order and degree, which will 
be of the form 

M + Nd^y = 0, 

M and N being functions of x and y. 

Every difiPerential equation of the first order and degree is 
either the direct derived equation of a primitive ; or it results 
from the combination of the derived equation with its pri- 
mitive, so as to eliminate a constant which enters in each only 
to the first power ; the former sort are called exact, the latter 
inexact. 

4. First, let u^f{w^ V)^^ be an equation between a 
and y, by virtue of which y is a function of ^; then, as is 
proved in the Difibrential Calculus, d^y is given by the equa- 
tion 

or, since d^g)U^ ^{3f)^9 ^^^ functions of w and y which we may 
represent by M and N^ . 

M + Nd^y « 0, 

a differential equation of the first order and degree, of which 
/(^» y)*^ ^* ^^ primitive or integral. Now if /(d?, y), 
besides other constants which are afiPected with w and y, 
contain a term + C independent of w and y ; this will not 
enter into M and N, having disappeared in differentiating; 
and if there be no such term, we may add it, and/(/p, y)^C^O 



is still a relation between m and y which satisfies the equa- 
tion 

under this form it is called the complete integral; and the 
constant C, which does not appear in the differenfial equation, 
is called the arbitrary constant ; if the integral did not contain 
such a term as + C, it would not be sufficiently general, and 
would be only a particular case of the complete integral. 

We shall presently give the test which every equation of 
this sort must satisfy, and the mode of integrating it. It is 
evident that no equation of the first order which is not of the 
first degree can be exact. 

5. Next, let Cx be another constant which enters to the 
first power in the equation 

then Cx will be affected with oo and y, and will consequently 
appear to the first power in 

M + Nd^y = ; 

and if a value of Cx be obtained from either of these equations 
and substituted in the other, the result will be 

Mx + JNTid^y^O, 

an equation of the first order and degree, involving all the 
constants which enter into /(a?, y) + C = 0, except Ci. Hence 
whilst the direct derived equation of 

/(*> y) + C = 0, viz. M + Nd^y * 0, 

does not involve the term C which is independent of uo and y, 
there will be as many other differential equations of the first 
order and degree that have f(po^ y) + C ^0 for their primi- 
tive, as it has independent constants entering only in the 
first power; if any constant enter in a dimension above the 
first, the differential equation obtained by eliminating it, will 
evidently not be of the first degree. 



There are two principal methods of integrating equations 
of this sort, which consist either in separating the variables, 
by substitution, or some algebraical process; or in restoring 
the factor which makes them exact. 



Exact Differential Equations of the First Order. 

6. Let yhea. function of <v determined by the equation 

«*=/(^5 ») = <>» 
then the equation which gives the value of d^y is 

rf(,) w + d^yU . d^y = 0, or JIf + Nd^y = ; 

the notation d^f^Uy d^^^Uy implying that these differential co- 
efficients are formed on the hypothesis that y is not a function 
of <r, i. e. that w and y are independent ; then, as proved in 
the differential Calculus, 

d^y^M = d^g^N. 

Conversely, an equation of the form M + Nd^y = being 
proposed in which M and N are functions of <v and y, if the 
condition 

d^)M = d(^)N 

(which is called the criterion of integrability) be satisfied, the 
equation results from the immediate differeptiation of an 
equation of the form /(a?, y)sO; if this condition be not 
satisfied, there exists no equation by the simple differentiation 
of which, the given equation can be produced. 

7* To integrate the exact differential equation 

M + Nd^y = 0. 

Let the equation from which it is derived be 

u^f{oB^ y) = 0; 



then d^^yu » M, d^^^^u » N^ and d^^^M « d^^^N\ 

denoting by F a function of y which may have disappeared, 
since M is the differential coefficient of u relative to a^ on the 
hypothesis that x and y are independent ; 

.-. d,Y^N-d^{J^.^M\, and F-/y|J\r-d,y)C(,,if)}, 

the complete integral involving one arbitrary constant. 

8* Obs. The equation F«j^f JV-rf^)(j[^jJlf)J will be 
absurd, unless the expression JV— d^f^iJ^M) be independent 
at w; therefore its differential coefficient with respect to w 
must vanish; 

must equal zero, which it does, since the criterion of inte- 
grability is supposed to be satisfied. 

9* As the simplest case of exact equations, we may first 
notice those in which the variables are separated ; they will 
be of the form 

X + Fd,y « 0, 

where JT denotes a function of w only, and F a function of 
y only; here the criterion of integrability is manifestly sa- 
tisfied, for ' 

djy)jr=rf^„F=0; 

and the complete integral is 

/,jr + /,Frf,y«C, or f,X^f,Y^a 

To this case may likewise be reduced the equation 

XY^^YXid^y^O, 



which becomes, when divided by ^|F], 

X r _, 

.*. sin * - + sin ' — =s sin — , 
a a a 



or sm 



Va a* a a?l a 



or ofy/cf -y* +y\/a*- af^-aCi 

at which we may also arrive, by multiplying the proposed 
equation by a;y, and integrating by parts, which gives 

or, since the part affected by the sign jC, viz. 

is equal to zero by the proposed, 

y \/c? - ^ + w^<^ - y^ = C 

Ex. 2. i+y4.y* + (i+ar + a?*) d,y « 
C (J7 + y + 1) e 2a?y + 07 + y - 1. 

10. The following are instances of the integration of 
exact differential equations by the method of Art. 7. 

Ex. 1. aa? + 6y + c + (bw + my + n) d,y =* 0. 



• • 



u s — + (by + c)w + Y9 



d^yj ur^hw •\' dyT^ bw + my + n ; 
/. dj^Y ^ my -i-n; 

log («w + \/a?* + y*) + C «= 0. 
2y - 9ixdgy (w + x/V-s^' 



Ex. 3. 






Ex.4, y* + «^ -K' - ^y) «<>y , p. 

y^ + Say it - a + o*^r 



Homogeneous Equations. 

11. We come next to the case of inexact equations, in 
which the variables are separable by substitution; of these 
the most important class is homogeneous equations. 

Let M + Nd^y » be a homogeneous equation, that is, 
one in which each of the functions ilf and N is or can be 
expressed by series of the form 

N = ayf'ar'f' + fiy'^x'"'' + yy'^cf''' + fee, 

the sum of the dimensions of w and y in each term of M 
and N being equal to r. 



8 



Let y " wx where x denotes a new function of x, then 

Hence, making these substitutions in the given equation, and 
dividing by af^y 

f{z) + <p{x) {x + wd^z\ « 0, 

1 d,» 

or - + ■ = 0, 

in which the variables are separated. Similarly, the variables 
may be separated by making a^yx; and the latter substi- 
tution will be more convenient when iNT is a more complicated 
expression than M. 

Hence it is easy to effect the separation of the variables 
in equations which are either homogeneous, or can be made 
homogeneous; besides these, the number of equations in 
which that separation is possible, is very limited. 

Ex. 1. sy^w + 2/»^ + y*rf,y = o. 



U Sy^a-^^a? 
Here -r^ = -^ — ; 

N tr 



y \yJ X )r* 



1 d,x 



1 x^d^x 
- ^ 1— 

w «r* + Sx^ + 2 



or - + o. ^ -0, 



or 



W V + 2 «*+!/ 



or »« + 2«* = C >/^* H- »*. 



S. (^-y)d,sr = a7 + y, tan"^ ? « log V -7^ ' 

5. wd^y - y = \/^-y*, sin-^ - = log 7; • 

6. ordyj? + y = \/a?* + y*, 2Cy-^C^-a^' 

7. \/y + (\/y - \/^) d*y = 0, 

12. In the following instances, the equations are not ho- 
mogeneous, but are made so by easy substitutions. 

Ex. 1. a<r + &y + c + (mw + wy + p) d,y = 0. 

Let aw -{-by + c :=^ z, mx + ny + p = t?, 

z and t9 denoting functions of a?; then 

a + 6d,y = d,«r, m + nd^y = d,t) ; 



10 
.*. m + ndfif = d,v . d^z « (o + bd^y) d,« ; 

J 

m - adgV , z 

or m» - «« + (6« - au) d,t> « ; 
this being homogeneous, assume v = xw^ w being a function 

of «y 

1 (an? -- 6)d,«7 

in which the variables are separated. 

Ex. 2. af^ (ay + bwd^y) « y" (ay + (i^d^y), 

ox (bar^fi^) — ^{af-aaf')^. 

y ^ 

Let a^^Xy y" = V, jir and v being functions of Xy 

m dgZ nd^y d^v nw ^ j . 

w X y ii my v 

.*. (bx — /3t>) - — d,v a= at) - a«, 
which is homogeneous. 

Ex. 3. d,y + aff^ai'^ + 6y*^ = will become homogene- 
neous by making y = x^-% the equation of condition between 
w, «, p, V, being (p + 1) (l - ?) = (m + 1) (l - w). 

Ex. 4. d^y^— ^ 



This may be written 



— + s 

y 

and therefore becomes homogeneous when x is written for -. 

y 
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Linear Equations of the First Order. 



yy 13. The next important class of inexact equations of 
^ the first order which admit of being integrated^ are linear 
equations, the general form of which is 

P and Q being functions of w ; they are called linear because 
they involve no power of y above the first. 

Assume y » vir, v and z being functions of Wy 

Now z being an indeterminate quantity, may be assumed so 
that the equation last written down may resolve itself into 
two others, each of which admits of the separation of its 
variables; to this end let 

vdgZ + Pxv = 0, 

* d X 

or, dividing by Vy d^z + P« = 0, or -^ + P == ; 

z 

.*. log z ^ — fgPj or » = e'f'^. 

The remaining part of the equation gives 

zdgV^^Qy or, substituting for «, rf^t? ^p Qc/***; 

and y = e-/'^UQe/-''+C}, 
the complete primitive involving one arbitrary constant. 

Obs. It is unnecessary to add a constant after perform- 
ing the integration indicated in- the equation z = e^f"^ ; for 
let z = e-/-^-*-^ « C^rf*^ ; then 

which is the same result as before, since CC^ is equivalent 
only to a single constant. 



12 
14. If we differentiate the result 

we get ef'^(d^y'\-Py)^ef'^Q, 

which shews that if we multiply the proposed equation by 
e^'^, each member is separately integrable; and this is the 
most convenient practical mode of integrating it. When it 
is once known that the factor which makes the equation in- 
tegrable is a function of w only, its value may be immediately 
found; for let it be denoted by JT; then 



Xdj^y + (Py - Q) ^ = is exact, 

X 

.-. logJr=j;P, or X=ef'^. 



.-. d,X^ d^^(Py - Q) X^ PX, or ~ = P, 



15. It must be observed that if the second member of 
the equation d^y + Py = Q be multiplied by any power of y, 
it is still reducible to the standard form of a linear equation 
of the first order. For suppose the equation to be 

d^y + Py^ Qy", 

then dividing both sides by y*, and multiplying by - (n - 1), 
we get 

'^.(-i-,)-^,(«-i)P=-Q(«-i). 

Hence the factor which makes both sides integrable is 
^-(•-i)/,p^ and the result is 



^--,-?3» /.P--ilog(l + «»)-log ' 



, + ^' •" ^ -- ' -' -^^v^T^' 
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.-. J- — = /t — ^+C=-T=== + C; 

.*. y =3 o J7 H- C \/l + a;*. 
Ex. 2. rf,y + y « ^y*, 






.-. -- e ^ + 1 + Cc**- 

5. d^y = a^y» - d?y, — = a;* + 1 + Ce/^. 

ft 



^ _ . , o sin 07 + cos a _ . 
o- d^y s a sm 07 + 6y, y = - o • — - + C/e*^. 



6 sin 07 + cos w 
7+^ 



Riccati*8 Equation. 

/ 16. There are certain cases of the equation 

d,y + 6y*= aa^ 

(called Riccati^s Equation, after the Mathematician who first 
considered it) in which the variables are separable. 
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First, letm = 0, then d^y^a-^h^y or — ^—— = 1, where 



the variables are separated. 



Wd„.l«„-,.^...„»o,.^, „W„. 



function of w\ 



1 , u hu^ a 

tP w Or Or 



.\ adgU =s a + M — 6«% 

where the variables are separated. 

1 u 
Thirdly, let «i = - 4, and assume y = i— + -^ > 

bw Or 



bar' a^ a^ baf^ a^ a?* a?** 



.•. oFdgU + 6i^* = a, 
where the variables are separated. 

17- Besides the above, the variables are likewise separable 

-4t 
m the cases when m = — : , i being any integer from to 

infinity; all which values of m evidently lie between and - 4. 

•^ 4l 
First, let m = — : , 

Assume y = 7 — V ^jr-l 

1 ^ d^u 

12 6 

^ 6^ arw a?*t** 
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.*. adding these together^ 



b d^u 
= aoT ^ , 

or a^d^fU + au^w^'^^ - 6 « o. 

1 m + 2 

Now, let^=««+3^ then d,w=:d,wd,« = (f» + 3)««+* d,M; 

a 6 *»+^ 

.•. d,w + w* = «'"'»+3, or d,w + fejw" a «,»"*». 

m + 3 m-i- 3 

— 4e m + 4 4i — 4 4(i — 1) 
But m B — ; , .•. — S3 = — — — sm.. 

2*-l w+3 2«-3 2(«-l)-l 

Hence by these substitutions the equation is transformed 
into another of exactly the same form, with i — I instead of i 
in the index of the variable in the second member. 

1 1 -i- 

Similarly, by substituting — + — -^ for w, and ss^^-^^ 

for %^ we shall transform the equation into another of the same 

form where nu « —r. ;^ — ^ ♦ and consequently, after i sub- 

2(i-2)-l ^ ^ 

stitutions the index of the variable in the second member will 
become zero, and the variables will be separated, 

— 4i 
Secondly, let m 



2i + l 



Assume y «= — » 



, dgU b 
then — + -7 = aaf% 



or - dgU + 6 = acfTiJ?^ 
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1 m 

let w « «*+* , then d^u = d,« d,^ t= (i» + i) »«+i d^u ; 



.". -(f» + !)«•+» d,w + 6«a«~+*w*. 



_f a ^ b — * 

or dgU + ;- u* = X "•+!, 

ii»4- 1 m + 1 



-4i -21 + 1 

but m = ; .-. m + 1 = — ; ; 

2«+r 2» + i 

m — 4i — 4-t 



m+1 -2i + l 2i-l 

Hence by these substitutions this case is reduced to the 

— isi 

former; and therefore when m = — : , the variables in the 

2i±l 

equation dj,y^bf^ = aar can be separated. It may be observed 
that the more general equation, d^y + 6j^ j/'"^ = aw^, is reducible 
to this form by putting w^ = «. 



/ • 



18. We shall now give some other instances of equations 
m which the variables are separable by particular substitu- 
tions. 



Ex. 1. a (wdt,y - y) = (4? + yd,y) y/ of + J^ - a*. 
When an equation contains the expressions 

yd^y H- a?, wd^y - y, V^-^^^ 

the introduction of polar co-ordinates will sometimes effect the 
separation of the variables ; that is, to assume 

^ s p cos d, y-p sin Oi 

p being supposed a function of 0, for then 

ydey + xdQW = pd^p^ aid^y - ydQW ^ p". 



mm^ammt^ 
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Hence the proposed equation, which considering w and p 
as functions of 0y may be written 



becomes a^ ^ p dep \/p* - 0% 

or o =s — '" V p* - a* ; 
P 

.\ aO « \/0* -a^ "-a sec"* ^ + C, 
' a 

or a tan"* - = \/ar ■\- fr - or —a sec"* ^ + C. 

Ex. 2, y - <»d,y -h d*y \/ {ax + 6y) - = 0. 



/ ^ 

or ydyOB - ^ + 'V («^ + 6y) - = ; 

let -* = », j» being a function of y, 
then ydyss - a? = y*dy«f ; 

where the variables are separated. 
Ex. 3. d^yfc ' "" ' ^ ^ 



" • - * 



Since this is satisfied by y = w, assume y = a? + i^'S 
.\ 1 - af"^cL» =! + —(- 2^;^"* - »"*), 

ot (<,jif =fe — , which is a linear ^quatibn ; 



18 



-5f 1 _** 
or 



or — --= - Le «*+ C. 
y — a? o* 

Ex. 4. (l - wy) dgy + y* + a^ « 0. 

Assume y = , so that x = 

\ -{- «% 1 - ^y 

, dgZ w 
then — + :; -^ = 0, 



where the variables are separated. 

. v^ n(l + y«)f 
Ex. 5. (y-^)d,y= '-^ . 

V 1 + ar 



assume y = ; 



• =0. 



,8 



(1 + s?) {% + ny/T\^) 1 + ^ 



Euler's Equation. 
19. To integrate the equation 

or, considering m and y as functions of a new variable t, 

dty y/^ + dt^\/F= 0. 

Let the function of t which expresses os be determined 
by the equation d^x = s/Xy and therefore that which ex- 
presses y by the equation d^y = -\/F; also let ^ + y = j99 
^ "" y = 9> i^ and 9 being functions of ^. 



^mfmm^ 
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Then since (d,«r)*«-A^; .-. 2dtX d^XM^diX^ or dfa^\d^X; 
similarly dfy^^dyY. 

Se 
« 6 + c (^ + y) + — {(a? + y)* + (a? - y)H 

4 

4 

Se 
„ 6 + cp+ — (p» + g*) + ^/p (p» + Sf)y 

and d|p. d,g=Jr-F=6(«-y)+c(»*-y)+e(«*-y*)+/(i»*-y*) 

«6g + epq + - ?(3p^+ 9*) + \fpq (p* + «*) ; 

1 ^ 

.-. d?p --dtp.dtqm-^^ -^fP^^ 
q % 



or d|{( — ) [»cd|p + 2/jpd,p; 

■• -^ = c' + ^i^+/p^ 

or \/^-\/F=(^-y)\/C + 6(a7 + y)+/(.i? + y)«, 

the integral required. The discovery of this integral, which 
is due to Euler, was of great importance, as being the first step 
towards the foundation of the Theory of Elliptic Functions. 

20. To integrate the equation 

\/l - c* sin* y\f + ^/l - c^ sin* <f> d^yf/ = o, 
or considering <j) and yff as functions of another variable /, 

\/l - c* sin* >// di<p + \/l - c* sin* d| >/r =s 0. 



so 

Let di0ai\/l-c*sin^0, and .-. d|>^a:-\/l- c*sin*>|f ; 
.-. (dj0)» = 1 - c* sin« ^, (df >/.)2 = 1 - c« sitf x/. ; 
/• <ii'^ + d|*>|/ ==' - i c* (sin 20 + sin 2 >/f), 
rfi*0-" dt*^/' "^ -^c*(8in20- sin2\/^). 

Let p«0 + >/f, q^^p-i^y 

.•. dj^p = - c* sin p cos g, 
di'g = - c*cosp sin ^, 

and dip.dtq^ i^t<t^y " O^t^Y = - c^ (sin^0 - sin* yff) 

= — (cos20 - G6s2\ff) s= - c* sinp sinf. 

/. dtp = C sin g ; similarly diq = Csmp; 

.'. \/l - c* 8in*0 =F \/l - c'sin*>/^ = Csin (0 =f >/^) 

is the integral of the proposed equation; which is only 
Euler^s equation under a different form. 

The equations y/l +jr* + v^l + a* d,f =* o, 

are immediately reducible to the above form, viz. 
y/l -^8in*^ + ^/l-^sin«0d<^>|^ = O; 

the former by making a? = tan ^ 0, y « tan ^ ^ ; the latter 
by making v^^cos0, \/y = cos\|^. 



21. The most natural way of obtaining the complete in- 
tegral of a differential equation of the first order, is to prepare 
it so that its first member may become an exact differential 
coefficient; for then we shall have only to integrate and add 
a constant. This preparation is always possible by means of 
a factor, when tlie equation is reduced to the form d,y + K=0. 
For let an equation /(at, y, C)~0 be resolved with respect 
to C, so that 

p 

.'. by differentiation, = P + Qd^y, or d,y + — - 0. 

Now the equation M+ Nd,y = may be put under the 
form d,y + K = 0, which agrees with the preceding, and may 
consequently be supposed to have arisen from the elimination 
of a constant between the primitive /(a?, y, C) = 0, and its 
immediately derived equation. On this supposition, therefore, 



„ P+Qd,y 
...d.y^K ^, 

or d,<p (.r, ff) = Q {d,y + K), identically. 
The second member therefore is an exact differential 
coefficient, which proves that there always exists a factor 
proper to render the expression d,y -I- K integrable. 

22. But although the existence of the factor in ei 
case is thus established, the investigation of it is usu 
attended with greater difficulties than the solution of 
original equation. 

For let P + ^d,y = be an exact differential equati 
and let », a function of x and y, be a common factoi 
i* and Q, so that P = Afz, Q = JV«, by the removal 
which, the equation is reduced to the inexact state, 
M + Nd^yOi 



then because P + Q.d,y =i is exact, 

*/,„/• = d,.jQ, or d,„(Jf«) = d,,,(JV«), 

or »rf(^jl/ + JfdjyjSf = adj,jA''+ A''d(,,«, 
or JVdj^a-JIfdj^^-ar^yjaf-rfj^Ar); 

an equation between as, y, «, and tbe partial differential co- 
efficients of «, for determining the factor x. The consideration 
of this equation in its general state must be reserved till we 
come to treat of partial differential equations of the first order ; 
but the foUowiog particular cases may be noticed. 

23. First, suppose that the factor is a function of only 
one of the variables ^i then d^y^x = o, and the equation be- 
comes 



which, being integrated, gives z ; for the hypothesis requires 
that the second member should be independent of y. 

Similarly, if the factor be a function of y only, it will 
result from the integration of 

of which the second member is independent of or. 
"ence, if in any equation M + Ifd^y s o we Gad 

ition of a only, or 

j{i<',.,N-d„M).r 

tion of y only ; the factors which may make it integrable 
jpectively e^-^, e^'^. 
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Ex. 1. djfff + (Py - Q) s= 0, the linear equation of the 
first order. 

This compared with M + Nd^y «* 0, gives 

.-. d(y)Jlf - d^^^N^Py and 1 {d^^^M -- d^,^N) « P, 

a function of a: only ; therefore the factor is e^'^. 

Ex. 2. ^ + (1 - ^y) d,y = 0, 

M^y\ JV = 1 - a?y, 
d^^^N - d^y^M =^ -3y; 

therefore the factor is , . 



/ ' 24. In the case of homogeneous equations, a factor 
proper to render them integrable, is readily discovered by 
means of the property that if t^ be a homogeneous function of 
^ dimensions of the independent quantities t and %y then 

WW = tdiU + xdgU. 

For suppose F, a homogeneous function of w and j^ of m 
dimensions, to be a factor which makes M + Nd^y an exact 
differential coefficient, M and N being homogeneous functions 
of w and y of r dimensions ; then if U denote the primitive, 
it will be homogeneous and of m + r + 1 dimensions, and we 
/shall have 

VM-\-VNd,y^d^L\ d^^^U^VM, d^y^U^VN, 

«FJIf + yFJV - (»i + r + 1) !;■; 

M + Nd,y 1 d,U 

JIf a? + Ny w + r + 1 ' f7 
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and as the seoond member is an exact diiFerential coefficient, 
it follows that the first is so likewise, and consequently, that 
M + Nd^y is made exact by means of the multiplier 

1 

Mm + Ny " 



/. 



/^/ 26. The property of homogeneous functions assumed 
above is easily proved. Let t^ be a homogeneous function of 
the independent quantities ^ and x of n dimensions ; then if we 
change t into ^ (l + A) and « into iir (l + A), u will become 

w (1 + A)* = tt + nuh + &c. 

But by Taylor^s theorem, u will also become 

u + d^u . ht + dgU . hx + &c. 
therefore, equating the coefficients of h, 

nu = tdiU + zdgU, 

I' And, generally, if t^ be a homogeneous function of n di- 
mensions of any number of independent quantities tj z^ fr, &c., 
and we change them into ^ (l + A), ij? (l + A), u; (l + A), &c., 
the new value of u will be equally expressed by w (l + A)* or 
by c(**^ + *^'^' + '")'*tt; and equating the coefficients of A** in 
these two identical expressions, we get, separating as above 
the symbols of operation from those of quantity, 

»(n-l)...(w-r + 1) w = (tdt + xd^'\' wd^ + ...yu. 

Ex. fl?y + y^ + {my - ^) d^y = 0. 

The factor is 

1 1 



..>^«.«-i^i<i«i— >. ■ ■ « I j_^i>^«>. 



{(Joy '{''^)a) ^{wy '-aF)y ^fm^ 

is an exact difierential coefficient, and gives the primitive by 
Art. 7. 

^ + ilog(^y) + C = o. 



S5 



■i 26. Whenever the variables can be separated in an equa- 
tion, a factor which makes it integrable can also be found. 

For suppose that M + Nd^y = 0, by the introduction of 
two other variables u and x, is transformed into R + Sd^u b o, 
so that 

M + Ndgy - -B + SdgU ; 

and suppose T to be a function of u and z^ such that if 
we divide R + Sd^u by it, the variables are separated, i. e. 

— contains x only, and — contains u only; 

1 R Sf 

is an exact difiPerential coefficient; and consequently — , which, 

upon restoring the values of u and x, becomes a function of 
w and ^, is a factor which makes M + Nd^y = integrable. 

Ex. 1. a + ba^^ + a^d^y « 0. 



Assuming y « — > we find 



a + ba^j^ + a^d^y = a + 6w' + wd,u - «, 
and dividing by w(a-u + bu^) we get 



d^u 



a + ba^^ + a^d^y ^ 1 



a? (a - w + bu^) aw - a^^y + 6^y* 
is a factor which makes the proposed equation integrable. 

Ex. 2. y* + «a? + (l - a7y) d,y =» 0. 

4 
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By assuming y = , it may be shewn that a factor 

which makes the proposed integrable is 

1 
y* + Sayw — a + a*a^ * 

Ex. 3. M + Nd^y » 0, a homogeneous equation. 

In this case we know, that making y b wXj we have 

M = arf{x), N « af'<p («), 
and M + Nd^y = arf{%) + /i?''^ («) {z + wd^x] ; 

« 

consequently, dividing by a^'^* {/W + »0 («)}= Jfa? + JNTy, 
we get 

M + iVd,y 1 («) d,« 

Mw + Ny w fix)-\-x(f){x)^ 

1 

ilf <r + Ny 

is a factor which makes the proposed equation integrable. 

Obs. That -TT r^ is an exact difFerential coe£Scient, 

MW'\' Ny 

provided M and N be homogeneous functions of w and y of 
the same dimensions, admits of an easy proof as follows. 

We must shew that 'i(«(^^) - ^i'>{j^^} 

Now putting - = «, we have 

w 

M 1111 



Mw + Ny w Ny x l + F{z) 
N 1 / if ^ 



Mw-^ Ny 



1 / ^ if ^ \ 1 1 1 

y V " if ^ + JV^y j "y'"yl + F(«)' 



"■"— ^•-*-~^''^ — ^^— ^^^i^-^. . ^ ,. 
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d,F(x).- 



^'^\lifx + Ny)~ w ■ {\ + Fix)Y' 

^'Amw + Ny)" y' {l^F{x)Y' 
which expressions are evidently equal to one another. 

27. There always exists an infinite number of factors 
which render an equation of the first order and degree inte- 
grable. 

For let % be the factor by means of which the equation 

M + Nd^y = 
is made integrable, and u^O its complete primitive, so that 

« (M + Nd^y) = dgU ; 

multiply both sides by F(u)y where F{u) denotes any func- 
tion of u^ and we have 

z F(u) (M + ATd^y) = F (u) d^u ; 

and as the second member is an exact differential coefficient, 
it follows that the first is so likewise; therefore xF{u) is 
a factor which makes the proposed equation integrable, what- 
ever form be assigned to F(u). 

/^ 28. The following geometrical problems are added to 
illustrate this part of the subject. 

I. To determine the trajectory of a given family of 
curves. 

Let AA\ BB>^ (fig. 1.) be two of a family of curves 
resulting from the equation / {X^ F, c) = 0, by giving par^^ 
ticular values to the constant c ; and let AB be a curve which 



^ 



cuts AA^ BB^j and all the curves resulting from the equation 
by giving all passible values to c, at the same angle; thai 
AS is called the trajectory of this family of curves. Let 
m^ iff be the co-ordinates of the point A in ABf between which 
we are required to find a relation ; AT', AT, tangents to the 
curve and trajectory at A^ tan TAT* ^ a; and let a value 
^ {Xj F) of ef^Fbe obtained from the equaticxi/(X, F, e) » 0, 
not involving c; then at the point A, tan ATN ^ yff {x, y), 
and im ATN^d,y, 

. ^^ d^y-^i/e, y) 

the differential equation to the curve AB ; and as it does not 
involve e, AB will cut every curve in the series at an angle 
whose tangent » a. The equation when integrated will in- 
volve an arbitrary constant, and consequently will represent 
a tystem of curves, every one of which cuts the former system 
at the same angle; the constant may be determined, if a 
point through which the trajectory is to pass, be given. If 
the angle TAT' be a right angle, or (a) infinite, the differ- 
ential equation to the trajectory, which is then called ortho- 
gonal, is 

^'\rd^y.y\f(a,y)^0. 

Ex. 1. To find the orthogonal trajectory to all curves 
resulting from the equation y* (e — j?) = ^, by giving all 
possible values to c. 

Here c - a? « ~ , .-. -!■:-_ j d*y ; 

y* it y" 

if + Safy 
.-. d,y^y\f(jiOyy)^—^ — , 

therefore, substituting for y\f (w, y) its value, the differential 
equation to the trajectory is 

2a;* + (y* + 8a^y) d,y = ; 
A homogeneous equation whose integral is 

of + y^ MS C y/^a^ + y*. 



^t^mmm^^f^l^^ 
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Similarly, let /(p, 0, c) « be the equation to the curve 
A A' referred to polar co-ordinates, and let it give for pd^O 
the value yfr (p, 6) independent of c ; then considering p and 
as co-ordinates of the point A in the curve AB, 

tan SAT = ^ (p, 0), tan SAT ^ pd^9^ 



a s 






which is the differential equation to the trajectory; or if it 
be orthogonal, 

1 +ylf(pj 0)pdp0^O. 

Ex. 2. Let the curves be a system of circles touching a 
straight line in the same point, then taking that point as the 
origin and measuring from the line, their equation is 

p » c sin 0, 

1 cos , ^ ^ _ sin . , ^v 

= -7-3dp0, or pd^0 5«>K/>, 0); 

smfl '^ '^ *^ COS0 ^ ^ 



.-. 1 + 



P 
sin 



cos 



-pdp0 = O, or cos0 + sin0pdp0»O, 



or do ( — ^] - ; .'. p « C cos 0, 
'^ \cos 0J '^ 

the equation to the orthogonal trajectory, which represents 
a system of circles passing through the given point and having 
the given line for their diameter. 

We may generalize this problem, by finding the orthogonal 
trajectory of all circles described through two given points. 

II. To determine a curve such, that the locus of the 
extremity of its polar subtangent shall be a straight line. 

The polar subtangent is a line drawn from the origin 
perpendicular to the radius vector to meet the tangent. 

Let p, 0, be the polar co-ordinates of any point P in 
the curve sought (fig. 2) ; then those of the extremity T of 
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its polar subtangent will be p*d.9 and 6 , which must 

satisfy the equation to a straight line, via. 

p' = c sec (d' - a) ; 

.-. p^d^O = csec(0-.a-^) « -r-p^ r • 

'^ '^ \ 2/ sin(0-a) 

c c 

.'. -5«sin(0-a)dp0, and - = co8(0 -a) + C; 

P P 



^"^ C+cos(0-a)' 

the equation to curves of the second order, having the pole 
for one of their foci. 

III. To find a curve in which SG varies as SP^ PG 
being a normal at P, and SG a fixed line through S. (fig. 2). 

Taking SG for the axis of ^r, the equation to the normal 
at P is 

therefore making 

FcO, X^SG^w-fryd^yj 

and SP ■■ \/^"+y*, .'. « + yd^y = e \/a?*+ y*, 

or \/a7^+ y*ac^+ C, 
the equation to a curve of the second order. 

^ IV. To find a curve which is always cut by its radius 
vector at an angle proportional to the corresponding angle of 
revolution; that is, z SPT oc z ASP9 (fig. 2). 

Let py 0, be the co-ordinates of any point in the curve, 
then the angle at which the radius vector cuts the curve, has 
for its tangent pdp0; 



.•. pd^O^taxinOf or 



dep cos n9 
o sin n9 ' 



.*. I^j =mnn9. 



SECTION II. 

DIFFERENTIAL EQUATIOKS OF THE FIBST ORDER, BUT KOT 

OF THE FIBST DEGREE. 



/^ 29. When a difierential equation of the first order is of 
a bigh^ degree than the first, we know that it is not obtained 
by the* direct di£Perentiation of its primitive, but results from 
eliminating a constant, (which enters into the primitive in a 
dimension above the first,) between the primitive and its de- 
rived equation ; the degree of the difPerential equation and the 
dimension of the constant eliminated above the lowest dimen- 
sion in which it appears, being always the same. The general 
form of such equations free from radicals, is 

{d^yY + pi (d,y)"-' + p^ (d^yy-'' + ... + Pn^id»y + ft = 0, 
the coefiicients being functions of w and y. 

If this can be resolved with respect to d^y into its (n) 
simple factors, it will assume the form 

(d*y + gi) (d,y + ^2) — (d»y + gj = ; 

then each of these factors put equal to zero, will be an equation 
of the first onler and degree, whose integral may be found by 
the methods of the preceding section ; and any one of these 
integrals, as well as the continued product of any number of 
them, will evidently satisfy the proposed equation. If, there- 
fore, we integrate the n equations, 

d^y + gi = 0, d^y + ga = 0, ... d^y + g^ » 0, 

and OHnplete them all with the same constant C, as the 
proposed equation is of the first order, we shall obtain the 
required primitive involving C in the »*^ power, by equating 
their continued product to zero. 
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207 
Ex. 1. {d^yy + — d,y - I = 0; 

y 



w / a^ 

y y^ 

or -■ ^ = ± 1 ; 

V^ + o?* 

.'. + -x/y* + /»* = d? + C, and - \/y* + o?* = a? + C ; 

.-. (v^»T^« C - a?) (v^y* + a^ + C + «?) = 0, 

or y* + a;* - (C + o;)' = 0, 

or y* = 2Ca7 + C®. 

Ex. 2. 07* (d,y)* - 2o7y d,y + y* - 07*y* - ot* = o. 
.'. 07d,y - y = ± 07 VVT^, 



^■0 



or . = 1. 



^M0 

.'. V 07* + y* + y = co7e*, changing the constant, 
and v^oT* + y* - y = cxef ; 
.'. (\/^ + y* - co7e' + y) (\/o?^ + y* - co^c* - y) = 0, 



or 



07/1 \ 

y SB - I ce* I . 

^ 2 Vc^ / 
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30. When the resolution of the proposed equation into 
its simple factors is impossible, there are still various forms 
for which the complete primitive can be determined, or its 
determination made to depend on elimination; this is done 
by means of substitution, or differentiation, or other analytical 
artifices, of which we shall now give some instances. 

Obs. The arbitrary constant in what follows is often 
reserved under the sign of integration. 

/^ 31. When the equation contains only one of the variables, 
^ 0? suppose, and can be solved with respect to that variable, 

so that ^=/(d,y); let d^y be denoted by p, then ^=f(p) ; 

and integrating the equation d^y = p by parts, we get 

y = a?p - f^wd^p = pf(p) - f^f(p) ; 

between which and the equation w ^f(p), eliminating p, we 
shall obtain the required integral. 

Similarly, if we have y^fip)^ since 

dpW = dyW . d^y = -dpf{p), 

we shall have to eliminate p between y-fip)^ 

and a?«= / -dpf(p), 
JpP 

Ex. 1. a? + ^{d^yy ^ 1, 

1 



W =s- J, 



l+P' 



-tan"^p + C; 



.•. y as \/x (1 '-•07) - tan~^ \/ — ■- — + C. 

Ex.2. y=ay^i +p2, J? + C = a log {y/f - a* + y). 

This is the solution of the problem in which it is required 
to find a curve such that the perpendicular on the tangent 
from the foot of the ordinate shall be constant. 
6 
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Ex. S. y'y/\-¥^^apt 



r a a ( P \ 



or J^ 

This is the solution of the problem to find a curve such 
that the tangent terminated at the axis of w shall be of a 
constant length. 

1/ 32. An equation not coming immediately under this 
case, may sometimes be reduced to it by putting p = wz^ or 

P^yx. 

Ex. {d»yy + €iwdgy + ^ = 0. Let p = wz^ 
then 0? (ar 4- 1) + asf = 0, or a? = - 



••. d*y « d,y . d,a? = - 






(ijf« + ly ' 

and iv must be eliminated between the integral of this, and the 

az 
equation w = -5 . 

1 -^ ST 

P 33. When the equation contains both the variables ^ 
and 9, provided it be homogeneous with respect to them, we 

may assume — ^m z; then the equation will take the form, 

(which is not solvable with respect to p by supposition,) 

/(«,p) = 0. 

Suppose this capable of being solved with respect to z, and 
let it give iv » (p) ; now y^wz gives p^z -^ xd^z^ 

1 d^z 
or - = 



w p^ z 

substitute the above value of z^ and integrate this equation ; 
then p must be eliminated between the result, which will be 
of the form log w = F{p)^ and y ^ w(^ (p). 



>*«ateiwt«iki 
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Ex. y - wd^y = wa? v^l-h (d^yf ; 

y J — 

.-. - = p + »v 1 +p^; 
.-. log 0? = |log (p + vTTp^) + w log \/l + p*} + logC; 

re 

Vi + p* ^ 

between which equations p must be eliminated. 

« 

This is the solution of the problem, to find a curve such 
that the perpendicular upon the tangent from the origin 
shall vary as the abscissa to the point of contact. 



Let«=i. — = 1 +|?2 +p\/l +p*, - = p + \/l +p*; 






or y* = 2 C« — ^. 

34. Another integrable form is y « ^^p^ -^fid^y^ which 
is called Clairaut^s form, after the Mathematician who first 
considered it. Substituting p for d^y^ and differentiating^ 
we get 

y^wp +f(p)y 
d»y^p^p ^- OBd^p + dpf(p) . d^p ; 

which resolves itself into the two 

a? + dpf(p) = 0, d,p = 0. 



3G 

The first of these gives p ^ t^ {ai) suppose ; this Talue 
substituted for p in the proposed equation, furnishes a relation 
between x and y which satisfies the proposed equation, but 
which involves no arbitrary constant, and cannot therefore 
be the complete primitive. The other equation must there- 
fore lead to the complete primitive ; but this gives p = C, 
and by substituting this value of p in the proposed we find 

y = Ca?+/(C). 

Hence Clairaut'^s form has the property, that the complete 
primitive is obtained by substituting the arbitrary constant C 
for p, in that form. If we integrate p=C, we find y^Cos-^C\ 
but the condition of the proposed equation being satisfied gives 
C = f(C)f the same result as before. 

We shall after\v«irds return to the consideration of the 
other solution, which is called the singular solution, and is 
not derivable from the complete integral. 

a (1 + 1?') 
Ex. 1. y = «»p + — ^^ ^-^. 

P 
Difierentiating, we get /? = p + ofdj) + a 1 1 — ^ j d,p, 

a (1 + (?) 
,\ d^p = gives p^ C, and y = Co? + — ^ 



or 



the complete integral, and 



a 
w •\- a i = gives 

P ' ' x-^-a 



p^ivczr, 



which, substituted in |?y = b + p^ (^ + a), gives 

^y\/ = 2a,. or 2^ = 4a (a? + a),. 



OB -\- a 
the singular solution. 
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Ex. S. y = <i?p + \/l? + ay, 

y ^Cx-^^/h^ •\- C^d^^ the complete integral. 

a*y* + 6W = a^l? the singular solution. 

This is the solution of the problem to find a curve, such 
that the product of the perpendiculars dropped from two 
given points upon the tangent may be invariable; for taking 
the line joining the two given points (whose distance 
suppose s 2 c) for the axis of w^ and their middle point for 
origin, and w^ y the co-ordinates of any point in the curve, 
the equation to the tangent at that point will be 

Y-y^d^y(2r-'W)t or F = d,y^+ (y-a?d,y), 

and the lengths of the perpendiculars dropped upon this 
line from the points (c, o), (- c, o) will be 

\ATp^ ' a/i + P^ 

and the product of these is 

(y - ^pY - c'p^ -2 «„«^«^. 
—. — 2 = o , suppose ; 

.-. y ss ap + vft^ + a^p% putting a^ ^b^ -{■ c'. 

C-^ 36. A still more general case is the equation 

y = a?/(p) + 0(jp), 

which by differentiation is reduced to a linear equation of 
the iirst order in a,, for we get 

p =f(p) + ^^pfip) ^»P + ^p^ ip) d»p ; 
••• \p -f{p)]<ip'" = <«dpf{p) +rfp^(p) ; 

which gives x ^ F (p) ; then p must be eliminated between 
this and the proposed equation. 



i^^'^^//A/-A /f>/-^ 
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Ex. 1. y » xp* + 2p, 

r 

or d„a7 + 



^ p - 1 P^ ^P^ 



which is made integrable by the factor (p- 1)*; 

.-. jr (p - 1)* x= - fJLlZ « - 2p + logp* + C. 

•> P 

But p = ± \/ - + -^ ; therefore substituting this 

W W Or 

in the preceding, we obtain the complete primitive between 
X and y. 

Ex. 2. y = d? wp + n v^l +p' ; 

f» ,••_*»» m— 1 

^jiTTi^ » r _P 



Ex. 3. y + p(a-a?)*nji \/l +p*; 

^ 2c»(w + I) 2 (7> - 1) (o - 0?)""^ 

This is the solution of the problem of finding the path of 
a point P which moves uniformly towards another point Q, 
also moving uniformly in a straight line. 

For taking A (fig. 3.) for the origiq, and AB, which is 
perpendicular to By the line in which Q moves, for the axis 
of w, we have^ supposing P and Q to start together from A 
and fi, BQ = nAP, or if AN » ^, NP = y, ^J? « a, 

y + (o - J?) p = n^v^l +p'. 



m "J ■■ ^m^^^f^^^^mmm^ 
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36. In the following examples the method of substi- 
tution succeeds. 

Ex. 1. (l-p*)^y«=p(ir'*-y*-c*), 

which expresses that the normal bisects the angle between the 
focal distances; 2c being the distance of the foci, the origin at 
the middle point between them, and the line joining them the 
axis of d7. 

wx % 

Let p = — , ,\ y* s= a^x - c* 



y 1 +« 

therefore, differentiating, 



( (1 + xyf ' ' 



this resolves itself into 

^ a- ± ^ which gives y* 4- (^ — c)* = 0, 

the singular solution; 

and dgX = 0, or «? = C, which gives y* = C { ** — ) , 

the complete integral. 

If we integrate p = — , we get y^ ^ Cjb^ + C , where C' 

if 

must be determined by the condition of the proposed equation 
being satisfied ; and by this condition, in general whenever the 
method of solution raises the order of the equation, must the 
number of constants be reduced. 

By the same substitution may be solved the more general 
form 

aayp* + p (/p* - ay* - 6) - /ry « o. 

Ex. 2. ^d,y-y = Jr \/(4y)«-^d,y+l, 
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y 

Let y = wx^ then d^y — x •\- xd^x = - + oid^x ; 

dgX X 



If -y= 1, we have sin"*»a= sec"**!? + C, 

or sin"* - « sec"*^ + C, 

X 

Ex. 3. and 4. 



VTT^ 



=/(a/^ + s^), 



y-xp f( /" V 

Introducing polar co-ordinates, we get for the first, 

/.^ ^' , =/(/>). or d,0 fM==^ 

Videpf + p* '^ pvp'-lfip)}* 

The second gives 

P fl^e) or ^ = Vl-{/(cosg)}' 

\/(rfep)' + p" •'' " p /(COS0) 

The former expresses that the perpendicular on the tangent 
from the origin is a given function of the radius vector ; the 
latter, that the sine of the angle at which the radius vector 
cuts the curve is a given function of the cosine of the angle 
at which it is inclined to the axis of <r. 



SECTION III. 

ON THE SINGULAR SOLUTIONS OF 1>IFFEEENTIAL 

EQUATIONS. 



//' 37. From the complete integral of a differential equation 
'^we can deduce as many particular integrals as we please, by 
giving to the arbitrary constant particular values. But some 
differential equations are satisfied by a relation between w and 
y not containing an arbitrary constant, and not deducible from 
the complete integral ; such a relation is called, as has been 
said, a singular solution of the differential equation. The 
existence of such solutions depends upon the fact, that when 
a solution of a differential equation has been obtained in any 
manner, it will still be a solution after a quantity of any kind 
has been introduced in any way, provided the same derived 
equation result. This is merely an extension of the principle 
on which the arbitrary constant is added. 

y-r^^ 38. Before entering upon the general theory, it may be 
useful to consider the following particular instance. 

Let the equation 

y = wd^y -^-a^-a {d^yf (l) 

be proposed, which, since it falls under Clairaut^s form, has 
for its complete integral 

y^CW'{-a + aC* ; (2) 

C being the arbitrary constant. If we now regard C, not as 
a constant, but as a function of a?, and differentiate, we get 

and if we eliminate C between this and y = Co? + a + aC*, we 
shall obtain a differential equation, but not the proposed one, 
6 
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Tor that arises by eliminating C by means of the equation 
d^y « C But if C be so determined as to make the co- 

efficient of d,C vanish, that is, if C «= , then the derived 

2a 

"equation will be d^y « C, and the result of the elimination 
of C will be the proposed equation. 

Substituting for C its value, we get 

a^ or of 

«« +a + — = ho, 

2a 4a 4a 

a relation which manifestly satisfies the proposed equation; 

for it gives d^y s , and these values of y and d^y^ being 

2a 

substituted in the proposed equation, make it identical. But 
this solution contains no arbitrary constant, and yet being the 
equation to a parabola it cannot, either by making C = 0, 
ar any other constant quantity, arise from the complete in- 
tegral which is the equation to a straight line ; it is conse- 
i|uently a singular solution, and arises from the complete 
integral by changing C into a function of so so determined as 
to make the term involving d^C disappear from the value 
of d^y. 

Thus we see how the singular solution arises from the. 
complete integral ; next, let us consider its geometrical signi- 
fication. The proposed differential equation expresses that 
the curves to which it belongs have the property, that the 
tangent at any point is intersected by a perpendicular upon 
it from a given point in a given straight line. 

Take the given point S (fig. 4.) for the origin, and TSy 
AS respectively parallel and perpendicular to the given line ACy 
for the axes of oo and y. Let TC be the tangent at a point 
whose co-ordinates are x and y, then its equation is 

«nd the equation to the perpendicular upon it from S is 
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and for their point of intersection 

but this point is always in AC for which Y «= a; 

the same as the proposed equation. 

The complete integral y = Cof + a (1 + C*), which repre- 
sents a series of straight lines, evidently satisfies the problem 
for all values of C; for let Ti be any one of these lines, 

then y a» — - is the equation to a line through S perpen- 

dicular to it; and combining the equations to get the co- 
ordinates of their point of intersection, we have 

y = - C*y + o(l + C*), or y = a; 

the intersection consequently falls in AC. And as a straight 
line is its own tangent at every point, the equation 

y = Co? + a (1 + C*), for all values of C 

represents a line such that the intersection of the tangent 
at any point, and a perpendicular upon it from S^ falls in 
the given line AC. Now the curve which is generated by 
the perpetual intersections of these lines will also satisfy the 
problem; for each of the lines will be a tangent to it, and 
therefore perpendiculars from S upon its tangents will inter- 
sect them in the various points of AC. To get the equation 
to this curve we must, according to the usual method, dif- 
ferentiate with respect to the parameter C, which gives 
= 07 + SCa, and eliminate C between this, and the equatioa 
y =r Cw + a (1 + C^), which gives 

y = + a, or 4a (a - y) « a?*, 

4a 

the equation to a parabola, vertex A^ focus Sj of which 
curve it is a well known property, that the perpendicular 
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from the focus intersects the tangent at any point, in the 
tangent at the vertex. 

This result being obtained by exactly the same process 
as the singular solution was obtained, of course coincides with 
it; hence it appears that the singular solution belongs to 
the curve which touches the family of curves resulting from 
the complete integral by making the arbitrary constant assume 
all possible values. The conclusions arrived at in this par- 
ticular instance, we shall now shew to hold generally. 

39. Having given the complete integral of a differential 
equation, to find its singular solution. 

Let /(/p, y, d^y) = be a proposed differential equation, 
and suppose it to result from the elimination of the arbitrary 
constant c, between the equation F (^, y^ c) «= 0, and its 
immediately derived equation M + Nd^y s o. Now change 
c into e', any function of w and y, then our equation becomes 
F(af9 y, e') «0; and its immediately derived equation 

Jir + N'd^y + Cd,c = (1), 

c' entering into JT and iVT just as c did into M and JV, 
and dgc' of course denoting d^^^c + d^^^c' . d,y ; now C is the 
differential coefficient of F (a?, y, c) with respect to c, regard- 
ing w and y as constant, and will therefore usually involve 
^, y, and c' \ and if put equal to zero, will give such a 
value for c as makes the last term of equation (l) disappear ; 
and then the elimination of c' must evidently produce the 
proposed equation /(zp, y, d^y) = 0. Let this value of c' be 
substituted in F (/r, y^ c') = ; then this equation is changed 
into (p (Wf y) « 0, and furnishes a relation between of and y 
which satisfies the equation /(/Vy y, d^y) » o, but contains no 
arbitr^ constant, and is not deducible from the complete 
integral by giving a particular value to the constant; since 
it results from the complete integral by substituting for c 
a variable value deduced from the equation d^ F (w^ y, c) » o. 
Consequently, the relation ^ (or, y) s o is the singular solution 
required. 
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40. To explain the geometrical signification of the sin- 
gular solution of a differential equation. 

Let F{po^ y, c) ^0 (1) be the complete integral of a dif- 
ferential equation between two variables; if we differentiate 
it with regard to c, we have d^Fiw, ^j c) = (2); and if 
between these equations we eliminate c, we get 0(a7, y)s=0 (3), 
where c does not appear, and which is a singular solution of the 
differential equation, of which (l) is the complete primitive. 
Suppose equation (l) to be the equation to a system of curves, 
in which the position and dimensions of any particular curve 
is defined by a particular value of the parameter c ; also let 
equation (3) be the equation to a curve referred to the same 
co-ordinate axes. Now equations (l) and (2), when e receives 
a certain value, are satisfied by the same values of w and y. 
Hence from the manner of its formation, equation (3) is satis- 
fied by the same values ; or the curves which are represented 
by (3) and (l), with a particular value of c, have a common 
point. But equations (3) and (l), being each a solution of 
the same differential equation, furnish the same value of d^y 
for the same values of a? and y; consequently the curves 
touch one another at their common point. The same thing 
happens for every one of the system of curves which equation 
(l) represents. Therefore the curve represented by the sin- 
gular solution touches in a point every curve represented by 
the complete primitive. This is the geometrical interpret- 
ation of the singular solution of a differential equation of the 
first order. 

41. Having given a solution of a differential equation, to 
find whether it is included in the complete integral or not. 

Let d,y =»/(^> y) be the proposed differential equation^ 
and y ^ F{w^ c) its complete integral, e being the arbitrary 
constant ; and . when c s c\ let this become y^u^u contain- 
ing no arbitrary constant ; then y ^u is a particular integral 
of the proposed. 

Hence, since ^ (a?, c) — w becomes zero when c « c , we 
have 

F(^, c) ':^ u:=i{c - c'y^. as = a%9 suppose, 



•■Sfc;'*, 
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^ being a function of x and c, which is neither infinite nor 
zero, when c^c\ or when a = 0; and m expressing the highest 
power of c - c which enters into every term of F (w, c) - u. 

Consequently the complete integral becomes 

■ 
y = w + aZf 

which being substituted in the proposed equation, 

gives dgU + adgZ =/(a?, u + ax) (l). 

Now since z is neither infinite nor zero when a = 0, we 
may expand it in a series of ascending powers of a, in the 
form 

)!? = JT + Ja« + So'' + &c., 

a, /3, &c. being increasing and positive, and iT, J, jS, Sec 
functions of ^; 

.*. dgU + adgX = d,w + d,Ka + d,-4a""*"^ + &c. 

Again, the development of /(a?, u + az) will be of the 
form 

f((Vj u + az) =/(a?, w) + Jif (a;^)** + N(azy + &c., 

»»» ^9 &c. being increasing and positive. Hence, by substi- 
tution in equation (l), observing that d«f^=/(<r, u)y we get 

d,jr . a + dgA . a^"^^ + &c. 
= if a~ (JT + ui a« + &c.)" + JV^a* (JT + Ja« + &c.)" + &c. 

Now unless this equation is identical, y = u cannot result 
^ s p/l, Ip y from ififel^M^ l>y changing c into c' ; and •*. y s t^ cannot 

be a paiticular integral of the proposed; and if it satisfies 
the proposed, it must be a singular solution. Now the indices 
^9 Tiy &c., are known, for they result from writing u + az 
for y in /(w, y), and expanding according to power of az ; 
and we must endeavour to determine a, /3, &c. so that the 
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two members of the equation may be identical. If m be > 1, 
this can easily be effected ; for we must have d, JT = 0, or 
K = a constant ; a+1 = niy and d,A =i MK^ ; and so on 
for the other terms. Consequently it will be possible to make 
the two members identical, and y = u will be a particular in- 
tegral. In the same way the identity may be established if 
m« 1. But if m< 1, there is no term on the first side cor- 
responding to MK^a^\ and since iT cannot be equal to zero, 
it is impossible to satisfy the identity ; and therefore y «= «^ is a 
singular solution. Hence to discover whether a given solution, 
y -Uf of a differential equation d^y =fi/xfy y), is a singular 
solution or not ; we must write u^h for y in the value of 
d^y^ and if the expansion in ascending powers of h involve 
a power of A, whose index is < 1, the solution in question is 
a singular solution. 

42. To deduce the singular solutions from the differential 
equation, without knowing its complete primitive. 

Let y ^u he B, singular solution of the equation 

then by the preceding article, substituting u •\' h for ^, we 
get 

f{w^U'\- h) =/(a?, u) + MK^ + NK" + &c. 

where th, n^ &c. are proper fractions ; 

•'• d^f(Wf w + A) = dj^f{Wi w + A) = mMhJ^"^ + nNK"'^ + &c.; 

consequently, when A « 0, d^f{pBj w) = 00* 

But d^f{w^ u) is what dyf{wj y) becomes when y = u; 
and therefore, conversely, every value, u of y, which satisfies 
d^y =/(a?, y), and makes d^f(wj y) =» 00, is a singular solu- 
tion of the equation d^y =s/(^» y). 



*-• •.* -i- 



43. It is not essential to give the equation the explicit 
form d^y^f(v^y). For let d^y-p, and let V=0 be the 
given relation between a?, y and p ; then we may regard p 
as a function of ^ and y determined by the equation r=0; 
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'• •^m'^+^(p)^<*vP = 0» ^^ ^vP"- 



'(jr) 



Hence the condition dyf{wjy)^a^ is equivalent to 
^f^^^FsO, provided d^^V remains finite; and therefore the 
singular solutions of the equation V^ F{xj y, p) « o, are 
determined by eliminating p between V^O and d^^^V — O^ 
provided always that these solutions do not make d^^V-0. 
It is evident that if we consider y as the independent vari- 
able, and put the equation under the form 

F= F (w, y, d^of) ^ Oy 

the same reasonings would shew that singular solutions may 
be obtained by eliminating p' « d^of^ between V — and 

rf'F-O. 

p 

Ex. 1. To find the singular solution of 

y - wd^y + ^ - j^ = «• 

Here rf^)F« -4? + — «0; .•. p* = -, and the proposed 

/^ 
becomes 

(tff + y-o)p = y + ^pS or (af + y-a)p^2y^ 
or (« + y - ay = 4a?y, 

which may be reduced to the form y/w + y/y = y/a. 

Ex. 2. To find the singular solution of 

(y - ofd^y) Iw - -^ J « a*. 4a?y = a^. 

\ dmyl 

Ex. 3« To find the singular solution of 



1 
I 
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These three examples determine respectively the curves 
which have the properties that OT^-OT is constant, that 
area of triangle 7^07^ is constant, and that TT^ is constant, 
TT^ bdng the tangent at any point meeting the axes 0<r, Oy 
in T and T. (fig. 5.) 

Ex. 4. y + Ja? - tan-^ (d,y)} d,y - i = 0. 



H«re d^^Y = « - tan"* p - ^ ^ « ; 






-1- . P 



a?=«tan"'p + 



l-^f 



y = ^-— ^> .*. a^-cos-^Vy + v/y-j^, 

which represents a cycloid whose base coincides with the axis 
of (By the origin being in the centre of the base. This is 
the solution of the problem to find a curve always touched 
by the same diameter of a circle rolling along a straight 
line. 

£x. 5. dgif 



y/a/^-^f-a^-y' 



.*. the relation cV^ + J^-a^-O makes d^^p infinite, and satisfies 
the proposed equation; it is consequently a singular solution 
of the proposed. 

Ex. 6. To find whether or not ^+j^-a'aO is a sin- 

Of 

gular solution of d^y « —^ = . The solution 

V «* + y* - «* - y 

gives y = v «* - «»*, therefore substituting \/a* - a;* + A for 
y in the value of d^y, we get 



when developed according to powers of A ; and as the index 
of A is a proper fraction, a^^ + j^-a' — oisa singular solution, 

7 
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44. Every factor pn^r to make a propoeed differential 
equation intqirable, is made infinite bj the singular solution. 

Let d^y +/(^9 y) » be the proposed equation, F{ay y) — c 
its complete integral, and z the factor which makes it in- 
tegrable, so that 

« {d^y + f{^y y)} =- d,F(wj y) ; 

also, let y^u he the singular solution; then this is not 
dedudble from the complete integral, and therefore if t« be 
written for y in F{wj y)^ the result will not be constant; if 
therefore we substitute u for y, in the preceding equation, 
since the second member will have a finite value, and the 
factor d,y + /(<r, y) of the first member will be zero, the 
value of X corresponding to this substitution must be in- 
finite. 

This property will sometimes lead to the discovery of the 
&ctor which makes an equation integrable ; as in the example 



(of - a^) d^y ■¥ wy^a ^cf + y* - «*> 

a singular solution of which is a;' + ^ — a^ = ; if we try a 
factor of the form 

we arrive at m«s - i, n «ft -• ^; and the factor whidi makes 
the proposed integrable is 



♦ . 
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SECTION IV. 

DIFFERENTIAL EaUATlONS OF THE SECOND ORDEft, AW OF 

HIGHER ORD]CRS. 



45. Every differentia! equation of the n^ order admits 
of a primitive with n arbitrary constants. 

Let /(^9 y9 C], Cs, ... c«) ^ o be an equation between the 
variables of and y^ containing n constants C|, Cs, •.. c.. Let 
the first n derived equaticms be 

/«(^j y, d,y, djy, c,, ... c,) = 0, 



/.(^» y» d,y, d^y, ... c^y, Ci, ... c„) = 0. 

Between these n equations and the original, the n con- 
stants may be eliminated, and the result will be 

F(a, y, d,y, djy, ... d^y) = 0, (l) 

a differential equation in which none of the constants enter. 

Conversely, a differential equation of the n^ order being 
proposed, it must admit of a primitive containing n arbitrary 
constants, because this number of constants, and no more, can 
be eliminated in its formation. Hence every differential equa- 
tion of the n}^ order admits of a primitive containing n arbitrary 
constants* 

46. Again, between the original equation and its first 
n - 1 derived equations, n - 1 of the constants may be eli- 
minated, and a differential equation of the (n - l)^ order with 
one constant will result* 
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Every gucb differential equation, having the same pri- 
mitive with equation (l), is a first integral of that equation ; 
hence a differential equation of the n^^ order has n first in- 
tegrals, each a differential equation of the (n - 1)^ order, 
and containing one constant. 

Also, between the original equation and the first (r) of its 
derived equations, (r) of the constants may be eliminated, and 
a differential equation of the r^^ order containing n-^r con- 
stants, will result, which is an integral of equation (1). Now 
(r) constants can be eliminated in a number of ways equal to 
the number of combinations of n things taken r together, or 

n (n — 1) ,.. (n - r + 1) 
1 .2.3...r 

This then is the number of integrals of equation (1) of 
the (n - r)^ order, each a differential equation of the r^ 
order, and containing n — r constants. 

47* Of the general equation of the second order 

F{dly, d,y, y, ^) = 0, 

we shall first of all consider the following particular cases, in 
which dlff is involved with only one, or two, of the other 
quantities or, y, d,y ; and which admit of integration, or rather 
of reduction to forms of the first order. 

I. F (dlyy a?) = 0. Let this by resolution, give 

and integrating again, and adding another constant, we obtain 
the complete integral. The same process applies to d!ly^f{pB). 
Also, if we have F(d]J"*'^y, djjy) = 0, and put d^y^Uj we 
get F (dgUf^ u)-0; and if this can be integrated, and gives 
u Bs f(ai)j it is reduced to the case just noticed. 

Ex. a^dly « o, •••<*!» = :i > rf#» = + C ; 

or X 

.'. y * o log - + C^ + C". 
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II. F(tllyy d,y) mO. Let this by resolution, give 
djy - /(d,y), or d,p =. f(p), putting d,y=pi 



J 1 r ^ 

and d,y^d^.d,a,^^, .-. y-jT^i 

p must be eliminated between these two equations. 

p w 

Ex. ad!y = d^y ; ••• d»p = - ; •'• log » « - + C, 
' a a 

and dpV^dgy.dpa^p.-^ a; .*. y ^ ap -^ C ; 

P 

.-. - + C = log ^^ . 

a ° a 

III. FCe^y, y) = O. Let this give <^y =/(y); 
.-. 2d,y<Cy « 2/(y) <*,». and (rf.yV - C + 2 J^/Cy) ; 

••• <*»» - VC + ifgf(y), and ar - f /- 

•'yVC + 2 



2/y/(y) 

Also F(dJ+*y, djjy) =0, putting d^y ^ u^ becomes 

F(dlu, w) = 0; 

and this, treated as above, gives co ^ Kp (u) ; and if <r = (t«) 
can be solved with respect to u^ it is brought under Case I. 

Ex. 1. a^dly + y = 0; .'. ^a^d^ydly + ^yd^y = 0; 
.-. «"(d,y)* + !^ = C; .-. -^L==i; 

/. a sin"* — y= « i!P + C. 
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Ex. 2. ^/ay dly = i. 

IV. /•(<?», d[,y, «)= 0. 

P for rf,y and d,p for <y ; let its integral be ^ (,; pf Q . o 
.J^ *"'*u^ "solution give p or d[.y =/(^), then y =//(*); 

2.^41 f 5^ ^ ^'''^- Q ' for it «» be solved « a linear 
equation of the first order. 

If it gives a =/(p), then 

y = f,p = wp- f^wd^p ^aip- f^f(p) ; 
and p must be eliminated between these equations. 

Ex. 1. dlj, + -d,j,^o; .: xd^,y + d,y^d,ixd,y)^0; 

■■' a>d,y^C\ .-. y = Cloga; + C'. 
Ex. 2. (1 + a^) <5y + J + (d^j,)» ^ Q. 

l+p' iT^* ' t«n-'p + tan-'a? = tan-'c; 

.-. p^^^i(ij^_,\: 

l-¥cw c \i^cw /' 
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• A ^ ^ J 77:^"^' suppose; 

.•. p c J =, and y=/— F==. 

This is the solution of the inverse problem of the radius 
of curvature, in which it is required to find the curve whose 
radius of curvature is any given function of the abscissa. 

Ex. 4. djy + (e* - 1) d^y = c**, y « ^ + Ce'"^ + C. 

Putting d^y = p, we get c^y ^ d^p .d^y^ P^yP \ ^nd the 
substitution of pd^p for c^y, and of p for d,y, will make 
the proposed of the first order in p and y ; let its integral be 

If this by resolution give P "=" d^y = f (y)^ then 
If it gives y=/(y)j then 
and p must be eliminated between these equations. 

Ex... ,dA.(..»y.o, 

or d, (yd^y) = ; .-. yd^y = C; .•. y* = 2Ca? + C'. 
Ex. 2. • yc^y + n (d^yY + n « 0. 

•'• yP^yP + ^P* + n = ; 
.-. ^ log (1 + p*) + » log y = ^ log C ; 
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•• (1+P')»**«C; 



P 



^y/W^h and a^^j^j£^=. 



V^C-j^ 



This is the solution of the problem in which it is required 
to find a curve whose radius of curvature shall vary as its 
normal; for this condition gives 

— or + according as the curve is convex or concave to the 
axis of ^. If n s 1 the curve is a circle, if n « 12 a cycloid, 
if n B - 1 a common catenary. 

Ex.3. i + (d,y)*-2ffc5y«Say{l + (d,y)"}i 

ayf+ C 



this gives d x 



v/y-(oyi + Cy 



Ex. 4. S^dJy+\/l-h(d,y)*=0, gives dyX^-y==l===== . 

VI. F(dly,y,ai)^0. 

As there is no substitution by which this can be generally 
reduced to an equation of the first order between two variables, 
the artifice to be employed in any case will depend upon the 
nature of the example proposed. Among other substitutions 
for dP^y^ the two foUowing may be noticed,. 



w 



d5y = d*(^d*- -)» and xdly^^^{wy)^Zd^y. 



Ex. 1. OT^dJy = 2y; 



••• ^ ipoy) - 2d.y = — , 

at 
or ^{wy) = - d, {wy) ; 



«7 

•'• d, (a^jf) « Car* ; 

C^ c 

^ 3 J? 

Ex. 2. (^ + j^y djy + o*y = 0, 



or 



V w) (or -^ yy 



^^UJ 



V Of) \ '\w)\ L^/y\T 



1+ 



/ ^ 

1 + J8r 

1 r vTT«^ 

or --= / y 

Ex. 3. d|^ = aa7 + &y; 

or c^ (oa? H- 6y) = 6 (aa? + fey), 

which becomes d^^sfei^r, putting aw + by = %^ and so falls 
under Case III. 

48. When the equation 
^(^> y* d,y, d^y) « 0, or F(cr, y, p, d,p) = 0, 

is homogeneous, reckoning the dimensions of p and d^p to 
be and — 1 respectively ; it may be reduced to an equation 
of the first order by putting 



q 
y = wxj and d,p = - . 



8 
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For each term, if r denote its dimeDsions, will consist of 
some function of p multiplied by a factor of the form 



g)"^-(d,p)% 



m and n being any numbers from to oo ; therefore upon 
making the substitutions stated above, every term will be 
divisible by ar% and the equation will assume the form 

and if this can be solved relative to q, we shall have 

d p 
? « (^j J>) ; but q = wd^p ^(p- x)-—^ (p-x) d^p ; 

••• («> P) + (« - p) d,p = ; 
let this give 

p = >/f («), then \\f («)=« + a?d,«, 

which will give the required integral 

Ex. n{l + (d.y)»}» = dJ»x/^T7; 

q 

putting d^y = p, y ^ wx^ d,p = - , we find 



V 1 + « 
which is the same as Ex. 5. Art. 18 
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Linear Equations of the Second and Higher Orders. 
"[ 49. The linear equation of the w*^ order is» 

all the coefficients being functions of cp^ and each term of the 
first member involving either y<, or one of its differential co- 
efficients, in the first power. The first step towards the 
integration of this equation is the establishment of the following 
theorem. 

If there be n particular values Ui^ U2y e^ ... t«„» functions 
of <a7, which, when substituted for y^ satisfy the equation 

its complete integral is 

y «= OiUi + Ofith + <h^ + ..• + «««^»> 
^19 ^89 <»-:• ^n being arbitrary constants. 

For let this value of y be substituted in the expression 

d'ly + Pid!!"'y + ... +p«y9 

and it becomes 

or, collecting the terms multiplied by the factors ^i, a^, ... a^. 
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Now, since t«i, u^, ,.. u^ satisfy the equation, each of the 
quantities within brackets is equal to zero, and therefore the 
whole is identically zero ; and therefore the assumed value of 
y satisfies the equation ; and it contains n arbitrary constants, 
consequently it is the complete integral of the equation. 

Thus the equation cfjy - n-y = is satisfied by y =* 6"', 
and since it is not altered by changing the sign of n, it is 
also satisfied by y = «""'; .*. y = a,c^ + flgC""* is the com- 
plete integral. 

Linear Equations with constant Coefficients. 
50* To integrate the equation 

all the coefficients and u being constant. 

First write y + — instead of y, and the equation becomes 

Let y = c"', .-. c'"*(m" + p,wi*"* + ... + j»,) is the value 
of the first member; now this will vanish if m be any root 
of the equation 

m* + p,m"~* + pam""^ + •.. + p„ = 0, or Jf = /(m) = 0, 

which is called the auxiliary equation of the proposed linear 
equation. 

Hence the n real or imaginary roots of this equation, 
mi, f?i2» nh •«• ^n9 provided they be all unequal^ will give n 
different particular values of y, e*"!*, c"**, ... «*"•*, which satisfy 
the proposed equation ; and therefore its complete integral is 

y = OiC^i' + Og^** + ... + a^e^*'. 

51. But if any of the roots are equal to one another, 
as, for instance, m^^ m^^ the value of y becomes 

y = («! + Og) ^^' + Ogees' + ... a„e«»', 
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which contains only n — 1 arbitrary constants (because Oj + dg 
can be reckoned only as a single constant), and therefore 
cannot be the complete integral of the proposed equation. 
In this case, in order to discover the complete integral, first 
suppose the two roots m^, m^, to be only very nearly equal 
to one another, so that nig » f9»i + A, where A is a very small 
known quantity ; then the part of the value of y correspond- 
ing to these roots is 

^ c"!' (ci + CoW -h CgA — + &c.), 

replacing the constants a^ + ag, and OjA, by Ci, Cg, respectively. i 
Now let A = 0, then this becomes «*"»' (cj + Cg^) ; and the ' 
complete integral consequently is 

y = (c, + c^cd) c"i' + <h^^ + ••• + <?««"•*• 

62. Generally, if we suppose r roots of the auxiliary 
equation to be tiearly equal to one another, and therefore 
to be represented by 

i»j + Aj, fHi + As, ••• ii>i + Ar> where A^, Ag, •••A,. 

are very small quantities, the complete integral takes the 
form 

or, expanding c*i*, e***, &c., 

y = c"i'|2(a) + 2(aA).^p+-^^ — ^af*+ ... 

1.2 



r - 1 r ' 



Lr 



»«»a.i* 



+ ar+i« '+» +'-. 
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2(oA'-') 
or, replacing the constants 2(a), 2 (a A), ... — j , by 

1 J ^2 • • • Cf 9 

y = c*i'(c^ + CgO? + ... + Cya?'"* + terms involving Aj, A,, &c.) 

Now let A| = s Ag « ... = A,., in which case the auxiliary 
equation has r roots each « m| ; then the solution becomes 

which contains 9» arbitrary constants, and is consequently the 
complete integral. 

53. Of the correctness of the above modification for the 
case of equal roots, we may assure ourselves by the following 
reverse process. Let y = e^'w, then since 

d;(wt?) ^dj^v.u -^ndl'^v.d^gU + di'^v . dlu -^ &c. 

d;(e"'w) =s e^\m*. u + wi»""^ d, w + —^^ ^ m""*djw + feci 

separating the symbol of operation from that of quantity. 

Hence the first member of the equation becomes 
c"^(fw+d,)*w+Pie"''(m+d,)""*w+P2e~'(^+d,)""*t^+&c.+p.e"'w, 

= 6^{/(w)w+V(m)d,w+-^/V)«iJw+... + p/"VK«^}. 

1 1.2 yn 

Now let f7t s mj, then since /(m) = has r roots equal to 
in,, each of the quantities /(m), /'(w), .../^''"^^(m) becomes 
equal to zero, and the first r terms vanish; and if 
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then €^t^, dl'^^Uy ... d*u likewise become zero, and all the 
remaining terms vanish. Hence for each group of r equal 
roots in the auxiliary equation, there will be in the value of 
y a term of the form 



54. Let A±A;\/-1 be a pair of imaginary roots in 
the auxiliary equation, then the corresponding terms in the 
value of y will he 

or e^ \C (cos A?a? + v - 1 sin kw) + C' (cos kw - v^ sin kw) ] 
= e*'{ciCOS&^+CgSinA?a?|, (changing the arbitrary constants); 



or 



= \/c? + c?e^< — p=L=cos kw + . „^ sin A:a?> 



= )3e*'cos(fta? + a), 
again changing the constants by putting 

\/cJ + C2 = i3, and - tan a = — . 

And if there should be r pairs of imaginary roots in the 
auxiliary equation, equal to A ^ & y/- 1, the corresponding 
part of the value of y will be 

€*'(ai + a^of + ... ar^*""^) (coskw + \/- 1 sin A^/v) 

+ c^' (&i + AgO? + ... brOJ^'^) (cosk{9 - v- 1 sin A?d7), 

or, changing the arbitrary constants, 

^'(oi + a2«2? + ..• ttr^*""^) co^kw 
+ ^' (A + /32^ + ... i3,a?'-0 sin A;^. 
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65. We shall now give some examples of integ^ting 
linear ecjuations with constant coefficients. 

Ex. 1. djy + d,y- 2y = 0. 

Let y = e"', .-. fii*e"' -h me"' - 2c*' « 0, 
or i»^ + m - 2 =a 0; .-. m «= 1, or -2; 

Ex. 2. a^d*y = cgy, y = C,c* + C^e'^ + Cs^? + C*, 

Ex. 3. d^y -^n^y « 0. 

Let y = 6"'; .". w^ + »- = 0, or f?» » :^ n\/- 1; 

= Cj (cos nw + \/— 1 sin n^) 

■ 

+ Cg (cos no? - y/- 1 sin no?) 
= (Ci -h Cg) cosn^ + (Ci - Cg) \/- 1 sin nw 
8 a, cos no? + as sinn^, (changing the constants) 
s /3 (cosna? cosa - sinn^ sina) «= j3cos(n^ + a)« 

This equation, d^y -^ n^y »= 0, of which the solution is 
y = O] cosn<r + a^ sin nofy or y ^ fi cos (na + a), 
is of very frequent occurrence. Hence also the solution of 

Ex. 4. d^y + n^y + aa; + 6 = 0, or 

/ aw b\ ( aw b\ 

aw b 
n TIT 
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Ex. 5. tP^y- Smd^y + (m* + n*)y » 0; 

••• y '^^'(C cos »a? + C' sin «a?). 
Ex. 6. i^y + 2md«y + n^y = 0; 

y « tf"*""/3 cos (a? -x/w* - m* + a), (« > m). 

Ex. 7. diy - adjy 4- 2rf> - 2d,y + y « ; 
tb? equation for determining t» is 

(m« + 1) («t - 1)* - ; 
y = c* (Oj + a^w) + /3 cos (^ + a). 

Ex. 8. d*y + Sw'c^y H- n*y = ; 

y « (a + a^a?) cos no? + (6 + 6, a?) sin wa?. 

Ex. 9. diy + 4>mtP^y + 2 (»^ + Sw*) djy 

+ 4fm (i»* + «*) d,y + (•»* + »*)* y as ; 

the equation for determining k^ when y a e^', is 

.•. y » e-«* (a + a,^) cosno^ + c""' (6 + fc^a?) sin nw. 

Ex. 10. d*y - 4dly + 6djy - 4d,y 4- y « 0; 

y = ^ (^ 4. fia? H- C/P* H- Da?*). 

/^^ 66. To integrate th^ linear equation of the fi^ order 
with variable coefficients, 

d^y -*- Pidl~'y -^ Pidl'^'y ... ^p^y^X. 
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Let Ui, 1^9 •-• ^» be the n particular values of u which 
satisfy the equation 

dlu ■\- p^dl'^u + pidl'^u -¥ ... +p»ti = 0, (1) 

which coincides with the proposed if we suppose its second 
member to become zero; then, as already proved. 

If we now divide both sides by u^ and differentiate, we 
shall eliminate c^ ; next dividing both sides by the coefficient 
of Cj, which suppose V|, and differentiating, we shall eliminate 
c^; again dividing by v.^, the coefficient of 03, and differen- 
tiating, we shall eliminate 03 ; and proceeding in this manner 
till all the constants are eliminated, our final result will be 
of the form (where each dg affects the whole of the expression 
after it) 

djg dg •»" dg — dg — = 0, 

in which, the coefficient of dlu is evidently ; 

therefore dividing by this, so as to make the coefficient of 
dlu unity, we get 

««-it>«.2--«^i«^i»# dg ...0,— o«— J 

which must be equivalent to the first member of equation (l). 
Therefore the same expression, only with y instead of u^ 
must be equivalent to the first member of the proposed 
equation, and therefore equal to JT. Hence equating these 
equals, and reversing the operations, we find the integral of 
the proposed equation, 

y = ^ifs^if*^2 - f'^n-if—— — — -; 

(each symbol of integration affecting the whole of the ex- 
pression which follows it), which is a general formula for 
deducing the solution of a complete linear equation of the 
n^^ order, from the solution of the same equation when the 
term independent of y is supposed to become zero. 
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^y 57* We shall now give some applicatioiis of this method. 
Ex. 1. djy + 4d,y + 4y = aijf. 

If we suppose u to denote the value of y in this equation 
when a = 0, since the roots of the auxiliary equation are 
-2, -2, 

ih which expression the coefficient of djf* is 6*', therefore 
dividing by this coefficient, and changing u into y, we get 

Ex. 2. dj y - 2 wd,y + (w* + »*) y = -y. 

If we suppose e« to be the value of y in this equation 
when X^Oy we have 

u = C16"' cos wa? + C2e™* sin nx ; 

-•. dividing by ^"•cosw^, and differentiating, we find 

d.e""' sec n^ . f^ = CgW sec" nx ; 

.*. d, cos* »a? d^e""** sec no? . t« » 0, 

in which expression the coefficient of ^u is e'"' cos n<r, there- 
fore dividing by this coefficient and changing u into y, we 
get 

c^'sec wd7d,cos*na?d,c^""secnar.y = X\ 

.'. y ^^'cmnw j^ sec? nxf^e'^' cos nxXf 

or, if we suppose the constants to be added after each inte- 
gration, 

y=c"'(C' cos wa; + Csin Wi'p)+e""cos«a?jr sec*«a?jic""'co8 wa^^V. 
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Henee the solution of diy + n^y » JT is 

y s Ccosnw + (f sinn/v-f cosnw fgSec^nxJ^coBnaX, 

Ex. 3. Let Jr s i4 COS (m^r ^^ a) ^^ A cos {nw + /S), 

then JgCo^nxX 

a — r i \m sin («>af -I- a) cosnar - n cos (WiT + a) sin «jr{ 

+ — < — sin (%nw h- j3) + a?cos/3>? • 
2 12« J 

••. jCsec^w^jLcoswarufss z ; 

■' "^ i»* - »* cos nw 

• 

gf ^8in(nay-f ^) sin/3 ^ ^ 

+ ~^{ — ;-tan»^>; 

,•. tf = C cos naf -f C sin «^ + — r ; cos (mx + a) 

nr ^mr 

-f — ^ sin (»<p + S)* 
2n ^ '^^ 

Ex. 4. rfjy - 10(Cy + 62d^y - 210d«y + 26ly = e*. 
The roots of the auxiliary equation are 

2±5\/^, 3^ 3; 

.-. u ss (a + Oioi) e^ + c** (6 cos 5a? + sin 5a?) ; 

•. d" (e^^'w) = -«-•{ (246 + lOc) cos 5a? + (24c^ 106) sin 5w] ; 

.'. d, co^bwdgtf sec 5a?<? (c'^w) = ; 

.'. e^BecSwdg cos^ 5a7 d^e* sec 5a?dJ {e'^y) = «• ; 

- cos 5w -{- 5 sin 5a? 



.% cos'5a;d«( ) = «-' 



26 
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•*^ ^' 26 V cos 5a? corf»5W 

.'. e* sec 5a?d? (c'^y) = -T ; 



.-. y s — ,- 4* (a + a,J7) «*• + e** (6cos 5^ + c sin 5af). 



Ex. 5. djy + p^dlly + |),cCy + |>,d,y + p^y = -Y; 
where the coefficients are such that the auxiliary equation 

has two pairs of imaginary roots. Then the solution of 

4 

d^u + Pjdjw + Pjrfjw + p^dgU -f p4tt e 0, is 
4^u^ (a 4- a*a?) cosn^+ (6 + A^o?) sin^a?; 

.'. cos'na7a«( ) es — (1 + cos2nd7) H — sin 2nar 

H- « (6 + 6*0?) ; 

•'• d*( ) *= " ^^' *^" 2^^ + w*^ (l + cos2na?) ; 
.*. sec'na7d«( ) » - 2f»a* tannor + 2ni'; 
.'. d^sec'fiJ^d'tf ( ) = - 2»^a* sec*«a?; 

.*. dgCo»^nwdg nec^nwdg( )«0; 
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and as the coefficient of d^u is ^' cosnw^ dividing by this 
coefficient and changing u into y^ we get 

e'^'secnwdgCOffniVdgSec^nwdgCOffnwdgSecnar^'y = JT; 

.•. 6*'y = cosnir^sec*na? f^cos^nw f^sec^nofJ^cosnafe^'Xi 

or, if we suppose a constant to be added after each integration^ 
we may add to this value of e^'jf, the terms 

(a + o^^) cos nor + (6 + b^ai) sin n^. i 

Ex. 6. Suppose the linear equation to have constant 
coefficients, and let — fli* — a,, ••• — a, be the roots of its 
auxiliary equation, so that 

w* + pjm»~^ + |),m»"* +..•+!>, «= (m + aj (w + Oj)... (m + a,), 

then u «CjC^*»*+C2e"^* + Cjer«8«4....-|.c.e"""«*; 

H- c, (flg - a.) (a, - a,) eC'*'^^-)', 

&C. 8 &C. 

therefore, dividing by e^*, which is the coefficient of d]Jw, 
and replacing u by y, we get 



■MaBM.ta^AMMM<.*M..^fe.«ate«^tari.kHlMM>-^ 
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Method of Parameters. 

58. There is also another mode, of integrating linear 
equations, which deserves to be mentioned, called the method 
of Parameters, which we shall now explain. It may be stated 
thus. The complete integral of the linear equation of the 
n*** order will be of the form 

where u^j t^g, ... u^ are the n particular integrals of the 
equation when the term independent of y becomea zero, and 
fj, i>29 ••* Vii A^^ functions of ^, determined by equations of 
the form 

Let the proposed equation be 

and suppose y = v^u^ + v^u^ -f- ... t>«tt„ = 2(t?w), 

and as we have made only one assumption respecting the n 
independent quantities v^, v^ •••^n? we may make n — 1 more; 

now cf^y s 2 (vdgU) + 2 (ud^v) ; 

.-. d^y = S (vd,w), 
putting for the first of our additional assumptions 

2(wd,tj) = 0, (1); 

similarly djy = 2(t>d;w), putting 2 (d,w d,t>) = 0, (2) 
&c> s &c. 
and d;;y « 2 («cCw) + JT, putting 2 {dr/^ud^v) = X, (w). 
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Substitute these values for y, d^y^ ^mVj ••• ^y^ in the 
proposed equation, and the first member becomes 

«,{d;«,+p,cC"^Vi+...+p,«i)+t^t(cCwj|+l>i<C''«*i+— + P.«i) 

which manifestly reduces itself to X^ since each of the quan- 
tities within brackets is equal to zero ; hence if the parameters 
^i» ^s) •*-^» ^ determined subject to the above n equations 
of condition, the assumed value of y will satisfy the equation ; 
and since there are n equations in which the n quantities 
4»Vi9 dg929 ... dgV^ enter to the first degree, each of the 
parameters Oj, v^y &c. will be determined by an equation of 
the form <l«v, »/j (^r), and so n arbitrary constants will be 
introduced. 

Ex.1, d^y + n*y = secn^. 
The solution of dS y + v?y «= is y^a^ cos nw + a, sin nx. 

Let .'. y «= V| cos nw + v^ sin nzi? ; 

•'• d^y «= - t>j« sin noB + «,« cos nw^ 

making d^v^ cos nzr + dgf)^ sin nor s o. 

djy s - ©jti* cosn^ - v^n* sin nw + sec fi«, 

making - dgV^n sin no? + d^v^n cos no? = sec n«v ; 

1 sin nw 1 

.-. d,«, = ^ , or tj, B — _ loff cos nw + C,. 

* n cosnA? «' ° 



1 0? 

dmVo = - , or t?8 as — + Cft ; 

n n 

^ sf~ log cosn/r + CA cosn/p + (- + Cj] sinna?. " 



II iBii ki > II «ni tm i^^^^^«ps^r-w«i^K^Mi«qaMW«av«aBP«BaaHHn9^BB^pmPPw 
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Ex. 2. d^y + p.dT'y + Pst^r*!^ + ... + p«y « A; 

where the coefficients are constant, and the auxiliary equation, 
Af=0, has none of its roots -aj, - «j> ... — «ii» equal to. 
one another. 

Here «i = e""»', W8 = c""«', ... «. = «"•••; 
and in order to determine the parameters Vp v^, &c. assume 



-^ -/ 



^ ^/ 



then, by substitution, the n equations of condition become 

and from these equations, we have to determine each of the 
constants o^, <?s^ ... c^, by diminating the rest; for this pur- 
pose multiply them in the reverse order by the indeterminate 
coefficients ], q^y 9,, ...9^.1, and take the sum; then 

+ ... + c, (ar^ + ?jar* + -. + 9-1) -^ ^ // • 

Now assume fj, q^f ••• 9«-i so that the coefficients of 
«29 ^39 ••« ^» ni&y ^ch equal zero; then iig, 03, ... a, are 
the roots of the equation 

therefore, replacing the coefficient of c^ by the product of its 
«imple factors, 

/ 1 / 1 " * 

:? /• "r^rV^X*^'^; similarly, «,=-p7-jjji«««*Jr, &c.; 



' 



7* 



•^ 



'^'•&:-^^^/;c.; 




This method eridentlj fails when if » has equal roots, 
rince one or more of the quantities d^a^M, d^^^^M^ be. 
in that case hecomes zero. 



Ex. S. 
Here 



rfty-y^e*. 



if = f»*-l =(m + l)(iii-^ l)(in + \/^)(m-V^) = 0; 
therefore the four values of f,^X are 



3 



£ 

-^^ 



^ and the four values of 



1+V-l 



are 



1-V-l 



^a. 




,— # 



I'Z-i 



g«v:5 



-4 









• + 



«»v-i /e»(i-viT) 



-1 /e»vi-v-i; \ 



g+c,.-'- 



— + CgC* + — + C3 COS a? + ^4 sin J? 

4 4 



Sc* ^pe* 



8 



+ Cj^"' + c,^ + C3 cos ^ + O4 sin w. 



Ex. 4. «?J y — y = cos ^ ; 

the four values of f:,e^X are 



I—* 



— (co8d; + sino^) -fc,, — (<^cosd7 + sin^^) + Cjj, 



^mmmtmt 



-' — >». 
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^ (^ + Bin a?e*^^) + Cg, ^{^ + sin a?^-* ^ + C4, 

and the values of yyK^ are the same as in the preceding _£ 



example; '^^c-^ 



1 . c 1 , c 

■••*y = - (cos w + sin a?) + — c~' + - (cos a? — sin ^) — -ef 

(^e-« V=i 4. gin Off) i-=r e-*'^ 



S^/^ 4\/-l 

+ p= (j?e*^ + sin ^) + — ^ e»>^» 

8\/-l 4\/-l 

coso? ^rsincV , , 

+ — 7=:(cosa? + v -1 sin«r) 7==(cosa7-v -1 sin a?) 

4\/-l 4\/-l 

<r sin <r 
=s ya^ff '\' a^e"" + czs cos x -v a^ sm a/. 

4 

/^ 59. When the auxiliary equation contains equal roots, 
we may, by the following theorem, depress ihe proposed 
equation to one whose auxiliary equation shall contain only 
unequal roots, and then, complete the solution by variation 
of parameters. 

If <Cy + PidJ"V + P2^"^y+ ••• +Pny = -y« be a linear 
equation of the n*^ order, with constant coefficients; and 
— fltj, — ttg' •••—«« the roots of its auxiliary equation, so 
that, 

m* + |>if»*'*+P2?»"~*'+...+p„=(f» + ai)(w + ag)...(«» + o,); 

then if we multiply both sides of the proposed by c"*, - a 
being any root, and integrate, the result will be a linear 
equation of the {n — 1)*'' order with the same auxiliary equa- 
tion, except that it wants the factor m + a. 
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For let m»-> + 9,«»*"* + — + 7«-i «= be this auxiliary 
equation. 

Now d^e'^idi'^y + gi«C"*y + — + gi.-i») 

the auxiliary equation being i»""* + gjm"~* -f ... + 5^,-1 =« 0, 
deared of the factor m + a ; and if this factor should occur 
r times, then repeating this process r times, we should get 

drry + hdrr'-'y + ... + /...» = 6-~//e"jr., 

whose auxiliary equation diflPers from JIf « 0, only in wanting 
the factor (m + a)**. 

60. If we know a particular integral of a linear equation 
of any order that has no term independent of jf, we may 
reduce it to another equation of the same kind, of the order 
immediately inferior. 

Let y B ti be a particular integral of /(d,) y » 0, 

or d^y + Pxdr^'^y + P^d^s'^y + — + PnV = 0, (l) 

and assume y^uf,Zj z being a new variable, a function 
of ^; then, separating the symbols of operation from those 
of quantity, we have 

dr,y^(d, + dlYuf^z, 
d, being understood to affect u only, and d^ to affect fj,z only. 

Hence, substituting for the differential coefficients in the 
proposed, by this formula, and again separating the symbols of 
operation from those of quantity, we get successively 

(d, + di)"wJ^«+Pi(d, + di)'-*ttjL« + &c. +PnUjifZ 



m^^wmm 



11 



1 1 • Z 

+ 1 dr^xf-^ (d,) u + dl-^% . u. 



n 



But since «^ is a particular value of y, fidj) «^ «= ; hence, 
reversing the order of the terms, and observing that /^'^(d,) 
denotes the same function of d,, that dj/(i^) does of v, we 
get for the depressed equation, 

udl'^x + (ndgU '{•PiU)dl'^x 

+ <-— — cgM+(»-l)p,d,w+P2«^>d;"'*«+...+/(d,)w.«r=0. (2) 

Similarly, if we know another particular solution, u^, of 

equation (l), then d, ( — | will be a value of z in equation (2), 

and we may depress this equation to another of the same 
form of the (n - 2)* order ; and if we know r particular 
solutions of equation (l), we may in this way depress it to 
an equation of the same form of the (n — r)^ order/ 

A' 61. It y =^u be a particular integral of the equation 

djy + -Pd,yH-Qy»0, 

the other particular integral will be y = w / — —. 

For if we denote the two values of y by t^ and v, we 
have 

c^w + Pd,u + Qw = 0, c^v + PdgV + Qv = ; 
.-. ud^v — V€^u + P(udj,v - vdgU) = 0; 
.". udgV -vdgU = Ce~^'^; 



... dJ-) =lce-/-^ 

.'. V =! Cu / ^;— 
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Cliange of Independent Variable. 

/2.^ 62. Besides linear equations with constant coefficients, 
very few equations of the second and superior orders admit 
of being integrated ; and those only by particular methods. 
Sometimes the integration of an equation may be facilitated 
by changing the independent variable. 

Ex. ^.y-^a (d,y)« + 6 j? {djyf = 0. 

Make y the independent variable, and consequently for 
^#y> ^y> write 



1 


<* 


dyW 


id,a,y ' 


9 


d?a? a hx 

■— 4 -4- 


• • 


{d,wy {d,wY {d,wf 




or d'oT - adyW — hx 



= 0, 



which is integrable. 

63. In the above example there is no difficulty in fixing 
upon the new independent variable. In other cases we must 
consider x and y as functions of a third variable t^ and sub- 
stitute the values of d,y, d*y, &c. corresponding to that 
supposition. When the equation is thus generalized, we noust 
assume for x or y some known function of t, according to the 
circumstances of the case, so that there may arise for deter- 
mining the function of t which expresses y or x^ a differential 
equation simpler than the proposed one ; between the integral 
of which and the assumed function, if we eliminate ^, we 
obtain the required relation between x and y, 

Ex. 1. (1 -^)d^y- {Jod^y + n^y = 0. 

The generalized equation is 

^^ diXdN -d^yd^x d^y _ 
(dixy dix 



Itamt^.-^- .i. » A. -.- .-3ri M« H ■■ Iril ■ 
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Let w = cos ^, then diso = - sin /, dj*w == - cos t \ 
therefore, by substitution, diy + vi?y = 0, which gives 

y- A cos n^ + -B sin n# ; 



.'. y = j!< cos(wcos-^/i^) + J? sin(n cos .a?). 

Ex. 2. ^^y + — -— d,y -h y = 0. 

The generalized equation is 

diw d^y — diy d^x A diy B 

■■■ ■ -|- ■ -^ y as 0* 

(d^a?)^ a + bw dios (a + 6^)^ 

Let d/^ = a + 6^, or 5^ = log {a + fca?), or ^^* = a + 6^ ; 

then d^x t=s bdtw; 
therefore, by substitution, 

d^y+{A-b)dty^By^O', 
the solution of which is 

y B Ci€^* + CgC"'* = c, (a + 6«r)^ + Cg (a + 6^?)^, 
m and m being the roots, supposed possible, of 

m* + (-4 - 6) m + S = 0. 



If the roots be imjpossible, and of the form m :^nv -1, 
the solution is 



m 



y = j3e*"* cos (»^ + a) = /3 (« + 6a?)* cos {log (a + 6<2?)'5 + a} . 
Ex. S. (a + fea?ydJy+il(a + 6^)^d5y4.jB(a + 6fl?)d,y+Cy 
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As in the preceding Example, 
let a + bw=^, .-. diX = ^, 

d;y-e-'"(«i«y-6d.y), 
«ey - e-^Ci^y- 86d«y+ iVd,y) ; 
.'. by substitution we get 

d?y + (1 - 36) c^y + (2^ - i<6 + JB) d,y + Cy =/ (— ^), 
a linear equation with constant coefficients. 

The same substitution of course succeeds for the equation, 

which, by putting a ^-bx = z^ may first of all be reduced 
to the form, 

«"d;y + a^s^-^a^-^y + &c. + a^.issd^y + a,y «/i (z) ; 

and then to a linear equation with constant coefficients, by 
the formida 

«"djy = d^y-p^d^^^y + &c. ±p,_iy, 
where «=^^^ and pj, p^t &c- ^^ such that the roots of 

k" - Pi*"-* + PgA;""^ - &c. •.. ± p,.i* = 0, 
are 0, 1, 2, 3 ... (w - l). 

Ex. 4. a?^cgy + wd^y - y = -- — -, 
' ar — 1 

. The substituticHi a ^^ gives 

<^y-y=« + p5^; 
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which agrees with the above, upon writing a for ^. 
Ex. 5. OT^dJy 4- adgy — y = aofy 



y = — + C^cB + 



a? " w* — 1 



Solutions expressed by Definite Integrals. 

/' 64. Sometimes the complete integral of a differential 
equation may be expressed by a definite integral, as in the 
following instances. 

Ex. 1. dj'^y B o^y + m. 

or* - — 

Let t^i >= Oi / f^^^^.e *, a^ being a constant quantity ; 

then cC"*«i ~ ®* / ^" «"'**« » = oj (1 + ^©i), 

change a^ into a,) as, &c., and let v,, f^s, &c. be the cor- 
responding values of Vi, 

.-. d;-'«j,— oJO + 4?©,), 



Now suppose aj, a,) •••a, to be the n*^ roots of unity; 
then multiplying each of these equations by an arbitrary 
constant, and taking the sum, 

*"* iP\^\ + CgVj + ... + c,©,) = c, + c, + ... 

let C| + Cjj + ... -^rc^-m^ 
11 
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then y =» OiV^ + CjU, + ... + c.», is the complete integral, con- 
taining only « - 1 arbitrary constants by reason of the equa^ 
tion of condition ; or it may be written 

/*a> 

a being a primitive root of ft" - 1 » 0. 

Let 9i» s 0, n « 3, then cJ^y « /ry, the solution of which is 

cos- 



/•oo f ^^^* 






Let m » 0, n B 2, then df,y « ^y, and 
ysC / e ^(g"-Hg )« y — / e*cos/i?i!f, putting ^=»v -1, 

= C\/27rc2, (Integ. Cal. Art. 112). 



which agrees with the result obtained by direct integration. 

If the more general form were proposed, 

cCy = (a + fiai) y, 

assume a + fiat ss fi» t; .*. d^jf s c^y d«^ s fi^d^y, 

•-1 11-1 

and d;"'y = j3 • dj-'y ^ fi * ty; 

and to this form may the still more general case 

di'^y « awy + 6^ + cy + c 

be reduced. (The above solution is taken from Crelle^s 
Journal, Vol. x.) 
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Ex.2. diy+fn^-^y^O*. 



Let y SB — where i» (m + l) « «> and suppose a to 
be such that m is a whole number; then 

w (d*w + n^u) = ^fndgU. (l). 

Differentiate this 2, 4, &c. times successively^ and add 
each result to the preceding multiplied by w', in doing which 
the symbols of operation and quantity may be conveniently 
separated, then 

and wff (dj + n*) w « 2mn^d«t^ ; 

.-. ^ (d|J + w«)» w « 2md, (d^ + w') w - 2d, (d^ + n*) w 

= 2 (m - 1) d, (d]J + n*) « ; 

similarly, a? (d|J -»- n*)* w = 2 (f» - 2) d, (d^ + n^Y u, 

and repeating the process m times, the second member becomes 
zero, and 

(dl + nT^'u = 0. 

Hence we have, for the determination of Uj a linear 
equation with constant coefficients, whose auxiliary equatiofi^ is 

.•. f« «B (a + Oi^ + «2^ + ••• + o«^*) cosn^ 

+ (6 + biW + b^a/^ + ... h^iif^) sinner, 



* This equation, when a = 6, presents itself in investigations relative to the 
figare of the Earth ; the solution in the forms here noticed was given by Mr Gaskin 
(Cambridge Problems, 1839). 
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this gives the form of the solution; the arbitrary constants 
must be reduced to the proper number, in any case, by 
actual substitution ; thus let m » 2, so that a^^S; then 

u B (a + a,/r + Os^ cosnw + (ft + biW + b^a^ mnaf, 

and substituting this value in (1) in order to reduce the con- 
stants, we find 

ary = u^ a [cosnat (1 1 + nx annw^ 



+ 6 {sinn^p ( 1 1 -nay cos na?|, 



for the complete integral of 

6 



<y+ (»*-^) » = ^- 



But when in the equation m (m + 1) = a, m is a frac- 
tion, 

make y = oT^^Uy then w (d/u + n*w) + 2 (m + l) d,tt = ; 
the complete integral of which is 

a and /3 being the arbitrary constants ; for this gives 






I3w 



.S /•» 



= -7^ r~ f (<» - n*)*+* COS («< + o) ; 



dlw-li' f (f - n*)" {- i*)cosimt + a) 



^^~ f - (<» - n«)- (<* - n* + n») cos («< + a) ; 



^b»^p r ■»>< 
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-• /•• 



2(m+l) 



Similarly, when m is an integer, it will appear by sub- 
stitution that 

f cos (^ y/t + a) \ 

u - /3C.. { ^ J 

is the complete integral of w(dlu -^n^u) ^Zmd^Uj the re- 
duced equation when «ii is an integer. 



Simultaneous Equations. 

65. In these equations, which sometimes occur in the 
higher parts of Dynamics, instead of one equation between 
Wf y, and the differential coefficients of y with respect to /r, 
being given to determine the relation between ^ and y; we 
have two equations containing ^, y, t (of which w and y are 
considered as functions) and the differential coefficients of a 
and y relative to t, to find that relation. 

66. To integrate the simultaneous equations of the first 
order, 

aw-^ by + dtw = T, a us + b'y + dty = T', 

T and T^ denoting functions of t Multiplying the latter 
by an indeterminate quantity m, and adding it to the former, 
we get 

oa? + fey + m (o'a? + b'y) + dj (a? + my) « T + mT', 

or d( (w 4- my) -i- (a + ma ) (or + j y) = T + mT". 

a + fna 

Let -^ = m, which will give two values of m, mj 

a + fna 

and mg ; then the equation, under this condition, becomes a 

linear equation of the first order; and we obtain by inte- 

gration 

(a? + my) «(•+««'>'« jj(r+ mT") «<*•+•"«'>•; 



■«9» 



.,j^ J I II— ^^>y^w>*nw^Bi^F— ■^^^^— *"* ■*" 
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and by substituting successively the two values of m, we 
obtain two primitive equations, which will furnish values of 
Of and y in terms of t^ and the relation between ^ and y^ 
if ^ be eliminated. If the two values of m are equal, we 
shall obtain only one equation between ^, y, and /; but if 
this can be solved with respect to w or y, and we substitute 
the value so found in one of the given equations, we shall 
obtain a second relation either between w and tf or between 
y and t; and then t may be eliminated as before. 

Ex.1. 5w — 2y + diW^ ^9 Gy-of-i-dty^^; 

.-. (5 - m) ^ + (- 2 + 6m) y + di (a? + my) = e* + wc^, 

-^2 + 6m 

let s m, or m* + m - 2 = 0, or m « 1, or - 2; 

5 — m 

.-. d| (a? + y) + 4 (a? + y) « e* + 6**; 

^56 

^ 2 
similarly 4?-2y« 6** + C\e"'', 

which determine ^ and y in terms of U 

Ex. 2. dj^ + 5d? + y = e*, dty + 3y -- w^ e^; 

.-. dj(a? + my) + (5 -m)a7+ (l + 3f»)y =6* + mc^; 
1 + 3m 



5 — m 



= m, or 1 - 2m + m* = (l - m)* = 0. 



Hence the values of m are each » 1, and integrating, 
we find 

a? + y a= Cif?"** + -e* + -e**. 
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By means of this, eliminate y from the first equatipO) and 
we get 

5 o 



... .V = ± c* - 4«** - Cite"* + C,c-^ 
25 36 ^ 



and y = — c« + -^e»'+Ci(l + Oe"*'-C,c-*'. 
^ 25 36 IV 7 s 



The more general form 

may evidently be reduced to the above by successively elimi- 
nating dtfff dtW. 

&J. To integrate the simultaneous equation of the second 
order, 

o/r + 6y + c + dfo? = 0, a'm + b'y -»- c' + d^y =• 0, 

Midtiplying the latter by an indeterminate quantity m, 

and adding it to the former, we get 

« 

(a + ma) a? + (6 + mb') y -^ c + me' + c^ (a? + my) « 0, 
or dj (47 + my + c^) + (a + ma) {at + my + Ci) = (l), 

if ,^m (2), > = Ci; 

a + ^a a-^ma 

therefore, integrating equation (1), and substituting succes- 
sively the two values of m given by equation (2), we obtain 
the two required primitives; or if the values of m be equal, 
we must proceed as in the former case. 



Ex. 1. dfw - (S4f + 4y - 5) « 0, dfy + (« + y + 5) « 0. 

Let Bin, or m*- 4» + 4 « (m-2)* = 0; 

m — 3 

.-. d^ (/r + 2y) - (jr + 2y - IS) = 0; 

4f + 2y — IS = c«* + c'c"', 

and eliiiiinatiiig w from the latter of the given equations, we 
find 

.-. y « 18 — (^ - i) c* + - /e"' + oc* + a'c"*, 

2 2 

and d7 = - 2S + c (^ + ^) 6* - c (^ - 1) «"* - 2ac* - Stae^K 

Ex.2, dlaf + adtOf ^ Oy d^ y -^ adtff + b ^ 0. 

68. If we have three variables a, y, itr, which are func- 
tions of tf and if 

diOf -^ aw -^by + c «^ Tj 
c^y + a,ar + biy + Ci» = Ti , 
di» + a,47 -I- 6,y + c,« « r,, 

then if we multiply the second and third by indeterminate 
constants m, m\ and add them to the first, and assume 

6 + b^m 4- b^fn » m (a + aim + Osm'), 

c + Citn + c^m! ^m {a-^ a^m + a^m^, 
we have 

dt {a + my + •»'«) + (o + mo, + m'a^ {x + my + m'x) 

= T + mr, + m'r,. 



— - — .« ■■Ly,^^—'" I «J I II Vm I I II i^B,^ ^-^ . * 



89 

To determine m, m, we obtain cubic equations; hence 
if wii, wig, W3, wi/, m^\ wig' be their roots, we have, by 
solving the linear equation, 

J? + may + m^'x = Fg (^), 
from which equations a?, y, % may be found in terms of t. 

Approximate Solutions of Differential Equations. 

69. When all the known methods of integrating a 
proposed differential equation fail, we must endeavour to 
resolve it approximately, that is, to obtain from it the value 
of y in terms of ^, in the form of a series. The first mode 
which presents itself of effecting this, is to assume for y a 
series arranged according to powers of a?, with both its co- 
efBcients and exponents undetermined; for in most cases it 
happens that the exponents do not follow the progression of 
the natural numbers, and that particular artifices are requisite 
for discovering their law. When the form of this series is 
known, we determine its coefficients by substituting it and its 
differential coefficients, for y, d^y^ &c. in the proposed equa- 
tion. The following application of the method to Riccati^s 
Equation will give an idea of the mode of obtaining both the 
exponents and the coefficients. 

Ex. djgy + bi^ - aa?" = 0, 

d X 
let fcy = -^; .'. dl% - ab% of" szOj or putting c for -ab^ 

dlx + cxaf = 0. 
Assume z - a}'^{A + Bafi + Ca!^^ + &c.) ; 

+ (a+/3)(a+)3-l)i?a?«+^-^+(o+2i3)(a+2)8-l)C^«+*^-*+&c.] 

+ cAx^-^"" + cSaf«+^+'^ + &c. 

= 0. 
12 
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Hence j3 must equal n + 2 ; and then to determine 
a. A, By &c. we have 

a (a - 1) ^ = 0, 

(a + n + 2) (a + n + 1) 5 -f c^ = 0, 
(a + 2n + 4) (a + 2n + 3) C + c5 = 0, 
(a + 3w -f 6) (a + 3n + 5) 2> + cC = 0, 



Hence a may be either zero or unity, and A remains 
undetermined; calling therefore A and A' the two values of 
A corresponding to a = 0, a = l, we get 

X sa A \i — + — ^ — JiC J 

^ (n + l)(n + 2) (w+l)(?i+2)(27i+3)(2n+4) '^ 

"^ ""(w+2)(w+3) (n + 2)(n+3)(2w+4)(2w+5)" ^*^' 

and substituting this value of % in the expression y = -^, 

we shall obtain the value of y, involving only one arbitrary 

A 
constant — . 
A 

As the terms of the above series have divisors of the 
forms (w + 2) i sf Ij where i is an integer ; if n be such that 
(n + 2) i sf i = 0, one or other of the series will be illusory, 
and we shall only obtain a particular value of %; and if 
n + 2 as 0, both' series become infinite, but in that case the 
equation may be exactly integrated by Art. l6. 

70. In the preceding instance we arrive immediately at 
a complete result; but it often happens that the solution we 
obtain by the method of indeterminate coefiicients involves 
no arbitrary constant. To supply this defect, we must in- 
troduce, instead of the arbitrary constant, a value of y cor- 



ft 
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responding to a given value of /v; that is, supposing these 
to be b and a, we must substitute in the given equation, 

y zsb + Uy w = a + t; 

then determine u in a, series all whose terms vanish when 
^ = ; and replace u and t by their values y — 6 and «a? — a ; 
in this way it is evident that the arbitrary constant will be 
involved implicitly; for, from the complete integral 

/(^, y, C) = 0, 

C may be expressed in terms of a and 6. 

Ex. d„y + y^g(xr, 

this becomes d^^u + 6 + w = ^ (a + ^)"* ; 
assume u = t^ (A + Bt^ + Cty + &c.) ; 

+ b + ^^ + Bt^^^ +&c. 

m(m-l) 

-ga"" - mgoT-^t ^^ ^ga"-^^-&c. 

1 • z 

.-. a = 1, /3 « 1, 7 = 2, &c. 

^ =^a"' -6, ^B^gmaJ^^^ - go^ + 6, 

6C = ^m (w - 1) a"-* - gmoT-^ + ^a" - fc, &c. 

71. The approximate solution of a differential equation 
may sometimes be obtained in the form of a continued fraction 
by assuming 

^'^TT 1 + 1 + &c. ' 

First, suppose. <i? to be very small, and for y substitute 
Aof^ in the given equation; then, retaining only the terms of 
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lowest dimensions in ^, A and a become known by equating 

coefficients and exponents. Next, write «« in the 

1 + » 

proposed equation, and in the result put % » Bafi^ and de- 
termine B^ jS, as before, by supposing ^ to be very small ; 

then in the transformed equation in x^ put z^ ; and 

SO on for the rest. 



Ex. my + (1 + a?) d^y •= 0. 

Let y^Jof^; .-. (f» + a)J^ + -i<a^"^ « 0, 
or Aacs^"^ =0; .-. a = 0, 
and A remains undetermined. Next, put 

y « , and » = Bafi^ 

and we get successively 

m (1 + «) ^ (I -¥ w) dgZ ; 

m + B^^(wi-)8) = i35a?/3-», or m^Bfiafi''; 

.-. /3=1, Bsitn; 
similarly, putting 

iif « , and t = Cay'y, 

we determine C and 7; and thus we obtain 
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Since the proposed equation when integrated gives 

y = J (1 + ^)-"', 
the above continued fraction is the development of A (l 4-/p)"". 

72. We may approximate to the integral of a diiflfer- 
ential equation by successive substitutions, in a manner 
similar to that invented by Newton for the solution of 
algebraical equations, as in the following instance. 

Ex. d^y + n^y + ay* + a « 0, where a is a very small 
quantity, 

we may assume 

y ^u •\- aUi + a^Ui + &c., 
which gives 

+ a^ (d^t^ + ^*% + SwWjj) + &c. = 0. 

Hence, equating the coefficient of each power of a to 

zero, we get 

» 

c^^u + n^u + a = 0, 
1*2+71^^2 + 2^^1 = 0, &c. (3.) 



a 

The first give u ^ r + c cos nof + c* sin nw ; 

n 

and this value substituted in the second reduces it to the 
form 

dlui + n^Ui = JTi, 

the integral of which by Art. (58.) is 
Ul'stcosnw(c^ fgX^ sin no) + sin nw (c/ + —fgXi cosnaf). 
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Similarly, these values of u and t^, substituted in (3) 
reduce it to the form 

which may be in like manner integrated; and in this way 
the coefficients of the powers of a may be deduced one from 
the other by a uniform process. 

73. We have seen (Art. 57.) that the solution of 
^iy + ^*y - ^ COS ^^ + a) + B cos (»a? + /3), 
is y = Ci cos nx + C2 sin no} 

+ — r cos (ma? + a) + — a? sin (nx + /3). 

Hence, if from the proposed equation we had to deter- 
mine y approximately, we could not neglect the term 

A cos {mw + a) 

even when A is exceedingly small, provided m and n are 
nearly equal to one another ; because in the value of y this 
term is divided by n^ — m^ which is very small. With re- 
spect to the last term in the value of y, we remark, that 
it is not periodical, but may increase indefinitely, as <r in- 
creases. The equations which present themselves for solution 
in physical Astronomy are usually of the above form; and 
upon the peculiarities just noticed depend some of the most 
interesting results in that subject. 



SECTION V. 

ON DIFFERENTIAL EQUATIONS INVOLVING TWO OE MORE 

INDEPENDENT VARIABLES. 



7^* I^ ^Iiis part of the subject, it is convenient to 
employ a system of notation which we have not hitherto found 
it necessary to notice, and which may be explained as follows. 

If t^ be a function of any number of independent variables 
^9 y^ %9 ty &c., and it u '\-^u be the value of u when these 
variables simultaneously become 

OD + ^x^ y H- 5y, % + 5^, &c., 

then ^u = dgU.^w + dyU» ly + d^u.Sz + ... 

+ terms of two dimensions in Sof^ Sy, &c. 

If we now agree to 'neglect all terms of higher dimensions 
than one in 5a?, Sy, &c , and denote the value of Su cor- 
responding to that supposition by dw, we have 

du c= dgU.Sof -i- dyU.Sy -{■ dgU.Sz + ... ; 

according to which definition, it appears that du denotes that 
part of the increment of m, as given by Taylor's Theorem, 
which involves only the first powers of the arbitrary incre- 
ments of the variables on which it depends ; hence, in con- 
formity with this definition, 5a?, 5y, Sxy &c., must be repre- 
sented by dwy dyy d%^ &c, ; for if f(jxi) = «i?, then 

f{w + hoi) = .r + Ina ; 
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consequently, that part of the increment of f{w) which in- 
volves the first power of Sat is, in this case, the whole of it, 
Sofy therefore, dof « S'.v ; and so on for the others ; 

.'. du = dgU . dw + dj^u . dy + d^u . d% + &c., 

where dw, dy, dz. Sec, are the arbitrary increments of the 
independent variables ^, ^, z^ &c., and are entirely inde- 
pendent of those variables; and du is that part of the 
corresponding increment of the dependent variable u, which 
involves only simple dimensions of dw^ dy, dxy Sec; and 
which approaches nearer and nearer to the value of the whole 
increment of u, the smaller do?, dy^ dz^ &c are taken. 

The quantities d^u.dwy dyU.dy^ Sec, are called the 
partial differentials of u with respect to <r, y, &c respectively ; 
and du h called the total or complete differential of u^ or the 
differential of w, merely. 

75. According to the above definition, the differential 
of duy or the second differential of u, will be, supposing it 
to involve only two independent variables, 

d?u = dg (du) dx + dy {du) dy 

m 

= dg{dgU . dx + dyU . dy) dw + d« {d^u . dx + dyU . dy) dy 

= d^w . dd^ + 2dgdyU . dwdy + (^^u . dyS 

because in this process dw and dy are entirely independent 
of X and y ; being, in fact, the arbitrary increments of those 
variables, which might have been denoted by h and A?, did 
not a due regard to the precision and symmetry of the 
notation require it otherwise. And in general, if t^ be a 
function of any number of independent variables ^, y, z^ t. 
See, the n* differential of u will be given by the formula, 
(where the symbols of operation are separated from those of 
quantity,) 

d"w = (da? . dp + dy , dy -^^ dz . d^ -{■ &c)"w. 



-v-^-< 
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76. Hence, we see that the differential of a function 
has no reference to one variable rather than to another; 
but in its formation, all the variables of which it is a func- 
tion are supposed to undergo simultaneous unconnected 
alterations; and the value of the differential depends upon 
the values of all those alterations. Whereas in forming a 
differential coefficient, one particular independent variable 
only is changed, the rest remaining unaffected ; and a quan- 
tity is produced, wholly independent of the value of the 
alteration which that variable may have received. 

77' It IS evident that the rules for deriving differential 
coefficients, suffice for finding the differentials of functions. 

Having, for instance, formed all the partial differential 
coefficients of the first order, if we multiply each by the 
arbitrary increment (or differential) of the corresponding in- 
dependent variable, and take the sum, we shall obtain the 
first differential of the function ; and similarly, the differ- 
entials of the second, and higher orders, may be formed by 
means of the formula at the end of Article 75. 

78. Total integration of a proposed differential, is the 
finding a function whose differential is the quantity pro- 
posed. This operation is denoted by the symbol / without 
any affixed variable; this process, like the former of differ- 
tiation, having reference to all, and not to one particular 
variable. 

Ex. If w = ci?2 + y% 

dw = 2 (jxdx + ydy) is the differential of u ; 

u 
also f(wdw^ydy)^—9 since the latter quantity differentiated 

produces the quantity under the symbol of total integration. 

When a quantity is presented for total integration, it 
must be reducible to the form f(v)dv; and then its integral 
= f^fip)' After this observation, the rules of ordinary in- 
13 
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tegration will enable us to perforin the process of total 
integration. 

Ex. To integrate — — j \ 



r ^- 

. y {wdy - ydw) j ^ 

1 ^v^T? -/ ^ 



(!) 



wdy — ydx 



r t? , - y ' 

= / J dv^ putting - « u ; 

= - Vd^ + y" + C. 



79. Differential Equations involving more than two vari- 
ables admit of division into two classes, Total and Partial. 

A total differential equation is one which expresses the 
differential of the dependent variable in terms of the other 
variables and their differentials, and sometimes also of the 
dependent variable itself; it is consequently equivalent to a 
system of equations in which each differential coefficient of 
the dependent variable is given explicitly. Thus, u being 
a function of the independent variables cOy y, z^ the total 
differential equation 

du = Pdw + Qdy + Rdz 
amounts to the same as the system of equations, 

dgU = Py dyU=zQ^ dgU^R. 
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Also, X being a function of the independent variables 
<r and y^ the total differential equation 

Pd/v + Qdy + Rd« = 

amounts to the same as the two equations, 

P + Rd^z = 0, Q + RdyZ = 0. 

A partial differential equation, on the contrary, is only 
a relation between all or certain of the partial differential 
coefficients of the dependent variable, and the variables; as 
in the instances, 

(w - a) dgX + (y - 6) d^x = » - c, 

o^^x + 2<ry d^d^z + y^dlx = 0, 

X in both being a function of the independent variables «r 
and y. 

Total Differential Equations. 

80. To integrate du » Pda^ + Qdy, u being a function 
of the independent variables w and y, and P and Q being 
functions of of and y. 

Since Pd^ + Qc^S^ must be identical with d,u dw + dyU dy, 
we have 

dgU = P, dyW = Q, with the condition dyP^ c?*Q; 

.-. u^f^P'¥ f(y), f{y) involving y only ; 

.-. Q^dyU^dy(f,P) + dyf(y); 

consequently u is known. 

Ex. du = 7 r-„ dof - 7 dy, 

dgU^-z 72' dyU=-- ; — — , 
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" - - jf, */<»' 



^ 1 y . ^i' X 



(a? + y)« J? + y (cT + y)' 



u = C- 



of + y 



81. To integrate dt^s= Pdw -^ Qdy -^ Rd%^ u being a 
function of the independent variables w^ y^ x. 

Since Pdof + Qdy -\- Rdx is identical with 

d,«* . dw + dy w . dy + dgU . d«f, 

we have d,u = P, dy w = Q, d,w = J? ; 

together with the equations of condition 

dyP^d^Q, d.Q^d^R, d^R^d^P, 

which are found by supposing each of the quantities «^ Wy y 
to be constant in succession (as we are evidently at liberty to 
do, since those variables are independent of one another), and 
then taking the corresponding equation of condition for du 
being the exact differential of a function of two independent 
variables. Hence 

u = fg P -k- Wy w being a function of y and « ; 
also Q = dy(f,P) + dyWj 

R^d,(f,P)+d,Wy 

which two equations giving the values of the partial differ- 
ential coefficients of «i;, its value may be found by the preceding 
Article; and so the value of u completely determined. 
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Ex. 1. du = yxdw + wzdy + jcydx ; 

rf*«^ ■= yXf dytt = af»^ dgU «= ^y, 

wz = wx + dyW^ or dyU) = 0, 
^By^xy-k-dgW^ or d,w = 0, 

.'. to = C; 

.-. w = wyx + C 



y cp wy 

Ex. 2. dw = — ^ — da* + dy + - — ^-— dx ; 

a - flf a — « (a - ijf)'* 

^ y ^ ^ ^ ^y 

dgU = — - — , dyU = , o,w = 



a- X ^ a - « (a - «?)'* 



w « -^ + C. 



a — X 



Ex. 3. du =^ a{offdw -^ ydy) + bdx; 

/ 82. We next come to the consideration of the equation 

Pdof + Qdy + fidijf -= 0, 

one of the three variables w^ y, x^ being a function of the 
other two, which are independent. 

Since the proposed equation may arise from combining the 
results of differentiating two separate equations, we have first 
to examine whether it can be satisfied by a single primitive 
relation between /r, y and x. If of^ y^ z he co-ordinates of a 
point, the cases will be distinguished according as the proposed 
equation is the differential equation to a series of surfaces, or 
a series of curves. 
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For example, the equation 

(^ - c) xdw + (» - c) ydy « {a^ (^ - o) + y (y - fc)} dx^ 
arises from combining the results of differentiating the equations 

for these equations give respectively 

ydgX - xdy% = 0, (^ - a) d,« + (y - 6) c^y« = « - c ; 

from which if d^^x and d^z be determined, and substituted in 

d% = dgZ . dx + cf^jsr . dy, 

we get the proposed equation ; which, consequently, cannot 
in general be satisfied by a single relation between x^ y, %. 
It is the analytical expression of the conditions of the problem, 
to find a surface belonging at the same time to conical surfaces, 
and surfaces of revolution about the axis of z. 

83. To find the equation of condition for 

Pdx + Qdy + Rdz = 0, 

admitting a solution of the form /(a?, y, z) = C 

If the proposed equation can be satisfied by a single 
relation between a?, y, and ar, 

P Q 

dz= - — dx-^dy. 



.-. d,«r=- — , dyZ—--^^ 



is the differential of a function of two independent variables ; 

with the condition ^y (■^) = ^* ("^) ' 

or, since P, Q, R may contain z, which is a function of x 
and y, 

rf(.) (f ) + d(0 (I) d.^ = ^<^) (!) + d,,, (I) d^z ; 
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therefore, substituting for d^x, dyZ, their values, 

««„, Q - ««,., (5) - «",., (I) - ^•'c (I) ^ 

and performing the differentiations, and reducing, we get 

the equation of condition for 

Pdx + Qdy + Rd% = ; 

admitting a solution of the form /(^, y, %) = C 

84. If the above equation of condition is not satisfied, 
then the equation 

Pdx + Qdy + Rd% = 

cannot, by being multiplied by any factor, become susceptible 
of a solution of the form 

f{w, y, ar, C) = 0. 

For suppose F to be a factor which renders 

Pdx + Qdy + Rdx 

the immediate differential of some function u of a?, y, isr, con- 
sidered as independent; 

or rd(,,P + Pd^,,V = Fd(„Q + Qd^,/ 

Kd(.jQ + Qd,,/ = Fd^^^JZ + Rd,^/ > (1). 

rd(„iZ + JZd(„F=Fd,,,P+Pd,,/ 

Hence, multiplying the first of these equations by ij, the 
second by P, and the third by Q, and taking their sum, the 
factor V disappears, and we find 

P (d^^^R - d^,Q) + Qid^,,P- d,,^R) + R (d^^^Q - d(,jP) = 0, (2) 

the same equation of condition as in the^ preceding Article. 



104 

85. If this equation be satisfied, which will be the case 
only when the proposed admits a primitive of the form 

equations (1) afford a means of determining F. Then 

du = VPdof + VQdy + VRdXy 

whence u can be found by the method of Art. 81, and 
u -^ C =iO is the required relation between j?, y, and z. 

Or, without determining F, we may integrate considering 
one of the variables constant, and add an arbitrary function of 
that variable ; then differentiate the result with respect to that 
variable and compare it with the proposed equation, and so 
the correction will become known. 

Obs. In the majority of cases which present themselves, 
the factor F is capable of being determined by inspection. 

Ex. {ay - b%) dw + {cz — aixi) dy + {hx — cy) dz = 0. 

Divide by {ay - h%) {bx - cy), and the resulting equation 

dx (cz — ax) dy dz 

+ . . X ./ . + 7~ = 0, 



bx — cy {ay - bz) {bx - cy) ay — bz 

satisfies the conditions, considering or, y, z^ as independent, 

dyP^d.Q, d.Q^dyR, d^R = d,P; 

and therefore the first member may be regarded as the dif- 
ferential of some function, «, of *!?, y, z considered as in- 
dependent ; 

.*. a„u = 



bx — cy 
1 



u 



= r log {bx - cy) + w ; 
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ay — bz 

ex ^ ax c 1 

(ay — hz) {ba — cy) b* bx — cy ^ ' 

or dyW^--. --; 

" b ay ^ bz 

1 ady — fcd« 

..*. dw — - - . ; — ; 

6 ay — bz 

.-. t/7 s ^ - log (ay - 6«) + C ; 

1 - /&a? — cy\ ^ bx -- cy ^ 

.\ w = - log ( -^) + C = 0, or -^ « C. 

b \ay — 6W ay — bz 

Obs. When in an equation of this sort, the differentials 
enter above the first degree, it is not integrable unless it can be 
resolved into rational factors of the form Pdx + Qdy + Rdz; 
for whatever be the integral, it must upon differentiation pro* 
duce a result of that form. 

86. If the equation 

Pdx + Qdy + Rdz « 

be susceptible of a primitive of the form f(x, y, z^ C) = 0, 
and be homogeneous and of n dimensions with respect to 
Xy yj z; then, putting x » vz, y = wz^ and dividing by ss*, 
it becomes 

S (vdz -i- zdv) +T (todz + zdw)+ Udz^Oy 
iS, T, U being functions of « and w ; 

dz Sdv 4- Tdw 



or — 



z Sv + Tw'\'U' 



hence the second member is an exact differential since the 
first is SO) , and it may be generally integrated by inspection, 
or by the method of Art. 80. 
14 



• • 
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Ex. 1. df -^ «) da -^ (a -^ x) dy + (x -^ y) dz » Oy 

put a s vx, y s w«^ 

dz (tr + 1) d© + (t> + 1) dw 

z V (fo + 1) + tr (v + 1) + V + tr * 

div . d (o + ti' + ««?) 
or — + i— ^ ^ = 0; 

Z * 1? + «7 + ©W 

Ex. 2. {^ -^yz-^s/^dw -^^ {a^ -^^ az -^ sif)dy 

+ (^ + a?y + y*) dj?r = 0. 

Make a » vjv, y « wzy 

dz {u^ + te? + 1) dt> + (r* 4- « + 1) d«^ 

' II i^ t ...■■■ ■ . ■ 111 ■■ 8B OL 

Z V («?• + IT + 1) + «I7 (t7* + V + 1) + 1?* + t^tr + tr* 



but of the latter fraction, the denominator 

aa i?tr (m^ + 1? + 1) + V (w + 1 + u) + w (t? + 1 + w) 
— (utr + «? + fi?) (w + fi; + 1), 

and the numerator 

= (1 + tj + w) d (tJ + w + VM?) - (» + w + ©») d (r + ») ; 

d«r d (t' + ur + t>w) d (1 + ti + «;) 

°l l+v + ti; j "^ 
.*. i?a7+»y + J7y « C(« + ^ + y). 
Ex. S. (oT* - y* + «*) d^ - «*dy 

+ (y - «^) («" + ^y + 0^) — -0, 
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putting w^VKj y^wx, this is reduced to 

dw tB (1 + «• - fc^) dv ; 






= Jje-*' + C, (Ex. S. Art. 18.) 



or « L6 «» + CXf 

z being constant under the sign of integration. 



Total DiffdrentiAl Equations that do not admit of 8 Single Primitive. 

87* We have seen that the equation 

Pdw + Qdy + Rdx » 0, 

when the equation of condition (Art. 83) ii$ not satisfied, does 
not admit of being derived from a single primitive equation 
involving two independent variables. The integral in this 
<sase will be exhibited by a system of two equations; dttd 
the proposed equation cannot be regarded as the diflTerential 
equation to a surface^ but to a system of curves in space, all 
endowed with some common property. 

Ex. 1. dx a aydw + bdy. 

Since the equation of condition ih this case is not satisfied, 
w and y cannot be independent, and we may assume y-f{^\ 

.% dx ^ af (ai) dw + bf {ai) dw ; 

.-. x^af^f(w) + bf{w)j 

which, with y^f(/xi)y constitutes the integral of the proposed 
equation. ^ 

In general, if F be a factor which makes Pdw + Qdy an 
exact differential, considering x as constant, and we^d - 

f{VPdx + VQdy) « w + («) ; 
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it is evident that 

w + (f}(x)':^0, together with d^w + 0' (z) - VR « o, 

satisfy the proposed equation, where <p (x) denotes any func- 
tion of z 

Ex. 2. zdw + ofdy + ydz ^ 0, 

y -^ z log ^ + (;y) s 0, log zp + d, ^ («)=»- . 

Ex. S. {^ (a? - a) + y (y — 6)} dar = (iif — c) {wda + ydy), 
a;* + y* + 20(i8) »= 0, or (a? - a) + y (y - 6) + (« - c) 0'(iif) « 0. 

P^ial Differential Equations. 

88. In partial differential equations of two independent 
variables, the differential coefficients of the first order d^z^ 
djfZf of the dependent variable z^ are usually denoted by the 
symbols p and q; and dlz, d^d^z, d^Zy the differential co- 
efficients of the second order, by r, «, t^ respectively. A 
partial differential equation is said to be of the n^^ order, when 
it involves one or more of the partial differential coefficients of 
the dependent variable of the n*^ order ; but none of a superior 
order. To be the general equation of the ri^ order, it ought 
to contain the independent variables, and the dependent va- 
riable together with all its partial differential coefficients from 
the first order to the n*^ order inclusive. To integrate a 
partia] differential equation, is to find for the dependent 
variable, an expression between the differential coefficients of 
which, that relation exists which is indicated by the proposed 
equation ; and under the most general form possible. 

Equations of the First Order. 

89. The complete integral of /(a?, y, «r, jp, q) = 0, the 
general equation of the first order, will involve one arbitrary 
or general function. 
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For let u^F\w, y^ %y 0(^)} ^^9 be an equation by 
•virtue of which jir is a function of the independent variaUes 
Of and ffy V being a known function of w^ y, and z. Then 
p and q are given by the equations d,u «s 0, dpU « 0, each 
of which will involve d«0(v) or (p^v), and may involve 
<j> (v) ; consequently, between the three equations f« « 0, 
dgU s 0, d^u » 0, it will be possible to eliminate (f> (v) and 
0'(t)), and there will result a relation /(^, y? ar, p, q) = 0, 
wholly independent of the form of the function (v) ; and it 
is evident that in general more than one function cannot thus 
be eliminated. Conversely, an equation of the form 

f(^j Vj ssy p^ q)^0 

being proposed, its complete integral to have all the gene- 
rality possible, must be of the form F \wy y, «y 0(v)} »0, 
where the form of 0(t)) is perfectly arbitrary. 

For example, let 

w «= « + fnw + «y - (©) a 0, 
where t> = (^ - a)* + (y — 6)* + (» — c)* ; 
then df,u « ]» + f» — <pl(v) {2 (ar — o) + 2 (« — c)p\ = 0, 
dyU -= 5^ + n - <tHp) {2 (y - 6) + 2 (« - c) g} = 0. 

Hence, transposing, and dividing one result by the other 
to eliminate 0'(t?), we find 

jp + w {m "- a) •¥{x-c)p 

or p {y — 6 - n («f - c)} — 9 {ar — a - w (« - c)} 

E3 n (^ * a) *- m (y ^ 6), . 

the partial difierential equation of which the complete inte- 
gral is 

X + mw + nysi <p {{/v - o)* + (y - 6)* + (« - cy}. 
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90. To integrate an equation in which (mly one of the 
differential coefficients of the first order enters with j% y, 
and «. 

Let the equation be 

Integrate it, considering y as a constant, and in place of 
the arbitrary constant C, add a function of y of arbitrary 
form. The resulting solution, containing one arbitrary func- 
tion, will have all the generality that can be attained. 

Ex. d,ir(«« + s^) ««* + »»; 



«« + y» «* + »• 



0; 



.*. tan"* tan~* — = tan"* ib (y), 

y y ^ 

or «-a? = (y* + ^»)0(y), 
<t>{y) being arbitrary in form. 

The equation / (<v, y, Zj 9) » is similarly integrated, 
the correction in this case being an arbitrary function of w. 

Ex. wydyX '\-nz = Oj 

91. To integrate the linear equation of the first order, 

P, Q, and R being functions of ^, y, and x. 

Let the primitive be P (^, y^ %)^0\ therefore, denoting 

^(#)^(^> y» ^) ^y '^\^)'> *°^ ^ ^'^ ^^'^ ^^® other co- 
efficients, we get 

P'(^) + P'(i»).p = 0, 

P'(y) + P'(«).g-0; 



Ill 

But dF(Wf y, z) = F'(w)dw + F'(y)dy-\' F'(x)dx^O; 

.'. PF" (a) d.v + PF' (y) dy + PF" (x) dz ^ Oj 

and PF(ai) dx + QF(y) dof + RF{x) d<r « ; 

A /^(y) {Pdy - Qd^} + F'Ciir) {Pdx - fida?} = 0, (l), 

which is satisfied by 

Pdy - Qda = 0, 



} 

0. j 



. (2)- 
Prf«f - Rdw 

Suppose that by integrating these equations, either se- 
parately or conjointly, we obtain M ^a^ N^by two relations 
between the three variables and the arbitrary constants a and 
6, which satisfy them. By means of these two equations, 
any two of the variables as y and x can be expressed in terms 
of a and b and the third variable w. The complete primitive 
becomes 

and the differential of (x, a, b) must by virtue of equations 
(2) be identically equal to zero, therefore (p {x^ a, b) cannot 
contain Xj and 

O = 0(a, 6) = 0(Jf, N), 

or M^fiN). 

The assumptions (2) satisfy equation (I) independently 
of the forms of F'(y)j ^(«), that is, independently of the 
form of jP; therefore the form of 0, and consequently of/, 
is arbitrary. Hence M^f{N) being a solution of the pro- 
posed equation and containing one arbitrary function, is the 
general primitive. 
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92. For .the success of the method, it is in general 
necessary that of the three equations 

Pdy - Qda m o, Pdx - Rdw = 0, Qdsg - Rdy = o, 

one at least should contain those two variables only whose 
increments or differentials appear in it. They are <;alled the 
reducing equations of Pp + Qq t^ R; and may be easily re- 
membered under the form 

dw dy dx 

T'^^'^'r' 

93. By a similar process, partial differential equations of 
three or more independent variables can sometimes be inte- 
grated. If jv be a function of a, y, t^ and 

Tn + Pp+ Qq^-R, 

where n a d^Xy we have 

R-Pp- Qq 

ns ^ -; 

T 

which, substituted in 

dx m ndt + pdw + qdy^ 

gives Tdx - Rdt » p {Tdw - Pdt) + q (Tdy - Qdt). 

- ... - . dw dy dx dt 

And if the equations -=r = "7r = ^;r=7;:5 or 

^ P Q R T^ 

Tdx - Rdt « 0, Tdx - Pdt « 0, Tdy - Qdt « 0, 

give L ^ a^ M tab, N t^Cj 

then the solution is 

0(L, Jf, JN0«O, or L^f{M,N). 

94. To explain the geometrical meaning of the solution 
of a partial differential equation of the first order. 
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If we regard w, y, x ba the co-ordinateB of a point, and 
the proposed differential equation as the equation to a system 
of surfaces, M=a, iV a 6, are the equations to two surfaces, 
which, being satisfied by values of /r, y, z which satisfy the 
differential equation, conjointly represent a line on a surface 
represented by the primitive. By giving any alterations to 
a and b we obtain other lines in space; but in order that 
these lines may always lie on the surface in question, these 
alterations must not be independent but connected, according 
to some law expressed by the relation 

a=/(6), or M=fiN). 

For example, the equation 

pw + qy ^» 

is the differential equation to a conical surface, and 

- = o, - « o, 

X X 

are the equations to any straight line through the vertex. 
If a and b undergo all possible variations consistently with 
the restriction imposed by the nature of the directrix, the 
assemblage of straight lines thus formed is the conical sur- 
face. 

We shall now apply the preceding method to examples. 

Ex. 1. pw + qx-^y^O^ 

dx^pdw -^qdy, 

mdx^pwdm -{-qady 

- ydw a pwdw + qxda ; 

.'. wdx + ydw « q {wdy - xdai) ; 

.-. wdx-^-ydw^Oj 

ady — xda « o, 

xdx + ydy =■ 0. 
16 
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From the last, V + y* = a*; 

.'. adz + s/c? —ffdw^ 0, 

dx dx 

+ — = 0; 



.'. sin''^ — + log^= logfc; 

Ex. 2. jp (^ - a) + 9 (2( - 6) B j2r - c. 

This compared with P|) + Qg = Jf, gives in the place of 

Pdx - Rda = 0, Qd» - jRdy = 0, 

the reducing equations 

{p '^a)dx^ Qs — c)dw ^ Oy (y - fc) dar - (» - c) dy = ; 

a?-a y-ft ^ 
a, = /3; 






«— c « — c 

« — c \«r-c/ 

Ex. 3. mp-^-nq^ 1; 

.'. md» — d^ «3 0, ndz - dy = ; 

.*. y — »«=/(a?- i»«). 

Ex. 4. (^ - mx) jp + (y - n») g = ; 

.•. (a - 1»«) d« B ; ,\ x^ Uy 

{w - m») dy - (y — niv) dor « 0, 
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or (w — ma) dy^(y^na)dwsaO; 



6; but a =s/(6), 



y — na 



.% X 






5. pw + qy^ nxy « = ^/ ( - ) . 

6. 5 + 1 = 1, «*..y+/(?). 

7. pa?«+gy«=«*, -«-+/( ). 

% a \y wj 

8. drd^x+ydyjv + ^OiJir + ttd^iVsni^, jjfss /•/(-, -, —I. 

10. (a^ + fty) p + (a'^ + }{y) q ^ c;ir. Assume z b e". 

95. The following examples require artifices in the com- 
bination of the reducing equations. 

Ex. 1. {y - fc - « (« - c)} p — {^ — o — m (» - c)| 5 

=s » (af • a) - HI (y — 6), 
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the equation to surfaces of revolution. Here the. reducing 
equations are 

{n (a? — a) — m (y — 6)} do; - {y — 6 — n (» — o)} d« = 0, (l). 

{y — &-«(«-«)} dy + {^-a-m(»-c)}dd7 = 0, (2). 

{a? - a — f» (jjf - c)} d» + {n (^ — a) - III (y - 6)} dy « 0^ (S). 

Multiply (l) by w — a^ and (3) by y — 6, and add^ 
.•. {n(af- a) — fi»(y — fc)} x 

{(a? - o) da? + (y - 6) dy -{-{%- c) d«\ = ; 
.•. {w -- a) dm -{- (jy ^ h) dy + {x -^ c) dz ^ 0% 
.•• (w - o)« + (y - 6)* + (« - c)* = o. 

Multiply (l) by m, and (3) by n, and add, 
.'. {n{w -^ti) - m (y — 6)} (iitdar + ndy + dx) « ; 

.*. mdx + ndy + div « o^ 

mo? + »y + » = /3 =/(a) ; 
.-. liT + mo? + ny «/{(a? - a)' + (y - 6)* + («f - c)*}. 

Ex. 2. (a? + y + iir) dt» + (^ -f y + iv) d«jv 

+ (^ + 0? + ») dyX SB ^ + tf? + y. 

The reducing equations are 

(a; + y + ir) dar e (^ + zr + y) d^, 
(a; + y + «) do? = (^ + y + «) d/, 
(ar + y + «) dy « (^ + 0? + «) d^; 
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•'• (a? + y + x) (dx - dt) «(<-«) d/, 

(d? + y + «) (dw - df) = (^ - /p) d^, 

(w + y + x) (dy - d^) « (^ - y) dt. 
Also, 
(ar + y + «) (dw -^^ dy -^ dx + d^) = S (f + or + y + «) d/; 

div — dt dt dw '\- dy •\- dx -{- dt 

t-x a-^y + X 3 (a? + y + » + ' 

••• (a? + y + » + t).(x - 0' = «> 
Similarly, (ar + y + ar + f).(^ - ^^ = /3, 

(or + y + » + 0-(y - 0^ = 7» 
.-. (w'^y + x + t).(x-ty 

= F {(or + y + « + 0.(^ - 0^ (^ + y + » + 0-(» - 0'}* 

96. To integrate F (^, y, iv, p, q) » 0, when it contains 
terms of more than one dimension in p and q. 

In order that 

dx =pdw + qdy 

may be a perfect differential, we must have 

^(y)P + 9d^z)P = d(*)? + P^w9' (2). 

If from the proposed equation we determine 

equation (2) becomes an equation for the determination of p 
of the form 
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the iategradoa of wHidi depends on the integration of one 
of the equations 

dp — Ldm s 0, dp-^ Mdm » 0, dx — Nd* « o. 

Let p ^f(^9 99 *9 o) be found from these equations^ a 
being an arUtrarj constant This value of p and the onrres- 
ponding value of q found from the proposed equaticn, being 
substituted in dssw^pdm + qdy render it an exact differential; 
and thus a yalue of s will be obtained involving two arbi- 
trary constants a and 6 ; and this will consequently be the 
complete primitive. The general primitive may be obtained 
by putting h^y\f (a), differentiating the equation with regard 
to a, and eliminating a. The result containing one arbitrary 
function is as general as any solution which the equation 
admits. 

Ex, 1. p* + 9*«l, 

dp -iO, 

p 

zdy — da^ 0, 



x/i--* 



P 



n dw 

^ dx J »0; 



\/i-p" \/T^ 

.'. p-a, q « %/l - o% 



jif - oar + v^l-a*y + 6, 
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or » as aof + V 1 — a*y + (a), 

ay 



a?- 



vT^ 



H-d„^(a). 



If a be eliminated between the two latter equations, the 
general primitive results. 

This is a simple case of the problem of finding surfaces 
of equivalent area to a given surface, that is, so that any 
cylindrical surface parallel to axis of x may always intercept 
equal areas in the required and given surfaces. 

If P, Q, be the values of d,», dyX in the given surface, 
the general condition is 

which, if the given surface be a plane, leads to the equation 
of this example. 

Ex. 2. X = pq^ 

« 



y + o 



^ + 0(a), 






Equations of the Second and Higher Orders, 

97- Here, besides the coe£Scients p and q of the first 
order which may be involved together with w, y, and jjf, the 
equation must contain one or more of the coefficients of the 
second order r, «, ^; so that in its most general form it 
will be 

F(a, y, Xj p, q, r, «, t) = 0. 

In partial differential equations of the second order, we 
cannot be certain of the; form of the solution^ nor pronounce 



190 

befiirdiaiid how many ar li itntf j fooctioiiB H on^it to oontam. 
For let 

«»/{ir, y, jr, 0(«), ^(»)} =0, 

be an equation cootainmg two arbitrary fanctioo. of v «id 
w two known functions of jr, y, and m ; then p and q will be 
giyen by die equations 

and r, #, and / by the equations 

These together with f» = make six equations, into which the 
NX quantities 

0(r), ^(«7), 0'(r), ^'(«7), ^"(r), ^"(«^), 

may enter. Consequently, it will not generally be possible to 
eliminate these six quantities and obtain a relation between 
^9 99 ^9 P9 99 ^9 ^9 ^9 independent o( the forms of the func- 
tions <p and yjf ; although particular cases do occur in which 
this can be effected. 

In general, it u^O contain n independoit Tariables, and 
m functions of the form <p (o), where o is a determinate func- 
tion of the yariables; these functions can certainly be eli- 
minated between the equations obtained after (r) differentia- 
tions, if 

1 in , ^(^ + ^) , . n(n + l)...{n + T.'l) ^ _^^ ^ ^^ 

1 + w^ -f -— T — + ••• H i > i» Ir + 11, 

1.2 Ir 

the former expressing the number of equations, and the latter 
the number of functions. 

98. As an example of forming a partial differential equa- 
tion of the second order, by the elimination of two arbitrary 
functions, we may take 

u = y -- w(l) («) -^ yff («) := ; 

.'. dgU^ -<p{%)- (0^\«)p^ y\/{z)p^Oj 
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•'• P + ^^(«) = 0; 

or ^r - 2pq8 + p*/ = 0. 

99. The equations 

rf> + Pd,« + Q = 0, 
d^» + Pdy» + Q = 0, 

where P and Q are functions of w^ y, and x, must be in- 
tegrated as equations between two variables, y being regarded 
as constant in the former, and w in the latter ; and arbitrary 
functions of those variables respectively, being introduced 
instead of constants. 

100. The equations 

dgdyX + Pd^x = Q, 
d,dyX + PdyX = Q, 

where P and Q do not contain x^ are reducible to the case 
of Art. 90, by considering d^x or dyX respectively as a single 
quantity v. 



Ex. dgdyX + ; — 7jd,» = ay*. 






This being a linear equation in d^x which is made in- 
tegrable by the factor . , we have 

16 
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101. Similarly, equations of the forms 

/(a?, y, isr, d,«, cgjjr, &c., dlx) = 0, 
/(^> y> «> dy«» dySB, &c., djar) « 0, 

may be treated as if they contained only two variables; 
arbitrary functions of y instead of constants being introduced 
into the solution of the first, and arbitrary functions of a? 
into the solution of the second. To this case may also be 
reduced the equation 

/(^, y. ^«, d^dlx, dld;z, &c., d^d;z)^0; 

for by putting d^x = «, it becomes 

/(^> y> «5 d^v, dlv, &c., c^tj) = 0, 

which will give a value of v containing m arbitrary functions 
of y ; and then d*x = v will give x involving n arbitrary 
functions of <r. 

102. To integrate the linear equation of the second order, 

Rr + Ss -^ Tt ^ F, 
where J?, *y, T, V are functions of a?, y, «, />, g. 

By means of the relations 

dp = rdiV -h «dy, 
dq ^ sdw + ^d^ ; 

eliminating two of the three coefficients, r and /, from the 
proposed, we get 

Rdpdy + Tdqdx - Fda7dy = « {i? (dy)* + ^(dci?)* - Sdwdy}^ 
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which is satisfied by 



Rdpdy + Tdqdw - Vdxdy = 
R {dyf + T{dxy - Sdwdy = 



(1) 



Let M = a^ JV = 6, be two relations between w^ y, %, p, q^ 
and the arbitrary constants a, 6, which satisfy these equations ; 
then M = (p {N) satisfies the proposed equation. This will 
be shewn by proving that it can reproduce the proposed 
equation. 

Let dy = mdw ; .-. Rm^ - Sm +7^ = 0. 
For each root of this equation, we have 

dy - mdof = 0, Rmdp + Tdq - Vmdx «= ; 



.". dy = mdof^ 

Vm , Rm , 
dq = — - dx — — tfp, 



dsis = pdoB H- ^dy. 



(2), 



Hence M ^ a gives on differentiation 

= dgM .djff + dyM .mdw + d^M {pdoo + qntdw) 

+ dpM . dp + d^M A — dw - -jT^p] » 

wherein all the known relations (2) having been introduced, 
dw and (2j> must be independent, 

Vm 
.-. = d^M + mdyM + d^M . {p + wg) + -=- d^M, 



Rm , ,^ 

= dpM — =- dqMi 
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Vvn 

.-. d^M «= - [mdyM + (jp + mq) d^M + -=- dqM], 

Vm . 
So d,A''=- {mdyAr+ (p + wig)d;,JNr+ -^a,iV}, 

d^N^^d.N. 

By differentiating the assumed equation Jf = (iST) we 
have 

Now 

Vvn 

dM=^- {mdj^M+ (p + mq)dgM -{' -z=^dqM}dw + dyMdy 



Rfn 
+ dgM (pdof + qdy) + -— d^M.dp + d^M.dq 



(dyM + qdgM) (dy -mdw) + -^=-(Rmdp + Tdq- Vmdai), 



and a similar value exists for dJV^; 

.'. (dj^M+qdgM) (dy- mdai) ^-^—(Rmdp + Tdq-Vrndw) 

= 0'(JV) {(dyJVT + qd,N) (dy - m(^<r) 

d N 



which may be put under the form 

Rmdp + Tdq - Vmdw = co (dy - md^) 
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or Rm(rdw + sdy) + T {sdw + tdy)- Fmda?= ^(dy-wda?), 
where dx and dy are independent; 

.*. Rmr + Ta - Vm = — ww, 

Rms + Tt = w; 

.-. flr-f 5'«+ r/== F--(i?w«-*yw+ r)= V. 

m 

Hence, M = <j> (N) satisfies the proposed equation. 

According as the roots of 

Rm?-'Sm-\-T = 

are unequal or equal, we are thus supplied with a total or 
partial differential equation for the determination of %. 

Obs. As the reducing equations (l) may contain 

^» y> ^» P^ Q9 

and as these together with dz » pdw + qdy will generally 
lead to an equation containing three variables, which will 
not always admit of a single primitive (Art. 83.), it may 
happen that the first integral of the proposed equation 
cannot be determined; but we must not thence conclude 
that the proposed equation does not admit of being solved* 

103. Hence, to integrate the linear equation of the 
second order 

Rr'^S8+Tt= F, 

the process is to obtain a value of m from the equation 

to substitute it in the system 

dy — mdof = 0^ Rmdp -{• Tdq^ Vmdw^ (l.) 
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to satisfy these, conjointly or separately, by two relations 
between .», y, ^, p, q^ 

then to put M ^ ff> (JV), and to integrate this equation of 
the first order. 

104. If R, Sj T9 be constant, and V a function of ,v 
and y only, then the values of m will be numerical, m and 
n suppose; and the integrals of equations (l) will be 

y - ma} = a, Rmp + Tq ^ m f^V -^ b^ 

where, previous to integration, mw + a is substituted for 
y in F, and after integration the value of a, viz. y — mof, is 
restored: consequently calling this value Fi, since Rmn — T, 
we have 

Rp + Rnq = Fi + (f>(y - mw). 

Next to integrate this equation of the first order, we have 
the reducing equations 

Rd% - {Fi + <f>(y - mai)\ dci? = 0; 
.*. y — nw = o, 

nw + a being substituted for y before the integration is per- 
formed, and afterwards the value of a, y — nx^ restored ; 
this gives jiFj= Fg, suppose, 

and /,0'(y - mx) = /.^' {(w - m) a? + a} = — — - 0(y - ma?) ; 

hence, including the constant multiplier under the sign of the 
function, 

Rz =^ Fg + (y - mx) + ^ (y - nx). 

Ex.1. dl% = a^dl%^ or r - a^/ = 0, 

or rdxdy - a^tdydx = 0, 
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or dpdy - s{dyf - a^ \dqdw - s{dxy} = ; 

.-. (dyy^a^'idwy, 
dpdy = a^dqdw^ 
dy = ^ adw. 

I . dy = adx^ y — aw == a^ 

dp^adq, p - aq ^ fi=f(y - aw), (l). 

2"^ dy = -- adWf y + aw ^ a, 

dp^" adq^ p-\-aq = fi = ^'{y + aw). (2). 
By adding and subtracting (l) and (2), 

2p =/'(y - cb^) + (i>{y + aw)y 

2aq=:: (p'{y + aw) -f\y - aw) ; 

^•. 2ad^ = 2apdw + ^aqdy = 0'(y + aw) (dy + ado?) 

-/'(y - «*^) (dy - adw) ; 
.-. ar = 0(y + ao?) +/(y - ao?). 

_ 2p p* 

Ex. 2. r --^« +~ / = 0, 

q q^ 

9,p T? 

rdwdy sdwdy + ~ tdydw = 0, 

qr q^ '' 

2© ©^ 

dpdy - 8 {dyf sdwdy + — {dgd<r - * (d<r)*} = ; 

.-. (dy)2 + 2 ^ dydo? + ^ (dci?)2 = 0, 

J? 
dpdy + — dgdiX^ = ; 

T 
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P 

q 

P 

dp - - dg = 0, 
9 



dz = pdw -k- qdy = , .•. » = a, 

qdp-pdq p 
5 = 0, - = 6 ; 

9 9 

To integrate this equation of the first order, we have 

dz = 0y or z = a, 
dy-^f{z)dw^O, or y -^0^/(0) = fi; 

this is the equation to the surface generated by a straight line, 
subjected to pass through three given fixed curves. 

Ex. 3. r — a^t = a?y, 

» = ^ o^y + f(y + ax) + (y - ao?). 

4. j;* r - y*# = 0. iif = \/^. ( -^ ) + "^ (^y)- 

6. o^r + 2fl7y« + y*^ •= 0, «f = ai(f> f - ) + ^ ( - ) • 

7. r H- (a + 6) « + aht = .jyy, 

« - ^ ar^y + —• (a + 6)«i?* =i (y - a^) + x|, (y - 60?). 

8. (l+p^ + g*)rH-«(9*-jp*) -(1 +pg+p^)^ = 0. 
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CALCULUS OF FINITE DIFFERENCES. 



SECTION L 

DIRECT METHOD OF DIFFERENCES. 



Definitions and Priuciplea. 

Art. 1. In the Differential Calculus it is the first term 
only of the series, arranged according to ascending powers 
of A, expressing «,^^-w, (where u, is any function of w 
and u^+f^ the same function of ^ + A), or rather the coe£Bcient 
of h in that term, with which we are principally concerned, 
and which we usually write hdgUg.. But in the Calculus 
of Finite Differences, it is the whcde of that series which 
forms the object of our investigations, and it is usually written 
At^j,, so that 

2. It is common, however, to suppose the finite incre- 
ment which the principal variable w receives, to be 1 instead 
of A, both for the sake of simplicity, and because that is the 
value of the increment when u^ is regarded as the general 
term of a series; and in that light it is by far the most 
frequently regarded in finite differences; so th'kt 

Aw, ««<^,+i -w,. 



There are, however, in this subject, several important 
theorems which it is advantageous to investigate on the hypo- 
thesis of an indeterminate increment h, for the principal vari- 
able, instead of unity, as the process is the same on either 
supposition, and the result one of greater generality. And 
if in other cases it should be desirable to introduce the same 
hypothesis, the expression must be prepared by first writing 
h» instead of o^, theii performing the operation on thie usual 

supposition o{ Ax^ 1, and in the result writing - for «. 

h 

3. By a Series is meant a regular progression of terms 
increasing or decreasing in magnitude according to a certain 
law; hence, when that law is given, and also the place of 
any term in the series, the magnitude of the term may be 
found, and thus the successive terms of the series may be 
produced in order. The place of any term in a series is 
assigned by giving the number of terms by which it is re- 
moved from some one which is considered as fixed. This 
nuknb^r is called the index of the term to which it belongs^. 
Thus in the series 

0, Ij 8i 27, 64, ... a?^... 

taking the first tenii as the point of departure, we have the 
ebrresponding series of indices 

0, 1, 2, 3, 4, ... a?, ... 

If the series be continued backwards, the indices must 
be considered as negative; thus the backward continuation 
df the above series gives the terms 

••• — dr, ••• — 27, —8,-1, 

with the corresponding indices 



w^^mtmti^^K^mmmmmmmmmmmmmm^mmm^^mmmmimmmm^mmmm 



4. Since the magnitude of every term is determined 
solely by its index and by the law of the series, it follows 
that any term is a certain function of its index, the form 
of which does not alter in passing from one term to another, 
but remains the same throughout the whole series. Thus 
in the above series every term is the cube of its index. 

This function analytically expressed is called the general 
term of the series; (in the above series the generiJ term 
is ^;) and it is evident that all the t^ms of the series 
will be produced from it in order, by substUutipg successively 
for the index Wy the progression of natural numbers 

The general term pf a aeries is usually deppt^d by u«, 
where u^ is a certain function pf m determined by the nature 
of the series. Thus, u, denoting the general term, the series 
will be 

5. The excess of any term Ug.^i above that which immedi- 
ately precedes it, or the function u,+i -t^« is called, as has 
been stated, the Difference of the function u^^ and is denoted 
by AUg. (In certain cases, which will however be expressly 
mentioned, we shall take Au, to mean t*,^^ --«*,•) 

It is obvious that u^^i - u^, is itself in general a certain 
function of w, the nature of which is entirely dependent on 
that of the original function t/«, from which it is derived, 
and susceptible of a difference. 

The difference, consequently, of the function Au, (which 
must be considered as having Au for its charaoteristio, in 
the same manner as u^ has u) is 

A (Aw,) = Au^^i - AUg, 

which is usually written A*ff«» 



tn like manner 

A (A*«.) = A»«. « A*«,„ - A»«., 

A (A»«,) » A*«, - A»«.+, - A»«., 

A"«, = A— •«,+, - A— •«,. 

6. Hence if in any function of # we change # into x + I^ 
and from the result subtract the proposed function, we obtain 
the first difference of the proposed function ^ and the second^ 
third, ice. differences are formed, each from the preceding, 
by a similar operation. To determine these differences of 
given functions, and to investigate the relations which hold 
between differences of any orders and the functions firom 
which they are derived, is the object of the direct method 
of differences. 

We shall now proceed to give instances of finding the 
differences of various functions, according to the above defi- 
nition. 

Differences of Explicit Functions. 

7. To find the difference of aUj, + c, u^ being any func- 
tion of Wf and a and c quantities independent of w. 

A (au, + c) ■ au,^i + c - (au, + c) =* a («**+i - w,) = a Ati,- 
Hence, making a^Oj AesO. 

8. To find the difference of the sum of any number 
of functions of a. 

A {u^ + 1?^ + w^) « w,+i + f?,+, + fir,+, - (u^ + ©^ + w,) 

*= Atf«4- A«^+ Aw,. 



9. To find the difference of the product of two functions. 

« «,Ai?, 4- t>,A«, + Aw,. At?, « w,At?, + »,+i Aw,. 

10. To find the difference of the quotient of two func- 
tions. 






v*+it?* 



11. To find the difference of the continued product of 
any number of successive values of a function. 

Hence in the particular case where w, « a + fto? 

A («,W,+ i ... U^+n) =* W*+l«*+2 ••• U»+n • (» + 1) 6. 

12. To find the difference of a fraction whose numerator 
and denominator are the continued products of any number 
of successive values of two functions w, and v, respectively. 

A « {Vs^s+n+l - W,1J,+^+i), 

Hence A ^ ■* — » 

and in the particular case where %^a + bwy 

1 (m+l)b 

^_ -r ^ 
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13. To find the differenoes at any ratkHial integral func- 
tion, and to shew that the n^ difference of a rational intend 
function of the n^ degree, is constant. 

Let «,« J«* + i'^~^ + — +Jrx-i-£, be a rational intq;ral 
function; then its first difference is 

A«,«^{(ar+l)"-«*}+J?{(jr+l)-*-af»-'} + ... + /(2df+l)+X 

which is a rational int^;ral function, one d^ree lower than 
the original function. In like manner, for the difference of 
this, or the second difference of u^j we have 

A*w, = n (n - 1) Jjr-« + B^af"^ + ... + /„ 

and so on; and for the n^ difference, we have 

A"ttr = n (» - 1) ... S . 2 . 1 . J. 

Hence the n^^ difference is constant, and the differences 
of all orders superior to the n^ vanish. 

Also A" (/»*) « 1 . 2 . 8 ... n. 

14. To find the differences of a'. 

^*a' = (a - 1) Aa* = a' (o - 1)% 
A"a' - a* (a - 1)". 
Also A a^' " a'''+^'^' - a?^' » a»- {a^^' - i). 

15. To find the difference of lpgo«. 

^ Oog^'*) - log^'.+i - log*'^ - log^^ - log f 1 + — !^V 



16. To find the differences of sini?^ and cost?,. 

• /' AX- ^ • ^^« / ^M 
A sm Vg = sm (v^ + Av,) — sm «, = 2 sm cos It?, + | 

= 2 sin ' sin {t?, + ^ (tt + At?,)}. 

A /AX o • ^^» • / ^M 

A cost?, e cos(t>, + At?,) - cost?, s» -2 sin—— sm It?, + 1 

= 2 sin — ' cos {t?, + 1 (tt + At?J}. 

2 

2sin-l sin {a76 + a + -(7r + d)}, 

2sin-| cos{a70 + a + -(ir+d)}. 



17- To find the differences of tan t>, and tan"^t?,. 

sin t?^+i cost?, - co8t?,a.i sin t?, 

Atant?^«tant?,.i - tant?, = ^^^ ^^ 

cos t?, 4.1 bos t?, 

sin(t?,+i — t?,) sin At?, 

cos t?,+ 1 cos t?, cos t?,^. 1 cos t?, * 

sin d 



Hence A tan ofO 



cos (<!P + 1) d cos w0 ' 



t?-xi — t?. 



Also Atan-^t?, «tan"^t?,+i - tan"*t?, = tan*^ -^^^^ — -^ 

' An 

= tan"* 



At?, 



Hence A tan"* a^O = tan"* 







1+ (a? + 1) 
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the MWCCMiTe Vaioe* and the Di g i ii i«i df a 



18. To expreM A*tf, bj «, and its it socoosTe TaliMiy 



*^#^* > ***+t* » ••* **#+•» • 



We here take h ioctead of unit j fiv the increment of the 
p ri ucip al Tariable, as the inTestigatioo is precisely the same 
on either supposition. 

Now suppose this hiw of the coefficients, which as far as 
we have gone is the same sb that of an expanded binomial 
whose index is the order of the difference, to hold for the 
n^ difference^ so that 

A"t*, - 
then A"'*''w, • 

- K + nA -Pl«** + (n-l)* + — ^ Pi^s + Sk T Pl^M + k =^ «**) 



r*^ 



which is the same alteration with regard to the coefficients as 
occurs in passing from (x - l)" to (« - l)"+\ If therefore, 
for any value of n supposed a positive integer, the coefficients 
of the expansions of A"w, and (x — l)" are the same, they will 
always be the same ; but these coefficients are identical when 
n = 1, 2, 3 ; therefore they are always the same ; 

or, supposing A « 1, 
A« n(n-l) 

Obs. For the sake of being easily remembered, this 
result may be written 

A"w, = (w, - 1)"; 

where, instead of the powers of u, from n to 0, the n + 1 
descending values ««,+«, t^*+»-i5 ... u^+i^ Ug are to be written. 

19* If in the series just investigated, we assign a par- 
ticular value to Ug^ we shall readily obtain an expression for 
its rf^ difference. Thus let u^^af*^ 

.\ A*(/p*) = (a? + w)'"-»(a? + n-l)"* + -^^^^^(<r+n-2)"-^c., 

in which equation, if n>m, since the former member vanishes, 
the second member is zero for every value of ^; and if m^n, 
so that A"^ = 1 . 2 . S ... n, we get 

1.2.S...n = (a7 + w)"-w(^ + n-lV+ -^^ ^(/p + w-2)"-&c.; 

\/v / 1,2 • 

and making w saO^ since the equation holds for all values of /r, 

w (w - 1) , 
l-2.S...n = w"-»(n-l)"+— ^ ^(»-2)"-&c. 

2 
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t If MmO^ Mad A*CP doHite the |Wfftifnkr value of 
A*ii^ vlieii # « 0, ire li«Te 

Of the nanben eompriaed in the Conn A*0^, we shaU 
make comideraUe tue in fatme iiiTestigatioos ; whenever 
n>m the ralue is zero, in other eases it may be computed 
bjr the above formuk; thus, 

A(r-1, 

AV- 2, A*0^ - 6, AV = 14, ... 

A*0^-6, A*0*-56, A*o»=l50, ... 

81. Reversing the order of the series in Art. 18, we 
find 

(- !)• A*«^ - t*^ - «t^^4.i + , . t*.+f - .-. * «,+.- 

Alio putting Ug^of^^ and then supposing msn, jrsl, 
we find successively, 

(-l)"A"/i^ -^-»(j?+l)«+ ^^^ ^ (ar+2)*-... ±(j?-i-n)"', 

1 • z 

1*2 

(-l)**L2.S..n- l"-n.2» + -^^ ^ . 8* - ... ± (n + l)*. 

1 «2 



t. Hence it may be proved that l . 2 . 8 ... (jp - i) -i- 1 
is divisible by p^ if p be a prime number. 

Let I*- ■■ «" - 1 ; 



i 
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.-. A"(cv"-1)= l.2.3...w«(^ + »)*-l-«{(ar+n-l)*-l} 

n(n ~1) f ^ , ' 

+ — ^ |(^ -^. w - 2)» - 1} - &c. 

Let X = 1, and n + I <= p, then a> + n^p, and 
1 . 2. 3 ... (p - 1) + 1 = p» - (p - l){(p - 1)'-' - 1} 



1 .2 



Now by Fermat'^s Theorem every term of the second 
member is divisible by p when p is a prime number; con- 
sequently 1 . 2 . S ... (p - 1) + 1 is divisible by p. 

23. To express t^,^„;^ by u, and its first n differences. 

We take h instead of unity for the increment of the 
principle variable, the investigation being precisely the same 
on either supposition. 

w,+3^=w,+2Aw,+ A*w,+A(w,+2A«. + A*w,) 
«w,+SAw, + SA"w,+ A^w,. 

Now suppose this law of coefficients, which as far as we 
have gone is the same as that of an expanded binomial, 
whose index is the number of increments which the principal 
variable has received, to hold for n increments, so that 

^*+nh = w,+PiA«.+P2A*w,+ ... +Pi A»-*«,+ A"w., 
then tt,+(^+i);k = w*+PiAw,+l>2A*w,+ ...+Pi A"~*w,+A"w, 

+ A(w*+PiA«,+ ...+P2A""*w+jPiA*~*tt,+A"w,) 
= w, + +Pi) Aw,+ (pi +Ps) A^w, + ... 
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whidi is the same alteration with regard to the ooeffidents 
at oocun in passing from (l + %y to (l + %y^^. 

Hence if the coeflBdents of the expansions of u^^^^ and 
(l + xy are the same for any yalue of it, supposed a positive 
integer, thej will always be the same ; but they are identical 
when n * 1, 2, 3 ; therefore they are always the same; 

or, supposing A"b i, 

ft Ci% ^ 1 ^ 

Obs. This result in order to be more easily remembered 
may be written 

each term of the development of (l + A)* being understood 
to be prefixed to u^. 

24. Let Ug as ^, then w,^, =■ (a? + n)" ; 

.•- (4r + n)"»(l + A)*^, 

and making jr « 0, n** = (l + A)'(r, 

each term of the development of (l + A)* being prefixed, 
as said above, to jr"* and 0"*, respectively. By the latter for- 
mula any power of a number is expressed by the numbers 
comprised in the form A"0". 

25. To deduce TaylorV Theorem from the formula 

n (« - 1) ^ , 
«**+n* - «*# + » Ati, + A*w, + &c. 

Let nhmtf and let h be infinitely diminished whilst t 
remains finite; therefore n is infinitely increased; and since 
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h is indefinitely diminished, we have, regarding the differ- 
ential coefficient as the limit of the ratio of the simultaneous 
increments of the function and the variable, 

A*w^ , /Aw 



= dgl~—^\^dlug; &c. 



Hence, preparing the formula as follows, 

Aw, nhinh-h) A^Wj 

IT TTi h^ 

and taking the limit of both sides by supposing h to be 
infinitely diminished and n infinitely increased, their product 
always remaining equal to a finite magnitude t^ we get 



u^+nh ^u,-\-nh — — + — zj- + &c. 



u»+t " ^« + td^Ug + dj Ug + &c. 

1 . 2 

The Differential Calculus is a particular case of that 
of Finite Differences; and the above investigation is. intro- 
duced to shew how, from results in Finite Differences ob- 
tained with an indeterminate increment for the principal 
variable, we may pass to the corresponding results in the 
Differential Calculus. 

26. The theorems of Arts. 18 and 23 have been proved 
by an inductive process; they may also be established by 
the theory of Generating Functions, the principles of which 
we shall now proceed to explain ; as it is a theory which, for 
its generality and power, especially merits our attention. 

Generating Functions. 

27* Let (f>{t) be a function of t susceptible of the de- 
velopment 

(pit) = ... + W.ji?"* + Wo + tti^ + ... 

+ w,_i^'"^ + uj" + u„+it"^^+ ... ; 

then u, may evidently represent any function of a? whatever, 
if we regard this equation as the definition of (f>{t). The 
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function ^(0 consequently by its devekqpment generates the 
coefficients Uo9 t«i9 *•• Ug annexed to their proper powers of /, 
and is therefore called the Generating Function fd ««, and 
is denoted by Gru,, so that 

<l>(t) « Gug. 
Thus since 

log (1 -/)-» =^ + A/* + 4/*+ ... + i /'+ ... 

log(l-0-' = Gi. 

Similarly, since 

t(l -/)"* = ^ + 2^ + 3/*+ ... +ar^'+ ... 

28. To determine the generating functions of Ug+^ and 
A*Ug from that of Uj,. 

Let 0(0 s Gt^«, 

then 0(0 « ... + u^f + «,+ir"*^* + ... + tt*+«^'*^*+ •.. 

.-. t"'^{t) « ... + w« + w«+i^ + ... + w»+.^ + ... 

/. ^""0(0 - G^w#+»i5 or t'^Gug^Gug^^. 

Again, f*0(O - ••. + «-« + «*-ii+i^ + ••• + «^*-«^ + ..• 
.•. ^0(0 =» Gtt,-«, or ^G«, = Gw,.„. 

Hence it follows that the g^ierating function of 



Aw« or«,^., -«^ 



«9 



is g - l) 0(0, 



for this function being developed will produce the difference of 
two series whose general terms are respectively Ugj^itf and Ugif\ 
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.-. G(Aw,)= (^-l)(?«,. 



Similarly, 



and G (A*tt,) = (7 - l) Gu^- 



Also 



GCA'-w..,)- (7-1)"^^'-"= (^-l)VGt^,= (l-0"Gw,. 



29' To investigate the expression for A"t^« in terms of 
^tf-i-fif ^«+ii-i9 &c., by Generating Functions. 



G(A-t*.)-Q-i)*Gw. 



= r»Gw, - «r »+'(?«, + ^^^ ^^ r"+«Gtt, - &c. 

'1.2 



= G^M+n - WGW,+„.i + -— 6rW,+,.g - &c. 



_ , n (n - 1) ^ ^ 

G (w,+» - »«,+».! + — — - — w,+«-8 - &c.) 

1 • z 



1 • z 



for, the generating functions of both being the same, the co- 
efficients of f in their developments must be identical, however 
those developments have been effected. 
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30. To inyestigate the expresskm taar u,^^ in tenns of u, 
and its first n differences, bj Generating Functions. 

- Gu, + n(;(A«.) + " ^" ~ ^^ G (A*«,) + &c. 

1.2 



G (tt, + n Au, + A*tt, + &c.) ; 

1.2 



» (n - 1) ^ , 
.% «,+«- «^ + n A«, + DkUg + &c. 

1*2 



31. It is obvious that by transforming the expressions 



[- - 1 j <tW„ and t-*Gu^y 



in different ways, we may obtain various other expressions 
for A^Ug and u,^^ besides the above. 

Thus to express A^u, in terms of A"w,.., A'+^w^.,,-!, 
&c., we have 



G(A*w,)- (--1) Gu,^^ — ~^ — -4 



(1 - tyGu, + n(l - ty^'Gus 
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+ !?^Lli} G (A«+*w....,) + &c. 



1.2 

n(« + 1) 

+ 1 

.-. A"w* = 



A-+*«,.,.« + &c| ; 



n(n 4- iV 
A-w,., + n A«+^t^,.„.i + \ ^ ^ A-+'f^,.,., + &c. 

e 

32. Again, to express a,^.^ in terms of w„ Aw,_r» 
A*«*,-ar5 &c., we must transform #"" into a series of powers 

of r ( — 1] , or develope t~^ in powers of m from the 

equation 



1 



•ar- 



which may be done by Lagrange^s Theorem ; and we find 

w(n + 2r-l) 

«*»+« *» «* + » At«,.^ 4- 7-r A'^a^^tr 

1 « « 

» (fi + 3r ~ 1) (/n- Sr - 2) , _ 

+ -— ^^ — A't*^^3r + &c. 

1.2*3 

(HerscheFs Examples). 

Obs. This method is obvioudy not confined to the 
function t^,^.i - w,; it is equally applicable to any other com- 
bination of the successive values u^^ t^«+i9 ^^+29 &c*9 ^f ^^^ 
first degree. If we take Aw, to mean aw^+g + 6w,+i + c«„ 
then the generating functions of Ae^« and A"s«« will evidently 

and the expression for A"t^« in terms of Ug and its successive 
values, might be obtained as in the preceding case. 
3 
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Sepantion of the SymboU of Operation from those of Quantity. 
33. We have seen (Art. 23) in the formuki 

an instance of the system of notation, which consists in sepa- 
rating the symbols of operation from those of quantity ; the 
use of which is not confined to simple cases like that just 
noticed, but may be extended with remarkable effect to a 
great variety of investigations connected with this subject. 

If the expression (l + A)", regarded as a functiou of a 
certain symbol Af be expanded in powers of Zi, it will pro- 
duce the series, 

n (n - 1) . 

1 -f n A + — ^^ ^^ A* -f &c, 

1.2 

and if the symbol A be looked upon merely as an instrument 
by means of which we are enabled to produce the numerical 
coefficients of the series affected with their proper powers 
of A, then the expression (1 + A)" must be considered as 
having no other meaning than as an abbreviated expression 
for its development; and when prefixed to the function u,, 
each term of this development is understood to be applied 
separately to that function ; to that the following expressions, 

(1+.A)»M„ 

(1 + n A + — ^ ^ A'' + ... ) «„ 

w(n-i) 
n, + fiAu, + — ^^ ~ A*u^ + ... , 

JL a Z 

may be used indifferently for one another, the two former 
being regarded merely as abbreviations of the latter. 

34. In general, if F (A) be a function of A capable of 
being developed in a series of powers of A, such as 

J A« + B A^ + &c. ; then for AA^'u^ + BA^u, + &c. 



6 
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the expression F(A)Uj, is used as an abbreTiation ; and the 
same notation is applicable to other characteristic letters, such 
as d,j jT, 2, Sg9 &c. 

Hence also, the successive performance of two or more 
series of operations represented by F (A), /" (A), &c. upon 
the same function t^^, is equivalent to the performance of 
that series of operations which is represented by their product. 

For example, we have seen that 

A (u^v,) « Aw, . <>,+ UgAv,^i « (A + Av) w,f?^ 

supposing A to affect u, only, and using At? for an operation 
affecting v, only, such that (A »)'«, = A'v,4.r; 

then A« (w.t?,) « ( A + A t>) A (u^v^) = ( A + A «)*«*,«. ; /^ I^a/z 
and, generally, < i 

A-(t^,t?,) = (A + At>)>.t?, 

_ ..) : 

1.2 - — ^ 

which may be also proved inductively by shewing, as in 
Art. 18, that the coefficients of the developments of A'*(t^,t?,) 
and (l + «?)*•, which are identical when » = 1, undergo the 
same changes in passing from » to » + K 

Using an accent not to imply that the operation denoted by 
A is altered at all, but merely that A' affects v« only, whilst 
A affects u^ only, the above results are sometimes written, 

A (w,t>,) = (A + A' + A A') w,t>,, 

A^ Kt^x) = (A + A' + AAju^v^, 

and, generally, A" (t^,t?,) *= (A + A' + A A')" «<,«, ; 

or, if there were more functions te^,, z^^ &c. ,i|nd we use 
A", A'", &c. to imply that these symbols only affect ic,, fs^^ &c. 
respectively, we have 

A* («,«,«?,«, ..,) 
= {(1 + A)(l + A') (1 + A")(l + A'") ... - I }" w,tj,ii;,*, ... 



n (» — l) 
A"«, • «?# + » A""*w, A«,4.i + A*"'*^w,A*v,+, + &c. 
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35. To shew that A"u, ^ (^' - l)"tt«9 in which the syn- 
bols of operation are separated from thpse of quantity. 

By Tayl<Nr*s Theorem, we have 

«#+• « «, + nd,u, + -— - cCtt, + &c., 

1 .2 

and separating the symbols of operation from those of quan- 
tity, we get 

1.2 1.2.3 ^ 

Similarly, u.^,.,-i^-^>'-u,. 



ti,+i-ee''i*,; 



1.2 

or, again separating the symbols of operation from those of 
quantity, 

A"tt. = («^' - ««^— ^>^' + ^ ^^ " ^^ «<— *>^- ^ &c.) w, 

1.2 

a celebrated theorem first given by Lagrange. 

36. To find a general expression for A*Ug in terms of 
u, and its differential coefficients. 

The development of the second member of the equation 

.A*w, = (c^-- l)X 
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will consist of a series of terms of the form 

or AoUg + AidgfUg + ... + A^^Ug + ... > 

and ^„ is evidently the coefficient of f^ in the expansion 
of (c« - i)\ 

Now (e* - !)• = €"* - W€(-^>« + ^^^"^^ e<-«>* - &c., 
^ ^ 1.2 

and the coefficients of f" in the developments of e**, c^""^)*, &c. 
are respectively 

w* (n - 1)*" 



&c.; 



• . ^x« 



7» (» -^ 1) 



.~ {»"• - n (n - 1)- + '^ ^ ■ (w - 2)- - &c. } 



A"0* 



(Art. 20). 



Now so long as m <ny this vanishes ; and when n^m, 
A"0"« |n; 

A"0"+* A*0"+* 



» + 1 



W + 2 



Ex. Let Ug c= 0?"*, then 
A*(0 « m(w - 1) ... (m - w 4- 1) 0?"* 



— « 



nna + l 



A"0 



n-\-l 



. m (m - 1) ... (m - n) a^ 



,BI — « — 1 



+ &c. + A"0". 



37. If in Lagrange^s Theorem for A"«^„ n = l, we have 
Aw,= i/' - 1) w„ or A = «''•- 1, the meaning of which is, 
that the operation denoted by A is equivalent to the series 
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of operations denoted by ^^ - 1. And, generally, the series 
of operations denoted by /(A) is equivalent to that denoted 
by f(e*' - 1). For let /(A) be developed in a series of the 
form 

/(A) = J A« + 5A^ + Cby + &€., 
then /(A) w, = u< A«t^, + Bb?u, + Cb^u, + &c. 

or, separating the symbols of operation from those of quantity, 
/(A) w, = {A(/' - 1)« + 5(e^- - 1)^ + C(6*- ^ l)y + &c.} «, 

38. Suppose /(A) = (l + A)% 

then /(c^' - 1) = (I + c^- - I)* «€**'; 
.*. (l + A)"War «* c*^'W* = w*+»> as already proved. 
Again, if /(A)-{log(l + A)}-, 

Sie'- - 1) = {log (1 + e'- . !)}• ^ (log e^-/ = «; 
.-. {log(l + A)}*w, = d;t^„ 

a formula by which the vl^ di£Ferential coefficient of a function 
is expressed by its differences. 

39- To find the general term of the expansion of f(/f) 
in a series ascending by powers of t. 

Writing /(c*) in the form /{i + (^- i)|, ^nd expand- 
ing by Taylor's Theorem, we find 

/(eO =/0) +/(!) (e*- 1) + ri/'O) (c'- 1)'+ ... 
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f 

Then taking, as in Art. 36, the coefficient of f" in each 
term of the second member, and observing that in/(l) it may 

be represented by — i , this quantity bein^^when msO, 7 1/ 

and 2ero in all other cases ; and that in 

1 /<'*>(l) (e«- 1)" it is T^/'kO ^^ 

we have for the coefficient of f* in the expansion of /(«*)> the 
value 

^--T^ 1/(1)0- +/(!) A0-+ ^/'(OA'^O-H- ...} 

- 1^/(1 + A) 0", 

a remarkable theorem first given by Herschel ; for the applica- 
tions of which, see hisi Collection of Examples. 



SECTION IL 



Ioteipatk« of EsfUkit F 

40. The Inrerwe Method of Differaioes lias fiir its olgect 
to detemiine the primitiTe ftmctioii fiom its giren £flSacnoe; 
or from giren relatione between it and its difle naices. We 
ibsll begin with the flmjdest case, 

in which it is required to determine a function whose differ- 
ence is given explicitly in terms of the principal 



41. Since /!iu, is the difference of w, + C, as well as of 
Ugf it will be neceusarj^ in passing from the given difference 
^t^, to the primitive function, to annex an arbitrary constant 
C9 in order to give the result all the generality of which it 
is capable. Also C may be a function of or as well as an 
arbitrary constant, provided its value remains unaltered whilst 

, w changes to ^ 4 1. For if C^ denote such a function of w 
that Cg^x ^ Cgf or ACg • 0, we shall have 

It i« evident that Cj,« 0(2X7r/p) has the property in 
question^ <p denoting any trigonometrical function, sine, co- 
linci bc.y and X any integer. 

42. The symbol 2 is used to denote the operation by 
which we pass from the difference Au^e to the primitive 
function; so that 

2 i^Ug) B f^j. + constant. 
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Also, Us the same function admits of successive diflFer- 
ences, so a function may be integrated any number of times ; 
the second integral of w^, or 2(Sw,), is written 2*w„ and 
the n^ integral 2"«,. 

We now proceed to deduce the integrals of various ex- 
pressions ; chiefly, by reversing the processes given in Section 
I. for finding the differences of functions. 

43. It is evident that 2 (u^ + »* + w„) = 2 w, + 2«, + 2 w, ; 
for if we take the difference of both sides, we get the same 
riesult, viz. w^^ + t>, + w,. And in the same manner it appears 
that 2 (au,) = a^u,. 

44. To find the integral of any rational integral function. 

Since the difference of a rational integral function is a 
function of the same kind one dimension lower, it follows 
that the integral of a function of that description is a similar 
function one dimension higher ; hence, to find the integral of 

we may assume it equal to 

then upon taking the difference of each side, and equating 
the coefficients of like powers of a?, there will arise » + 1 
simple equations to determine the n + 1 quantities a, 6, c,...ft; 
the last term I will remain indeterminate, being in fact the 
arbitrary constant which must be added to mfike the integral 
complete. 

Ex. To find 2(a?*+ 1). 

Assume 2 (^ + 1) = ajfi + 6^* + car^ + da^ + ew ; 

.'. a?* + l = a(5a?*+ lOa^-^lOai^ + 5d?+l) + 6(4.r^ + 6a'* + 4a?+l) 

+ ©(Sa?* + Sci? + 1) + d(2.r + !) + «; 
4 
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.'. l»5a, Omsioa+^b, 0«lOa + 6&-f 3c, 0»5a + 4&+3c + 2c{, 

l^a-^-b + c + d + e^ 

1 11 29 

.". ass-, 6 = — , c = -, desO, c = — , 

or* 0?* .r* 29 

^ 5 2 3 30 

45. To find the integral of the product of consecutive 
terms of an arithmetic progression, we must annex one more 
factor at the beginning, and divide by the number of factors 
:So increased and by the common difference. 

For let u^^ a-h bw, then we have seen, (Art. 11.) that 

'therefore, taking the integrals of both sides, and writing w-^l 
for J?, we get 

(» + 1} O 

Which proves the rule stated above. 

„l(2.-£)(2. + l)(2. + |)(2 



^ + -) + c. 



If one or more factors be deficient in a factorial of this 
kind, it may be resolved into others which are complete, as 
in the following instance; 

(2a? + 1) (2a? + 5) (2a + 7) = (207 + 3 - 2) {2af + 5) (2a? + 7) 

= (2a? + 3) (2a? 4- 5) (2a? + 7) - 2 (2a? + 5) (2a? + 7). 
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46. A rational ifitegral function may often be resolved 
into factorials of the above form, and in this way its integral 
more conveniently found, than by the method of Art. 44. 

Ex. 1. ^ + a?* = a;* (^ + 1) = (^7 - 1 + 1) /p (a? + 1) 

= (a? - 1) a; (a? + 1) + a? (a? + 1) ; 

4 3 

And in general^ any quantity of the form 

aa^ + 6a?""" ^ + ca^""* + &c. 

may be resolved into factorials, by the method of indetermi- 
nate coefficients; thus, if we assume 

oa?^ + 6a? + c = -4 (a? + 1) (a? + 2) + S (a? + 1) + C, 

making a? = - 1, we get a - 6 + c «: C; 

.•. o (a?2 - 1) + 6 (a? + 1) = J (a? + 1) (a? + 2) + jB (a? + l), 

or a (a? - 1) + 6 «= J (a? + 2) + J? ; 

make a? = - 2, .-. - 3a + 6 = 5 ; 

.-. a (a? - 1) + So = -4 (a? + 2), .-. A ^ a. 

In practice, however, it is generally easier to resolve a 
function by inspection, as in Ex. 1, than by this method, 
which is theoretically certain. 

47. To find the integral of a fraction whose denominator 
is the product of consecutive terms of an arithmetic progres*- 
sion, and numerator constant, we must efface the last factor, 
divide by the number of factors remaining and by the common 
difference, and prefix a negative sign. 



9^ >f ^^^^^'ig^ ttioft mt aawe mm. \JkefL liL^ 



«hi«tH pr/yrii* dM: ra!« j-MC iCaent 

«!, If time ^fiymd frtedea, iattad rf k*«iiig ia ■>. 

iiih*' A^.p^M tA ih^ ncimerat^^ ^^^^^ ^ k^>rt k'V'cr br two anhs 
ttfM ihfti f4 i\f^. denominator,) we most reduce the nunientor 
Uf h ii4ff}p% of U^rmn eddi of which is the product of eooaecotiTe 
trn'im^ r^4*,hfniffff tnnn the hepmnng cf the daMmnnator ; that 

Af^ '^ ^liff * -f ,,, + /for + X. - ^' + ^w, + Cu,u,^^ + ... 

th^if i1irt^li;};inf{: th(; nec^md member, and equating coefficients 
of H)ii« )i<nr(^ii <rf'/p, we obtain n^\ equations for determining 
A\ //', f/, ,., ^' ; and the fraction resolves itself into the fol- 
lowing, <>fM5tt {iH which iw integrable^ 



If tlip dp^rpt* of the numerator were only lower by one 
nttit. timtt timt of thr donominator^ we should arrive at a term 

^^ ' I of which we arc able to find the integral, only ap- 
{imHlmatn^Iyi 
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Hence also, if any of the factors of the denominator of the 
fraction in Art. 47 be wanting, they may be supplied by in* 
trodudng them into the numerator and denominator at the 
same time; and then the resulting fraction may be treated 
as in the present Article. 



Ex. 



1 (a? - 1) a? (a? + 1) 



-4 (a?- 2) ... (^H- 2) 
1 3 



(w + 1) (a? + 2) Of (off + 1) (jv + 2) (w- 1) (w + 1) (a? + 2) 

6 

(a? - 2) (a? - 1) 0? (a? + 1) (a? + 2) ' 

which is got by assuming 

(^-1) a? {off+1) « a (a?-2) (^-1) off + b (w-2) (jo-l) + c (a-2) + d, 

and making <v s 2, 1, 0, successively ; taking care to reject 
the factor common to both sides, after each substitution. 

/ 49. To find the integral of a'. 

We have seen (Art. 14.) that Aa'« (a-l)a'; 



a' ^ ., •^ - a' 



... 2a' = +C. Also 2"a' = 



suppressing the part introduced by the constants, which would 
be a rational integral function of the (n - l)*^ degree. 

50. To find the integral of log Vg. 

If u^ = log (t?it>ar8...t?,_i), 
then Au, = log(vir2...t?,) - log(t>ii?2---v,-i) = logt?,; 
.'. 2 log Vg = Ug + log C = log (C t>it?2 ... tJ,--i) = log CPvg^i, 

using Pt?, to denote the product of all the successive values of 
the function «,, from some fixed term tj,, (or more generally 
«», n being independent of <r) to v^ inclusive. 



A*«ii (4r# ^a) - [tAiSi'\ lAn {*tf +a +^(ir+tf)}, (Art. Id), 



inieffriiiing haih fidet n tiiiiet, and refdadng a bj a- ^(«-+0)9 



w« gi;t 



rrfn(/r^ + a) 






H-D 



Again, A rfn 4^0 - 2 sin - co« (/p + ^) d, 

g 

.'. A »in (or - 4) « 2 sin - coso^d ; 

2 



^ ^ Sin (4? - i) ^ 

2 sin - 
2 

COS {070 + a {'JT + O)} 

and X" cos (<rd -fa)- — — 



(*"'"!)" 
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. 52. The preceding expressions may also be integrated by 
substituting for them their exponential values; as in the 
following instance. 

Sa'cosa?0 = ^ S (a'c'^-^ + a'e-*^^ 

* 

^ i" o« - a (e^>^ + 6-^^0 + 1 

. acos(a?- 1)0 - COS/P0 ^ 

= a* .^ ^ ^ + C- 

a* - 2o cos + 1 

Hence, putting a^ ^ 2a cos 6 + 1 » c, and denoting 
afcoswO by t*„ we have c2w, ^a^w^-i - w,, 

and generally, 

1 .2 
Exactly the same formulae hold for Ug^^af' sinwO. 

In the , same manner the integrals of a* (sin /v^)"*, 
a'(cos<v0)'" may be obtained. 



53. To find the integral of 



cos wd cos {w + 1)0 



-,, . ^ sind , . 

Since A tan wd = ^^ -t-^ , (Art. 17.) 

cos wd cos (a? + 1) 

^ 1 tan 070 ^ 

we have 2 2 7 -r-5 = . ^ + C, 

cos /r0 cos (<r + 1) sm 
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54. To find the int^^ of tan'' 



p-^qs-^TJ^ 



Since Atan'H<'+6«)*tan'' ; — — , (Art. 17.) 



we may aMume Stan"* i«= tan~*(a + 6«), 

and take the difference of both sides ; then if the proposed 
function is capable of being integrated, the indeterminate 
coefficients a and b will become known. 



55. The integral of «'w,w,+i... w,+«_i, where Ug=pa^-^q, 
may be determined by assuming it equal to the same ex- 
pression (only with another factor at the beginning instead of 
of) multiplied by an indeterminate coefficient; for 

= (fUgUg^i ... w*+«_i .p (a" — a*"^) ; 



a 



a* 
Similarly, to find 2 , for the assumption 

we must efface the last factor in the denominator, and write 
instead of a' an indeterminate coefficient; for 

1 w*fn"-w. po'(a"-l) 



— — — -^ » 



.% s 



o" 



W,«,+ , ... W,+„ p (O" - 1) W,W,4i ... W,+«-i 
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56, In like manner, if u^,^ a -{- bw^ expressions of the 
forms 

(p H- qw) f (p + ga? + ra?*) if 



can sometimes be integrated, by assuming their integrals equal 
to expressions of the same form, except that the last factor in 
the denominator is effaced, and the polynomial in the nu- 
merator is replaced by another one dimension lower with 
indeterminate coefficients. It is of course only when a certain 
equation of condition between the quantities a^ 6, jp, q^ t is 
satisfied, that this method succeeds. 

Ex.I. Lets, J±l__- „ riHr ; 

(2a' - 1) (2j? + 1) 3' 9.x - 1 



(2 a? 



-1) (2a? +1)3* ^^' V2cr + 1 2a?- 1/ 

= ^(i\*+i _^Li<^±i)_. 

^^^ (2a? -1) (2a? +1)' 



-•. -i = - f , and 2w, = C* - f 



Ex. 2. Let 2 



(2a? - 1) 3' 
a7* H- 6«'i? + 12 -4 + Bx 



a? (af + 1) (a? + 2) 2' a? (a? + 1) 2* * 

X ■\'Z 



then ^ =. - 6, 5 « - 2, and 2w, = C - 



a?(a?+ 1)2'-* 

57. Since A («^xV-p) = w^At?, + «, ^Aw-rj we have 

2(t^^Av,) = w^w^r - 2(i?^4.iAw,), 

the formula for integration by parts, corresponding to the 
formula 

fxUdgV « wf? - fgVdgU. 

58. Change v^ into 2i)„ and consequently A«, into «„ 
then '2 (u,vj) « w,2i?, - 2 (A ^^, 2t?,+i) ; 
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hence, hf pneeemiwe tolMtitotioiis, we get 



2(A*»y2-«,^.) - AX2-^'tT,^.-2(A»^'«,r'+'fr.^..-,i); 
/. 2 («,u,) - «,2i?, - Att,2*«,^,, + A*«,y p,^., 

Ob». In order to be more eaulj remembered, this for- 
mula maj be written 

2(«,uJ « (1 + A2i>)"'w,2i?,, 

where A affects «^, only, and 2f affects 20, ooly, and implies 
an operation such that 

(2i?y2r,-2"^'iJ,+^. 
Hence^ 

V(u,v^) -(!-»■ A2t?)-^^ 2(«,2«,) - (1 + A2t?)-*tt,2?f>,; 

and) generally, 
2" (u^v^) - (1 4 A 2v)"''«,2*u, - «.2-«, - n At*,2»+'r,^., 

n(n 4- 1 ) 

4 \ ; ^ A*«,2»*«t>,^, - Stc, (1) 

which may be also proved inductively, by shewing that the 
eoefflcients of the developments of fi^iu^vj) and (l-^xy^ 
which are identical when n ■ l, undergo the same changes in 
passing from n to n - l . This appears by differencing both 
sides of (1). 

80. The above formula always enables us to find the 
integrals of functions made up of two factors, one of which 
leads to sero, as the value of one of its successive differences, 
and the other admits of successive integrations. Suppose, 
for example, that ti, is a rational integral function of the 
n^^ degree, and that %m a*; then 



Si). 



o* ... o* ^ ' 



- ^, r<)*^.i"7 775, &c,; AX«const, A*+*tt,«0; 
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„. Ufa' Am,o'+' A*«^o'+* A*tt;a'+" 

... :£(„x)-— - —^ . _^_^-...._d_,^c. 

Again, suppose Ug to be a rational integral function of 
the n*^ degree, and u, =»cosa?0; then 

co^\wO (ir + 6)\ 

2X = ^ — 1^-^; (Art. 51.); 



(« «n -) 



2 («*, cos 0)9) = ' ^ 

2 sin - 

2 



. cos|(^p+l)0-(7r + e)} ^, co8{ (j;-H2)e-|(7r^e)} 

(2s,n-) (r 

- &c., 



2 sin - 1 



the series terminating with A*w,. Similarly, if v, « a' cos J7©; 
or if v^ « o' cos"* ofQ sin* ^0, since the product cos* wd sin* w9 
may be replaced by simple dimensions of sines and cosines of 
multiples of ^0; and it will be noticed that the fraction of 
Art. 48 may be brought under this case. 

60* Since the performance of the operation 2 upon any 
series of terms -iA^w, -*- 5A*Wr + •••> reduces it to 

A^sT'^u^ + 5A"**w, + ... ; 

it appears that prefixing 2 to {JA"* + £A* + ...\u^ has the 
same effect as prefixing A~^; in other words, 2 is equivalent 
to A'^ And in like manner, since integrating A'"t«« n times, 
reduces it to A"*"* War? 2* must be equivalent to A"*. 

The same reasoning is applicable to the symbols f^^ d~"; 
whenever therefore, in separating the symbols of operation 
from those of quantity, as in the expression F(^)u^y f{d^Ug^ 
terms containing negative powers of A and d^ occur, they 




ft; T^v /ii*f>*rm;n<^ the jrsneracfn^ fincnra» if 

2'v,. ^w^. ;••., 



6'(A*«,). (i.ij"ci 



>f 






mid 



OL^u.) - (7 - i) Gt^,. 



Oil. Again, ilnce 

/. ... k (m,^^ -w^)r + ... - U-ij Gw, 

A" 

- 1 1 - A log < + (log <)" - &c. - 1} Gu^\ 

thi^l^ftiri), dividing both sides by h and then making A»0, 
whti'h wi» ar0 at lilx^rty to do since h is indeterminate. 
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or 



... + d^u,, f' -f ... = log-. G^w, ; 
.: G (d,u^)^ log -.Gu,, 

and ©(eCtt,) = flog") ^»#- 



Again, log - G (J,u,) = G (d^f,u,) = Gt*. ; 

.-. G! (/,«,) = (log 1] e«„ and (?(/>.) -(log 7) ^«- 

63. To investigate a general series for 2««, involving 
only j^tt^, «,, and the differential coefficients of u,. 

G(2«,)-= iX-\\Gu, = (e'°*^-i)-'G«, 



Gu, 



1 



log- 



log^- 

^ ^ - 1 



- ('<« ;) 'c«.{" - 1 ""g J * 1^ ('-W j)'- »•«■ 



t) 



assuming, as will be proved in the next Art., that -^ can 



e^-1 



be expanded in a series of the same form as that within 
brackets, and denoting by 

-^, -^^^. &c., (-!)"+> ^ 
1.2 1.2.3.4 \in 



l,^ V#<* rfH »%Xw*l* ^ ^ ^ ^^ 









' //^y - V-^^ \^^^ ^ </V, <«, ^ _ 



^ <^:r 






m 



^^-i^^l^B^my* 






A % a M) 












6«) 
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Numbers of Bernouilli. 



64. The numbers J?i, Sj,... which are required in the 

general value of ^u^ in the preceding Art., are called the 

numbers of Bernouilli, and are of great importance in the 
theory of Series. They are defined by the equation 



c*-l * 1.2 1.2.3.4 






(«l)»+i_l!L_L^«igjC., 



and their values may be computed in the following manner. 

The development of -^ can involve no odd powers of 

t 
t above the first. For if (b{t) «= , , we have 



« 

.-. 0(/) - 0(- "= - *» and dt'0(O - **0(- - o» 

or d?<j>{f)^di^{^t). 

Hence, when the development is twice difibrentiated, it will 
not be altered by changing t into (- /), and therefore con- 
tains only even powers of t ; therefore^ {f) can only contain 
one odd power of t viz. the first. Also since 



6*-l / f ' 

1 + - + -f &c. 

2 1.2.3 

we may assume for ^(^) the form of development given 
above, viz. 

«*-l * 1.2 [4 ^ V / |^2n 



m 



2Ji+ 1 

«nd it A^^ be the coeflieieDt of ^ fn this derelopiiient of die 
ieeond member. 



1 (2 ny - 2 n (2 II -!)*■+..: . 

-»- . j— T«C- 

2« + I [2n^ 

[2n ' * * 2n + 1 ' 

all the termn after A^'O'" yanishing, since A*0^" is zero when 

m > 2w. 

Hence, 
/?...! -(-irM-i^0*-+^A'0«--&c.... + ^^A«-O^«}. 

By this formula if|, J?,, i^s, &c. may be readily com- 
putedy supposing the numbers comprised in the form A^O'^" 
to be known; we find 

11115 

^»"6' ^-"i5' ^*"=ii' ^'^^Vo' ^'^66'^'' 
68. Also, since 
A"*©** - »•■• - m (m - 1)»" -I- ^ ^ ^ (w-g)*" - &c. (Art. 20.) 

1.2 

we may, if we please, eliminate the numbers A'*0^ from the 
expression for Bun~\\ and we find 

|2n 
(.ir»Bt..^--i4i(««"-«)-i(S^-S.2*+5) + &C.+ I— 



2fl+l 



41 

Besides serving to express the general value of Stc^,, the 
numbers of Bemouilli have various other uses, of which we 
shall now give one or two of the most remarkable. 

66. To find the general terms of the expansions of 
cot and tan in powers of 0. 

Now the general term of the expansion of 

or i 

[2n 

.*. the general term of the expansion of cot is 

and cot = B^0 ?— ^ - &c. 

1.2 1.2 .3.4 

Also since tan d s cot d ~ 2 cot 2d, 

the general term of the series for tan is 

4- 2 . 



2w * Wn 



2^(2^-l)g,,., ^„, 

[2n 

^ 4.3^^ l6.15 ^ ^, 2«(2*-l)„^ „ 

.-. tan = S,e + 5305 + —5^ / Sftfl* + &c. 

1.2 1.2.3.4 1.2... 6 
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Hence, by differentiating the above expressions for cot 6 
and tand, we may deduce the general terms of the expan- 
sions of cot* and tan* ; and by integrating them, the 
general terms of the expansions of log sin d, log cos 0. 

67. To find the sum of the series 

j«n q2m gin ^ -" 

or the value of ^^J* -— . 

S.e...«..{,-(?)]{,-(|)](,-(A)]... 
changing into ird, we get 

-—{-(?)*}{'-©]{-©]••• 

therefore, differentiating the logarithm of each member, 

«i -1 



-i{'-(i)r'-- 

But the coeflBcient of 6*"~* in the expansion of tp cot nrd is 

[2n ' 
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JB^" 1 .2. 3 ...2w * 

„ 111 TT* 1 1 1 TT* 

Hence .--4.--+ — +... =--, — + -. + -- + ...=—. 
1* 2* 3* 6 1* 2* 3* 90 

Obs. Calling S2n the sum of this series, we have 

S^, ^ (2n + 1) (2n + 2) gg,,! 
02114. 2 ^^ "211+1 

Now suppose » very great, then — = — , which 

proves the divergency of the series formed by the numbers 
of Bernouilli ; these numbers increase very rapidly, beginning 
with -B13. 

68. To find the sum of the series 



or the value of ^S 



00 



(2j? + 1) 



9n 



Since cos 



-{-©"K-^H-d)}- 



changing 6 into — , we get 



COS 



"i'Hm-m-it))- 



therefore, differentiating the logarithm of each member, 

_/l r%/\ e //lv21-l an f /i3\ 21-1 

TT iru 
— tan — 
2 2 



?i-(f)r"^F{'-m 



f{-©} 



81-1 

+ &c. ; 
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and equating the coefficients of 0**'* in each member, which 
are respectively, 

we get •^•- ^^Ig,- /''^^"^^ - 

'^ (2jr+ !)«• ^ 'l.2.3...2n 

111 11 Stt^ ir^ 

Hence -r + -; + -; + &c. = - . - . « -- . 

1* 32 52 2 6 1.2 8 



i.^ At 1 ^ 1 ^ 

69. Also T-+ -T r + --- 

j2« 22. 32. ^2» ^ 

-^+ ii;i + ^ + — "i5;(ji;^+^'»- 32.+ — j 



12*»-1 . _ 1 g**-*-** 






2n 



2 [2» "" " 2^ 

2» 



70. Since log (l + e"') = c"* - ^ e"'* + J e-"* - &c., 

integrating this 2»-l times between the limits « = 0, ^ = 00, 
we find for result the series just summed ; 

(2^-^-l)^"5«,., 



... cry-4ogo+o = 



\2n 



71. To find an approximate value of 

r(^ + 1) « 1 . 2 . S..,a?5 when w is very large. 

Making w, «loga? in the formula 

1 • z 
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we find, (Art. 50.) adding logo? to both sides, 

log {1.2. 3,.. (j?- l)a?} =C + a? log a? -07 -^log^ 

B, I B^ \ 

+ -—+... + log J? 

1 .2 w 3.4 a^ ^ 

= C + (a? + ^) loga? - a? + log (1 + A), 
putting log (1 + A) = :^ - ^-^ + ^ - &c.. 

SO that A is a quantity continually approaching to zero as a in- 
creases. 

Now to determine C, suppose x very large so that h may 
he neglected, and change x into 2^7, then 

log (1 . 2 . 8 ... 2a?) = C + (2 J? + ^) log 2ci? - %w 

= C+(2a? + j^) (loga?+log2)-2a?, 
and log (2 . 4 . 6 ... 2ar) = log (2*. 1 . 2 . 3 ... ar) 

= a? log 2 + C + (/F + "l^) log X - w'y 

.-. log {1.3. 5. ..(2a?- 1)} =a7loga?+ (a7 + |^)log2 -a?; 

= C + ^ log2a7 - log2; 

^ , , 2.4. 6. ..(2a?) - 

.'. 2C - 2 log 2 « 2 log p— ^ — ^ - log 2a? 

1 . 3 . 5... (2a? — 1) 

« . 2 . 4 . 4 . 6 . 6 . . . (2 a? - 2 ) 2 af 
^^ 1 .3.3.5.5.7... (2af-l)(2a?- 1) 

= log — , by Wallis's Theorem, 
since a? is indefinitely large; 
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.-. C = ^ log 27r; 
.•. log(1.2.S...a?) ~^log 27r + (^ + ^) logo? - a? + log(l-i- A) 



= log \/^irw + log ( - J 4- log (1 + A) ; 



.-. 1 . 2 . 3 ...a; s y/ 2irw • ( - ) •(!+*), 

where A is to be calculated from the series 

fij 1 ^3 1 B^ I 

log(l+A) = — ..- — .- + — .-.&c., 

and in general it will suffice to take the first term only, 
which gives 

1 . 2 . 3 ... a? = v^27ra7 afe '^^ 

Obs. The preceding series, even for large values of a?, 
becomes divergent after a certain number of terms ; this will 
happen after n terms if 

-"211 + 1 1 -0211-1 ^ 

(2» + 1) (2» + 2) ' a?^+» (2« - 1) 2W a?*«-i ' 

g2n + i ^ (2n -t-l) (2n + 2) ^ 
or — > ar^ 

^2.-1 (2n-l)2» 

but the first member of this inequality never exceeds 

.'. (2w - 1) 2w > 471^0?*, or w > irw. 

It can, however, be proved that an approximate value 
of the series will be obtained by taking the aggregate of 
the convergent terms only. 
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SECTION III. 



EQUATIONS OF UIFFEUENCES. 



72. We now come to the case in which the relation 
between the principal variable and any function of it is to 
be determined by means of an equation between ^, u^^ and 
one or more of the successive values, or differences, of u^\ 
that is, from equations of the form 

F (wy «„ «,+!> ... w,+«) = 0, 

or, /(w, w„ AWx, ... A" w,) = 0, 

since, by the theorems of Arts. 18 and 23, these forms are 
convertible one into the other. An Equation of Differences 
is said to be of the »*^ order when the successive value, or 
the difference, of the highest order which it involves is the 

73. The complete integral of an equation of differences 
of the n^^ order will contain n arbitrary constants. 

Let ... «,, w,+fc) ^x + 2h^ ••• ^^ * series of terms corre- 
sponding to the successive values a?, ^h-^j ^ + 2A, ...; and 
let 

be the equation by which the general term is determined as 
a function of a? and a, or the equation of the series, a being 
an arbitrary constant. Since this equation must hold for all 
the succeeding terms, we shall have 
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Eliminating (a) between these two equations, we get an 
equation between a?, u^j and u^^f^i or, substituting t^, + Af^^^ 
for u^^^9 an equation between a?, w,, and Aw.,, which is the 
equation of differences of the first order whose primitive 
equation is 

F {w^ Ugf a) = 0. 

In like manner if the equation of the general term con- 
tained two arbitrary constants a and 6, as 

F (<r, w„ a, 6) = 0, 

we might eliminate a and b by means of the two succeeding 
equations, 

F(a;+h,u^^^, o, &) = 0, J' (a? + 2 A, w,+2;k, a, 6)=0, 

and thus get an equation between a?, «,, w,+a> ^x+2*^ ^**> 
substituting 

w, + A«, for t^,^^, and w, + 2 At«, + A^w, for w,+2fc» 

an equation between ^, Ug, At^«, A^u^, without the constants 
a and 6, which is an equation of differences of the second 
order, having for its complete primitive the equation 

F(Xy Ua.9 «» b) =0. 

Hence it appears that every equation of differences of the 
first order, or between two successive terms, will introduce 
one arbitrary constant into the equation of the series; every 
equation of differences of the second order, or between three 
successive terms, will introduce two arbitrary constants into 
the equation of the series ; and, generally, every equation of 
differences df the n^^ order will introduce n arbitrary con- 
stants into the equation of the series. 

Linear Equation of Differences of the First Order. 

74}. The general equation of the first order and de- 
gree is 
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Ag and Bg being functions of w. To integrate it, assume 

and in order that this equation may resolve itself into two 
others, each of which admits of being integrated, assume 
(as we are at liberty to do, having made only one supposition 
respecting v, and wj) 

Vg^.1 Wg — Ag o, fOg » 0, 
or, dividing by «,w„ 



But Alogt^^alog — i, .•. Alogv^slogJ,, 

Vg 

. \ log V, = S log Ag = log P J,_ , , (Art. 50.) 

or Vg= PAg^if the constant being unnecessary. 
The other part of the equation gives 

and «. - ^^-1 {2 (^) + cj , 

the complete integral, involving one arbitrary constant. 

7 
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Obi. Taking the difference of the result 

PA,., ^ \paJ + ^' ^^ «** PA. PA. ' 

which shews that -— r- is a factor which makes each side of 

PAg 

the proposed equation integrable; and it is generally the 
most convenient waj of integrating the equation to multiply 
it by this factor. 

Ex. 1. «,+! — aUg = ^. 
Here Ag = a, PA^ = o' ; 

.-. ^n = 2 \^^ - 2 (pi^a% puttmg - = a, 

« ^^ -*• 5 4-C; (Art. 59.) 



a - 1 (a - 1)* (a - l) 



• • 



Ex. 2. w^+i - aug = cos ZP0. 

cos (j^ — 1)0- a cos /rd ^ x 

«*' ^ — T— ^-^^ — 5—; + ^ «*• 

o*-2acos0+l 

Ex. 3. Two vessels which hold (a) and (6) gallons re- 
spectively are filled, the one with proof spirit, the other 
with water ; (c) gaUons are taken from each and poured into 
the other; and this is repeated such a number of times as 
to make their contents of the same strength ; find the number 
of times. 



w 
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Indirect Integrals of Equations of Difierenoes. 

^5. Since the equation of differences of the first order> 

/(a?, u^y Aw,) = 0, 
is formed by eliminating the constant (a) between the equations 

u,^F(afj a), t^,^^ - F(af + A, a), 

it follows that we shall arrive at the same equation of differ^ 
ences, whether a be constant, or be a function of ^, as a^y. 
provided it be such that 

i^(a? + A, a^+j)--F(w^hy a,). (1.) 

Now this equation is satisfied by a;^^^ = a„ which givea 
Aa, aO, and a, = a, a constant, and leads to the ordinary 
or direct equation to the series, 

Ug = F(Wy a). 

Also equation (l) will be divisible by o,^^ — a„ be- 
cause fij^ is a value of o^+^^y which satisfies equation (l); 
and if aimensions of a^^^ and a, superior to the first are 
involved in it, the result of this division will be an equation 
involving a^^j^ and o,; i. e., an equation of differences of 
the first order with respect to a«, the solution of which will 
give one or more values of o, in terms of w and arbitrary 
constants ; and these being substituted for a, in the equation 

will furnish equations of series, which are primitive equa- 
tions of 

/(^, «,» Aw,) = 0, 

and each involves an arbitrary constant. 
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If equation (l) does not involve higher dimensions of a, 
and a^^j^ than the first, they will disappear from the result 
when it is divided hy a^^^^ - o,. In this case a^ will have 
only one value, viz, a^ = a, and there will be only one equa- 
tion of a series corresponding to the proposed equation of 
differences. The mode in which the indirect solutions just 
treated of are obtained, is analogous to that in which the 
singular solutions of differential equations are obtained; but 
whereas the latter can contain no arbitrary constant, indirect 
solutions of equations of differences may contain as many 
arbitrary constants as the complete integral itself from which 
they are deduced. 

Ex. w, -»^Aw, + F(A«,). 
Taking the difference, we find 

A«, = Au^ + (a? + 1) A*w^ + AF(Aw,), 

or = (a? + 1) A^u^ + AF ( Att,), 
which is evidently satisfied by Au^ — a, a constant; 

Ug^ aw + a' ; 
and substituting in the proposed equation 

aw ^ a' ^ aw + F (a), .•. a' ^ F (a) ; 
.'. u,^wa + F(a)y 
the complete integral, containing one arbitrary constant. 

For the indirect solutions we shall have 

u^ « wa, + F (aj, 

a^ being determined from the equation 

i'V+ l)a, + i-f i^(a,+i)« (w + l)a^^F(a^). 
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Suppose, for instance, that F (a,) « a*, 

A (x + 1) (a,4.x - o,) + «;+, - «; « 0, 

or, rejecting the factor a^^^ - a^, 

a 2a? + 1 . 



or «,« +C(-1)'; 






Linear Equations of Differences of all Orders. 

76. The linear equation of Differences of the n^ order 
is 

all the coefficients being functions of a ; the first step towards 
its integration is to establish the following theorem. 

77* If there be n particular values 
which, when substituted for u^^ satisfy the equation 

that has no term independent of »„ its complete integral ifi 

^19 ^y*<^ii being arbitrary constants. 



S4 

For let thfs^ valued be substituted in the expression 

and it becomes, (collecting the terms multiplied by the factors 
Oi, 02 9... a,) 

«i(X+«+l>i^«^+— 1 + ... + p«X) + «»(^'+» + Pi *«**+«-! + ••• 

+ P» *«*,) + — + «« C«^*+« + PiX+— 1 + .. . + p«X)- 

Now since *«,, *e*,, ...*«*,, satisfy the proposed equation^ 
each of the quantities included within brackets is equal to 
zero, therefore the whole is identically zero ; consequently the 
assumed value of u, satisfies the proposed equation, and it 
contains n arbitrary constants, therefore it is the complete 
integral of that equation. 

78. To integrate the equation of differences, 

^M+n + Pl^s+n-l + P2^M+m^2 + f + Pji«* « 9> 

all the coefficients and q being constants. 

Assume u^^v^-^k; then by substitution we get 

^s+n +Pl^*+n-l +P2«>*+n-2 + ••• 

+ p„«, + k(l +pi +P2 + — +Pn) -q^O. 

q 
Let k = , then the equation becomes 

I+P1+P2+.- +Pii 

«»+« +Plt'4-+«-l +P2»#+«-8 + -. +Pn^s ^ 0. (1.) 

Let Vg = a% then a'{a* + pie»*"* +P2«""* + '•• +Pn) is the 
value of the first member ; now this will vanish if a be any 
root of the equation (called the auxiliary equation), 

/(a) = a" + pitt*"^ + P8«"~* + ••. + P«-ia +Pn^O. 
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Hence the (n) roots of this equation a^^ a^y o^ymm.Un 
will give n particular values of t;,, a'» af> a^.^^al which 
satisfy equation (l) ; therefore its complete integral is 

and the complete integral of the proposed equation is 

q 



1+Pi+Pa+ ...+p, 



79* If the auxiliary equation have equal roots, the 
above ceases to be the form of the complete solution ; because, 
in that case, it does not involve the due number of arbitrary 
constants ; and it must be modified as follows. Suppose two 
roots a^f ^2, to be very nearly equal to one another, so that 
€ij^=ai + hf h being a very small known quantity; 

then, 
C\a\ + c^al = (cj + Cg) a' + c^ {ara' "' A + a' "*A* + &c. } 

= Cxa\ + C2 \x€l[ + — ^ ^af-^A + &c. J 

replacing the constants 

Ci + Cg by Ci, and — ? by C% ; 

now this continues true however small h be taken, and 
therefore when A s 0, when it becomes 

(Ci + C^x) a\ ; 
.\ Ug - {C^ + dai) at + c^al + &c. 

Similarly, if the auxiliary equation have r roots equal 
to O], the complete solution will be 
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of the correctness of which, we may be assured by Ae fol- 
lowing reverse process. Assume Ug^a'Vgy then 

and the first side of the equation becomes, when divided 
by a', 

a» (I + Ay«, + p, a— > (1 + A)*-*«, + &c + p,t>, 

- {{a + oA)" +Pi (a + oA)"-' + &c. +p,\ v, 

-/(a + o A) ». = {/(o) +/(o) . o A +/' (o) ^ + bc.]v, 

1 • Z 



o ^ . . - a* 



«/(o) •ffg-k-- *f (a) . Ar, + /' (a) . A*«, + &c. + a»A*»^. 

1 1*2 

Now suppose a = Oi, and /(a) » to have r roots equal 
to Gi ; this makes the terms as far as /^'"^^(a) vanish ; and if 

t?^asCb + Ci/r+ ... +c,.i^^~*, 

then A^t?^ = 0, A'+*r, = 0, &c., 

and all the remaining terms vanish; and consequently the 
equation is satisfied by 

Ug s (co + CiOf + ... + c,_iaf'"0 »'• 

80. Also if the auxiliary equation have a pair of ima- 
ginary roots, 

mAn\/^ = /o(cos0 ±\/^sin0), 
putting p « V m* + «*, tan = — , 

the corresponding terms in the value of u, will be 

Cp'(cosd + v^ 8in^)' + Cy (cos0 - v^^sin^)' 

x3 p' (C| cos ord + Ct sin ^d), 
changing the arbitrary constants. 
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And if there be r pairs of imaginary roots, the corres- 
ponding terms in the value of «, will be 

(ao+ «ia? + a2^+... +ar-i^''"Op' (cosd + v-l sin 9)' 
+ (60 + *i^ + ••• + *r-.i^'"0 P' (cos 9 - yZ-lsinOy ; 
or, changing the arbitrary constants, 

(co + Cia?+ ... + Cy_i<r'"^)p'cosa?0 
H- (c'o + c\a? + ... + Cr_iaof''~^) p' sin 476. 

Ob 3. The Differential Calculus being a particular case 
of that of Finite Differences, a strict analogy exists between 
the methods and results in the two subjects, as the reader 
cannot fail to have observed; indeed whenever in the latter 
a result is obtained with an indeterminate increment for the 
principal variable, it is possible to pass to the corresponding 
result in the former, by a method similar to that pursued 
in Art. 25 ; and which we shall further illustrate by the fol- 
lowing instance. 

In the equation w,+2fc "" ^^^x+h + (^^ + ^0 ^* == ^> 

i_ 

putting w, = a', we find a = (m ± w\/— l)*; 

Now the proposed equation is 

A^u^ - 2 (wi - 1) Au, + {(m - 1)* + n^} u^ = 0, 

or if we replace the known quantities w — 1 and n by other 
known quantities m^h and n^A, the equation and its solution 
take the forms 






8 
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u, « {(1 + m,hy + w?A«}^ |c^ cos {^ tan-' fT^) 



+ Co sin - tan"^ r > 

' \h 1 ■¥m^hj f 



Now take the limits of these expressions when A = 0, and 
we find the well known results, since (l+2wiA)** becomes e""*, 

djw, - 2wiid,w, + (mj + n\) w, = 0, 
Ug » e"*!* (Cj cos WjO? -f- Cg sin w^ a?). 

81. To integrate the linear equation of difierences of 
the n^ order, 

the coefficients being functions of w. 

If 'w,, *w,, ... "^u, be 71 particular values of v^ in the 
equation 

with which the proposed coincides when its second member 
is zero, we have (Art. 77.) 

If we now divide both sides by 'w,, and take the difierence, 
we shall eliminate c^; next dividing both sides by the co- 
efficient of Ca, which suppose 't?,, and taking the difference, 
we shall eliminate Cg; again dividing by *«,, the coefficient 
of C3, and taking the difference, we shalL eliminate c^; and 
proceeding in this mannec till all the constants are eliminated, 
pur final result wiU be of the form (each A affecting the 
whole of the expression that follows it) 

11 1 «, 
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in which expression the coefficient of i>,+» is evidently 



^»+l' ^x+2'" ''*+»-l ^*+ii 



therefore, dividing by this coefficient, we get the expression 

1 1 I , V, 



tt-i 



»*+i* \+2'*'^^»+n-i^^s+n^-;^zr' ^^i^^^r~'"^~ ^r^» 



«?* '^J, v^ «* 



which must be equivalent to the first member of equation (l). 
Therefore the same expression, only with u, instead of v«, 
must be equivalent to the first member of the proposed 
equation, and consequently equal to JT; hence, equating these 
equals, and integrating, we get 

w, = X2^»,2V .-. 2""V, S^jn ;n 1 1- — » 

(each 2 afiecting the whole of the expression which follows it) 
which is a general formula for the integration of any linear 
equation of differences whatever. 

82. As a first exemplification of this method, suppose 
the coefficients of the equation to be constant, and ai, a^y ••• a^ 
to be the n roots of its auxiliary equation. 

Then «, = c^ a' H- Cg al + Cial+ ... + c^ a* , 

changing, both in this, and in the similar succeeding steps, 
the arbitrary constants; 

ApyAp)-Ai;.e.p)%...+<,(5.y 



\ «„ / \««-l/ \^2/ «! 
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in which expression the coefficient of «,+» is evidently 

1 1 1 



< ««a«-i-..a2ai (-1)"P»« 



n 



therefore, dividing by this coefficient, and replacing v, by u,y 
we get 



( 



or, if we choose to introduce the arbitrary constant after 
each integration, 

«, = Cj a* H- Cg ttg + ... 



83. If the auxiliary equation contain equal, or imaginary 
roots, this method is still applicable; it is only necessary ta 
assume for v^ the value belonging to the case of equal or 
imaginary roots, as will be seen in the following instances. 

Ex. 1. To solve the equation 

Here u^ = (c + c^w) — ; 



«• 



.*. A^(a'«,) = 0, 



* For this method, and the corresponding one of solving the general linear 
differential equation, I am indebted to Mr Gaskin. 
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in which expression, the coefficient of »,+, is a'+*; therefore, 
dividing by this quantity, and replacing v, by u^, we get 



a*+* 



a 



Let Jr = a?, then 2^(a?a') = a?Fa' - 22V+^ +c'^ + c, 

(Art. 58.) 

1 r war 2a'+^ , 1 

a'-^Xia-iy (a -If j 



1 a* J? 2a^ 

= (C + C^47)- + 7 ;^2 - 3 ,.j 

a (a — 1)^ (a - 1/ 



Ex. 2. To solve the equation 
Here t?, = Cip* sin a?0 + Cgp' cos w9. 



where p^ ^m^ -{- n\ tan = — ; 



«* 



= Ci tan 070 4- Cs ; 



'. A 



p'^ cos w6 

Ci sind 



\p* cos w9j 



p* cos a?0/ cos wQ cos (a; + 1) ' 



.•. A cosci?0 cos(a?+l)0A (— — ^— -r|=0, 

^ Kp'coswOj 

in which expression, the coefficient of t?,^.8 is g » 



» 
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therefore^ dividing by this coefficient, and replacing v^ by u^, 
we get 



COS 



^ —-. A cos wd cos (a? + 1) A ( — — "^—pA = X; 

(^ + l)fl V/o cos/r0/ 

/»^ /. ^ V /.^ f-X^cosfayH- 1)01 
.-. w, =// cos 0702 sec d?0 sec (jf + 1)02 < ^^j^ — ^>; 

to which we may add the terms Cip'sina7 + Cg/o'cos^S, if 
we suppose a constant to be added after each integration. 

Ex. 3. To solve the equation 

^*+4 + Pi Ux+3 + P» ^x+2 + 1^ W,+ i + P4 W, « JT, 

where the auxiliary equation has two pairs of imaginary roots, 
so that 

Here t?, = (a + o' <r) ^o'cos 07 + (6 + 6' w) p^ sin a? 0, 



s: a + a* Of + (b + b'of) tan a? 0, 



.-. A 



p* cos 07 

— -: « a + 6 tan 070 + (6 + 6 0?) -i r-:, 

p* COS *P cos 0? COS (07 4- 1) 



.'. COS Of COS (O? + 1) A ;: = ia'jCOs(2 07+ 1) + COS 91 

p'coso?0 » ^ ^ ^ ^ 

+ ^6'{sin(2o? + l)0-sin0} + (6 + 6'o7)sin0, 
.-. A cos 07 cos (o7 + 1) A -r ;: = - o' sitt 2 (o? + l) sin 

p COSW0 ' 

+ h' COS 2 (o7 + 1) sin + 6' sin 0, 



.*. sec'^ (o7 + 1) A cos 07 cos (o7 + 1) A 



^' cos 07 

= - 2 o' tan (o? + 1) sin + 2 6' sin 0, 



J 
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2a'sin*e 
^ ' ^ ' cos(^+ l)0cos(a? + 2)0 

.'. A cos ((JO + 1) ^ cos (a? + 2) A sec* (af + 1) x 

A COS J? cos (a? + 1) A ;: = 0, 

p* COS X 

in which expression, the coefficient of t?,^.^ is ^ ; 

therefore, dividing by this coefficient, and replacing v^ by w„ 
we get an expression equivalent to the first member of the 
proposed equation ; hence, equating it to JT^ and integrating, 
we find 

w, = |o' cos <r S sec J? sec (a? + 1) 2 cos* (w + 1)0 x 

2 sec (ci? + 1) 0sec (^ + 2) 2 1 \:^^ — —\. 

84. Having given a particular integral of the equation 

to find its complete solution. 

Let Vg be a particular integral, then 

Hence, eliminating Ag, we find, putting — ?»jir,, 

A»,+i B^v^ 

T-- « = tii„ suppose ; 

A x„ t?,+2 

or A log (A xj) = log ««7, , 

.-. A «, = CPw^.^^ (Art. 50.) 
.'. «, = C«, 2(Pw,_i). 
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85. If we know a particular integral of a linear equation 
of any order that has no term independent of w„ we may 
reduce it to another equation of the same kind of the order 
immediately inferior. 

Let Ug = t7, be a particular integral of a linear equation 
of differences of the n^^ order reduced to the form 

A"w, + 5fiA""^w, + 5^2^""*^, + ... H-gn«**=/(A)w,«0 (i). 

Assume w, = t?,Sw, ; 

then A"w,= {A + (1 + A)A'}»t?,2w,, 

where A affects «, only, and A' affects 2m?, only, (Art. 34.), 
and the proposed equation becomes 

/{A + (1 + A) A'} «J,2ti7, =0. 

Now /{a + (l + A) A'} is a rational integral function 
of A and A\ which we wish to arrange according to powers 
of a'; and this may be done at once by Taylor^s theorem, 
which gives • 

/(A) + A7.(A)(l+A) + ^V.(A)(l + Ay+...+A'"(l + Ar. 

Hence, observing that /(A) «x = 0, since «, substituted 
for u^ satisfies the equation, and that (1 + A)'*?, = »*+r> we 
get the depressed equation 

w,/i(A)tJ,+i+-- Aw,/g(A)t?,+2+ ... + A"-'<.i>,+,= 0; 

or, reversing the order of the terms, (since f^ (A) means the 
same function of A, that dlfijxi) does, of a?), 

+ /i(A)«,^.i.w, = o, (2.) 



1 
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a linear equation of the (n — l)^ order, of the same form 
as the original one. Similarly, if we know another particular 

value of u^y «r^, Ihen A (^ j wiU be a value of w^ in equa- 
tion (2), which may be depressed to another of the same 
form of the (n — 2)**^ order ; and if we know r particular 
solutions of equation (l), we may in this way depress it to 
an equation of the same form of the (n - r)*^ order. As 
a linear equation of differences of the first order and degree 
can always be solved, it appears that to obtain the complete 
integral of a linear equation of the n*^ order, we must know 
» * 1 particular solutions. 

86. Besides linear equations of differences, and such 
equations as can be deduced to that form, very little is 
known of equations of differences of the second and higher 
orders. The following are instances of equations which admit 
of reduction to linear Equations with constant coefficients. 

where pi^ j^g, &c. are constants, and v^ is any function of 
w. Assume u^^iWtgPv^^ then the equation becomes divisible 
by Pvg+n9 ^^^ is reduced to the linear equation with constant 
coefficients, 

^»+n -t- Pl^x^n-1 + P2^s+n-2 + ••. + Ph^» = 0- 



2, w,+. + pia'u^^^.y + Piioru„^n-% + ... + p«a"'w, = 0. 

Assume w, «t?,a^***"^' '*^'^, then any term p^ a'* «,+,«, 
becomes 

therefore the equation becomes divisible by a*'('+^»+*), and is 
reduced to the linear equation with constant coefficients, 

9 
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3. Uj,^t + {a + ft (- 1)'} «,+ J + cu, - 0. 
Let u^ as r, y/a-k- 6(-l)% 
then «,+, = t?,,.! V a - ft (- 1)% 
«#+f « «^*+8 \/« + ft (- 1)' ; 

Assume «, + a = -^^ , then «,+i -♦- a « -^— , 






+ c = 0, 



or t>,+2 - (o — 6) i?,+i + (c - aft) t?, = 0. 

The two arbitrary constants which will appear in the 
value of Vg9 must be reduced to a single constant, bj the 
condition of the proposed equation being satisfied. 

5. t*,4.i - 2wJ+ 1 =0. 

Assume «, « cos «,, then w, ■= ^ (c** + c"**). 



Prob. 1. To find the value of 

V 2 - Vg-VVr- .,. to zr terms, 
let it =a M,, then u^^^ «= vi — «J; 
.•. ttj+i + te, ■= 2 ; let w, = 2 cos t?„ 

then 4 cos* t?,4., + 2 cost>, =s 2 ; 
.'. 2 cos* tJ,+, = 1 - cos r, a 2 sin* ^ v, ; 

cos Vg^i » sin -^ V, « cos f ^ v*) » 
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W 



•*• '^jr+l + "g"^* ^-p » 



•• »-="c(-i)'+|; 



U.S2.COS 



{o(-i)-^f}; 






+ - «= V 2 = 2 cos — , .'. 

.-. t*.«.2cos|^(-^)' + j|^28in|{l-:(-i)-}. 

Prob. 2. A swan breeds three cygnets in its second year, 
and four every succeeding year; and the young ones all 
breed according to the same law; required the number at 
the end of the w^^ year. 

Let u^ = number at end of J7*^ year ; 

then in (^ + l)^ year, ^here will be bred three each by those 
in their second year, and four each by all the others ; 

••• t**+8 - «**+2 - Stt,+i - w, = ; 
.-. «»'-w*-Sm- 1 « f»(i»'-l) - (w+ iy«0; 

.-. f»=-i, and l±\/2; 
.-. «, = a(-l)' + 6(l+\/2)* + c(l-\/2)', 
t^Q s s a + & + c*, 

t*j = l = -a + 6(l+ -y/i) + c (l - \/2)j 
t^s s 4 := a + 6 (3 + 2 \/2) + c (3 - 2 y/i) ; 



es 



a»l, 6 



•■• «.-(-!)• ti(^-.)(i+v^)--i(^+i)(i->Ar. 

If the law of increase be only one the second, and every 
succeeding year, 

^-.-(■-^T-(^r 



87- The following problems relative to continued frac- 
tions, give rise to equations of differences, which can be in- 
tegrated as linear equations. 

I. To determine the value of the continued fraction 

u, = to w n-actional terms. 

o -»- a + ... a 

c 
.-. u,+^ « — — , or u,^i (u, + a)^c; 
a + Ug 



assume 



tt, + a = , then ( a 1 = c, 



or t>,+2-ai>,+i-ci>;, = 0; 
,*, t), ■= Ci a' + Cj /S*, a and ft being roots of A*-oAp-c = 0; 

*^- — ;r^^ /ai > putting C8«=-. 
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a a +€iP a 



a •-/3'+^' 



.*+i 



If a and /3 be imaginary, so that 

, sin^0 
then Ug^ - p -7—7 r-;: . 

^. To find the value of the continued fraction 

c c c c 

to w fractional terms» the last term being -- . 

b 

Proceeding as in the last example, we find 

a'-' + Cs/J*-' 



«, = c. 



a' + Csfi' ' 



, c l+Cs 6-a 

and «i»7=c-— — -Qj ••^s'^-T — 5> 
& a + e^p 6 — p 

3. To find the value of the continued fraction to a; 
fractional terms, 

1111 



«# = 



a+ 6+ a^- 6+... 



1 1 1 

Here tt,+2 = , - ^TTTT^ 

0+ r 

6 + «, 

or 8 6 4- »«• 

1 - a«,+. 
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Assume 1 - a w, « -^ , .-. 1 - a «,+» « -^* , 

... lf,-:^A^.6^ifi_!=nY 

or tJ,^3 - (2 + a 6) t?,+, + r,_i = 0, 
of which the auxiliary equation is 

Let a, , - a, - , be the roots of this equation, then 

a a 

tJ, = a' {c, + a, (- 1)'} + ^ 1^3 + C4 (- 1)'}, 
where o «= ^3 + <^4 (- l)^"' ^ ^ , .,.i . 



a 6 a + j92 a 



-1 



aWjsa -s= 1 - 



.8_i_ *» ^-3* 



.3 . ^ _-3 



a +©2 a 1 / — 

.-. a6 + l = a*-l+ o"* = r-y '.' a — = \/ab ; 

a + p2 «" o 

.-. Pa = 1 = C - C, .'. C = 0, C' = - 1, and p,» (- 1)', 

a'-'4-(-l)'a-'+' 



/. a w^ = 1 - 



a'+' +<-!)• a—'* 
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4. Suppose 6 to be essentially negative and « - c, so that 
the continued fraction is 



«+• -c+ a + - c + ... 



and let \/ah « \/^\ y/ac = 2\/- 1 sin 0, 

then ft* - 2\/^ sin 0& - 1 = 0, 



. .% a = cos0+\/-l sin0; 



sin (a? - 1)0 

.% at^« a 1 — ; — ^ — , 

sin (/r + 1) 

according as iv is odd, or even. 



or 1 - 



cos {OD — 1)9 
cos(^r +1)0' 



SECTION IV. 

SUMMATION OF SERIES. 



Integration of the General Term. 

88. One of the most direct and important applications 
of the Calculus of Finite Differences, is the general method 
which it furnishes of assigning the sum of any number of 
Terms of a Series of the general term of which we are 
able to take the integral. There will be two cases to con- 
sider, according as the general term is given explicitly in 
terms of the index, or is only given by means of an equation 
of differences. We shall begin with the former, by shewing 
that the sum of any number of terms ending with the 
general term u^^ is equal to the integral of the following 
term, plus a constant. 

Let S, denote the sum of the w first terms of a series 
whose general term is u,^ 

then iS, = «i + W2 + ••• + *^*» 

*5'*+i = «i + «^ + ... + «^ + «,+i ; 

.*. iS, = 2m,+i + C. 

Making a? = 0, So = 0^ 2«,^i + C; 
or, lis h is uisually written, 

.'. Sg^ 2w,+i - 2wi. 

In general the arbitrary constant will be determined by 
the term with which we make the series commence. Thus, 
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if w^ make it begin with Ur instead of Uu this amounts to 
suppose the sum of the terms preceding u^, that is, 2t^^ to 
be zero; and therefore to determine the constant, we have 

Sr-i - ^U,^r + C ^ 0. 

89. As the summation of series, where the general 
term is given explicitly as a function of the index, thus 
resolves itself into the cases of integration treated of in 
Section II., it will only be necessary to give a few numerical 
examples; every expression integrated in that Section gives 
the sum of a series of which it is the general term. 

1. To find the sum of the first a? terms of any pro- 
gression of figurate numbers. In the r^ order the general 
term is 

_ a? (ar + 1) (^ + 2) .>. (^ + r - 2) 
1.2.3...(r-l) 

.'.^u.,^^^''^'"''^-^''^'''Kc. (Art. 45.) 

a?(a7 + l),..(a? +r- I) 



• '• Sff s= 



1 .2.3...r 



Similarly, for the series of inverse figurate numbers, ex- 

1 
cept when r = 2, that is, for the series 1 + -]^ + ^ + :^ + ••• + —• 

„ , 1 .2.3...(r- 1) 

For we have u^ = 



w(a?+ l)...(a? + r- 2)' 

_ 1 . 2.S...(r- 1) r A X 

.-. 2tt,+i « C - 7 -p , ^ ^ -r , (Art. 47.) 

(r-2)(j? + l) ...(a7 + r-2) 

2wi« C =»o, 

r -2 

10 
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^ (r lJ ^ l.2.5,.,(r-3) I 

and iS » . 

• r-2 

2. To find the sums of the squares and cubes of the 
natural numbers, 

t«, s /r* c= (^ — 1) or + ^ » 
.-. 2w,+i = ^ + ^ + C, 

2wi»0+C = 0; 

/r (ar + 1) (2^ + 1) 
.-. S, = - . 

Again, w, « a?* = a? (/r* — 1) + a? ; 

(j?-l)a?(a?-H)(a? + 2) ^(^ + 1) ^ ^ 

.-. 2tt,+i = ^^ + ^ + c» 

/P(J7+1) f(^-l)(^+2) 1 ^ f ^(^H-l) i' 

SimUarly, 1« + 8« + 5* + ... + (2a? - 1)* = "^^^^"^^ . 

1' + S* + 5' + ... + (2a7 - ly = 2a?* - a?*. 

3. To sum the series 

l«-2» + S^-4* + ... ±a?*; 



76 



Similarly^ 

4. To find the sum of of terms of the series 

111 



+ li— ^2 + ^— S + &C. 



2« - 3« 4» - 3« 6« - 3 



*** =* 7 — C5 3 =■ "7 r-7 ; > which falls under Art. 48. 

(2^)* -3* (2a? - 3) (2a? + 3) ' 

1 J? + l 1 6ar + 5 



18 (20? + 1) (2a? + 3) 6 * (2j? - 1) (2a? + l) (2a? + 3)* 

1 



and «y« = — * 
18 



Similarly, 



5 8 11 3a?+4 

+ H + &c. to a? terms = 2 - 



1.2.3 2.3.4 3.4.5 (a?:H)(a? + 2) 

5. 2. 4. 7 + 4. 7. 13 + 8. 13. 25 + 16. 25. 49 + &C. to a? terms. 
w,« 2* (3 . 2'"^ + 1) (3 . 2' + 1). 

Assume 2w,+i = J(3.2*-* + l)(3.2'+ l)(3.2'+i+l); (Art. 55.) 

.-. W,+i = ^(3.2' + l)(3.2'+^ + l){3.2'+*+l-(3.2-* + l)} 

21 
^ — A. 2*+^ (3.2*+ 1) (3 . 2'+* + 1) ; 
4 
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4 

.-. 2t#,^., = — (8.2'-i+ l)(3.2'+ l)(».2'+*+ 1) + C; 

4 5 
2^1 =» — .-.4.7 + C = 0; 
21 2 

4 « 40 

.-. *S, = — (3.2'-* + 1)(8.2»+ l)(S.2'+^ + 1) ---. 



2 4 8 i6 ^ 

6. -^ + + &c. 

3.3 3,9 9.15 15.33 



w, = 



(-2y 



K-2r-l}|(-2r^-l}' 
.'. aS, = - + -.7 — —r, — , (Art. 55.)^ and Sa.^- . 



16 21 2 26 /2\* „ 

7. . - + X • - 1 — &c. to a? terms. 

2.3.4 3.4.5 3 4.5.6 \3j 

16 + 5 (a; - 1) e 

then J ss - 3, and 



■• *'"*"(*+ 2) (^ + 1)^"^^' 



Similarly, 



19 1 28 1 39 1 52 1 



1.2.3 4 2.3.4 8 3.4.5 I6 4.5.6 32 
to cT terms, 
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(a? + 1) (^ + 2) ^^' 

1 lie 

cosdcos2 cos 2 decs 3d COS 3 d COS 4 d 

to w terms. 

1 



u^^ 



cos W COS (^ +1)0' 



^ tan (^ + 1)0 ^ 

2 «i,+i « . ^ + C, (Art. 53). 

sin 

^ tan ^ 

sin0 

tan (^ + 1) - tan 
smd 

9* 1 cos + 2 cos 20 + 3 cos 30+ &c., to w terms. 

«^*+i = (^ + 1) cos (<» + 1) 0, 

2 w,+i = (^ + 1) — ^ •^ 0^^^^ (^ + 1)^ (Art. 59.) 

2 sin - 2 sin - 

2 2 

_ (gy + 1) sin (^ + ^) cos (a? + 1) 

2 sm - I 2 sin - 

2 V 2/ 

COS0 



2^1= i +-^ 5r8+^ = ^' 



2 sm - 
2/ 



1 
• p— 



{'^-^) 
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2 8in- 28in'- 

2 2 



^ 2 I 

sm- / 



2 

Hence also if *y, = 1 co8*0 + 2co8*2d + Scos'S + &c. 
then 2 aS'^ » 1 (1 + cos 2 0) + 2 (1 + C08 4 0) + &c. 
a?(a?+l) ^ sin (2 ^ + 1)0 , /sina?0\* 



.-. S.^ 



4 



^ Sin (2 ^ + 1) , /sin a? C^y 
"^^ 4rii^^ * V8ine j ' 



10. tan""* i + tan"^ 1 + tan~^ ; + &c. 

1 + 1 + 1* 1 + 2 + 2* 1 + 3 + S* 

to w terms. 

1 , A^ 

u» = tan"* r = tan"* 



'« 



l + a? + d?* l+a?(a? + l)' • 

.*. 2 1^« = tan"^ A7 + C, (Art. 54.) 

.'. Sg^ S^x+i - 2^1 = tan"* («+ i) . 

4 



Recurring Series. 

90. We next come to the case where the general term 
is not given explicitly in terms of its index, but only certain 
relations between the consecutive terms, or these and their 
indices are expressed ; of this sort of series the most remark- 
able are Recurring Series. 

A recurring series is a series in which an equation of the 
first degree with constant coefiicients, holds good between a 
certain definite number of consecutive terms, in whatever part 
of the series they be taken. 
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For example, in the series 

3 + 5 + 9 + 17 + S3 + &c. 
we have 9 ^S.5 -^ 2,3, 17 = 3.9 - 2.5, &c. ; and in general 

The general equation of every recurring series is 

The series of coefficients which connects any term with 
the preceding ones is called the Scale of Relation. Thus 

is the scale of relation of the recurring series whose equation is 

91. A recurring series may generally be resolved into 
two or more geometric progressions. 

For if Oij 02, Os, •«•<{„ be the roots of the equation 

the complete integral of the equation of the series is 

Ug = Cj o* + Cg af + ... + Cn 0^9 (Art. 78.) 

Wj as Ci tti + ^2 ^2 + C3 Os + . . • + C„ 0„ 
t^ =s Ci Oj + C2 «2 + ^3 ®3 + • • • + ^n ^«> 



w, « Ci a" + C2 aj + C3 oS + ••• + Cn ajjj 
and the series consequently is transformed into 

Ci (oi + fl^i + ... + «J) + C2 (flfg + «2 + ... + aj) + &c. 
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92. In the particular cases in which equation (2) has 
equal, or impossible roots, the recurring series can no longer 
be resolved into geometric progressions ; for the complete in- 
tegral of the equation of the series becomes in those two 
cases, respectively, 

w# ■« (co + <?! ^ + ... + Cr-i of^"^) p' cos /r ft 

93. Hence, to find the general term of a recurring series, 
we must integrate the equation expressing the relation be- 
tween its successive terms, and determine the arbitrary con- 
stants by making the general term u, coincide with a sufficient 
number of given terms. When a series is known to be 
recurring, its equation may be determined by assuming it to 
be of the form (l); and .then forming a sufficient number 
of equations for finding the coefficients p^, p^^ &c. j?., by 
substituing the given terms in order. 

94. To find the sum of w terms of a recurring series. 
First, let its general term be of the form 

Ug « Cjflj + C2tt2 + ••• + Cj^a^y 

.*• <3, =s 2iW,+ i — ZWi = Ci . + C2 + ... 

flfi — I o^ — l 



+ c.. 



a^^-«n 



Secondly, let the general term be 

then since 

2a^o' = -7 r^ + ;t- - &c. (Art. 59.) 

o?+» K+'(a?+l) <+« 1 

2t^,+, = Co + Ci \-^ — i— -J + &c. ; 



•«■ 
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t 

.*. *y^ as'^Sw^+i — Swi, is determined. 
Thirdly, let 

Ug = p* (cq cosofO + Co sin w0) + p'.r (cj cos^d I gtjl^fl'fl^ ^ .^^^ x ^ ) 

+ p'a;* (cg cos a?0 + c/ sin 070) + &c., 
then each term to be integrated will be of the form 

af'p* co^ijeQ -{• a), 

the integral of which may be found by Art. 58, because 

2"{/)'cos(^0+a)} 
is always assignable. 

£x. 1. To find the sum of w terms of the series 

1 + 5 + 17 + 53 + &c. 

Let the equation be u^^^ +i^^«+i + 9^* = 0; 
.'. l7 + 5jp + g = 0, 53 + 17p + 5gf«0, 
which give jp = — 4, ^ = 3 ; 

Let Ug " o*; .-. a* - 4a + 3 = ; .-. a = 3 or 1 ; 

1 =: 3ei + Cg, 5 = 9Ci + Cg, 

.*. <?l ** ■§■> ^3 ~ "" ^> 

.-. w, = 2.3'-^-l; 

3* 
3 "■ 1 

.-. *y, = 3* - ^ - 1. 

Ex. 2. 2 - a - a* + 2o* - a* - a* + 2o* - a' - &c. 

_i (2a? - 2) 9r 



w- = 2a'" cos 



11 
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Ex, S. i+4+i» + 80 + S56+ ke^ to m tams. 
Here ii,+,-4«,+, -2«,-0; 

.-. 2>/6«,-(2+v^'-(2-\/6)', 
and 

1 ^y/e l-\/6 



Aptplicatkm of tbe Integial Caleohis to the SwmnMitinn of Series. 

95. As integrals are often expressed by Series, so, 
conversely, the latter may be represented by integrals; and 
it is often desirable to find the integral of which a proposed 
series is one of the developments, in order to subject it to 
the methods which we possess for calculating, at least ap- 
proximately, the value of any integral taken between assigned 
limits. We proceed therefore to notice one or two processes 
given by Euler for effecting this; they consist chiefly in 
performing certain operations on the series, by which it is 
transformed into another series which we are able to sum, 
or which is similar to the proposed one. 

96. Series which proceed according to the powers of 
some quantity U affected with coefficients consisting of factors 
in arithmetic progression either in the numerator or deno- 
minator, may be summed by the aid of the Integral Calculus, 
the denominators being taken away by diffecentiation and the 
numerators by integration. 

Ex. 1. Let ««- + — 4- — + &c. (to 00 ), 

2 3 4 ^ ^' . 



Ji t 



1 -t 
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d,(«0 1 



t (1-0*' 



= + 1 + C; 



• • o ^ — ~— »— ^— — . .^ 



1 -I? 



T, a + 6 + 26^ 
Ex.2. ^= -^-1- ;:^-f-.. 

a + /5 a + 2/3 



a + ^6 
+ 7i f + ... to «» terms. 



.-. fidt(8tP) = ... + (a + »6) ^+g"^ + ... 

Ex.3. « = 1+ ^t + ^_L^L_^/8 

a + i3 (a + /3)(a + 2/3) 

a (a + 6) (a + 26) ^ « ,. x 



?-i 



6^'** " ^Tfe + ^7^ + (a + /3) (a + 2)3) 

■ « (« + fe) 1^' . 

+ (a + /3) (« + 2/3) (a + 3/3) ' +*'^' 
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a (a + 6) 5+« 



+ 7 m-T — —7^ tP + &c. (1.) 



a 5-1 



tfi ^tPs; 



a — h 
hence, performing the di£Ferentiation, we find 

(fit''hf)dt8-^(a- at) s^a^ 
a linear equation of the first order for finding s. 

For the sum of m terms, we shall evidently have from 
equation (l.) 

1 (a + i3)(«+2^)... {« + (»»-!) ^1*^ /• 

Ex.4. ,.l+-.^.;l(^±i)M±l)^ 

1 7 1.2 7(7 + 1) 

. «(«+l)(a + 2)/3(/3 + l)(/3 + 2) ^ 

+ 7 T-;^ r r* + &c. (to 00 1. 

1.2.3 7(7+l)(7 + 2) ^ ^ "'^ 

(/* - ^?« + K« + i3 + 1) < - 7} dj« + a)3« = 0. 

97» When consecutive denominators have only a single 
factor in common, or none at all, the following is a conve- 
nient mode of summing the series. 



Ex. 



o (a + 6) (a + 26) (a + 26) (a + 36) (a + 4fe) 

+ &c. (to CO ). 



c::-.»'A^<^i'>«>:>UUBM^ ^ ■•fh3>_ 
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(a + ofb) {a + (a? + 1) ft} {« + (a? + 2) 6} a + wb 

2 1 



a + (a? + 1) 6 o + (a? + 2) 6 

resolving all the terms by means of this formula, by putting 
^ = 0, 1, 2, &c., we find 

12 2 2 2 

25^*5 = + r r + 77 - &c. 

a a-vh a + 2o a + So a + 4o 






Similarly, 



« (a + 6) (o + 26) (« + 36) 



+ &c. (to 00 ). 



(a + 36) (a + 46) (a + 56) (a + 66) 

98. We may obtain an expression for the sum of the 
series, 

l+^.d,0(«) + ^(?40(«r)p + ^^«?,|0(«)}« + &c., 

where («) is any function of ss^ by Lagrange^s Theorem. 
For if ytmz + af(j) (y), then 

y = « + -0 (») + — d, {0 («)}« + &c. 

where d <f)(y) must be obtained in terms of a^ and at, from 
the equation y — x-^of^ (y). Suppose, for example, 

0(y) = i(2^-Oj then d^<l>(y)=-y, 
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and y = ar+ -(y*- 1), ••• 1 -a^y = \/l -2«a? + ^; 

1* 1.2 2' vl-2»^r? + jr 

Hence if (l -2»a? + a?*)"i=l +Ziar+Z2^+ ... +Z,a^ + ..., 

then Z« = -^d; («*-!)". 

The function Z. has some remarkable properties, of which 
the following may be noticed. 

From the formula 



•/^V(a + 6a7)(c + ca?) vc6 [ e ) 

it is easily shewn that 
-V+* . P «x 1 / 2a a'\-i 1, /l+a\ 

.-. /*^(1 + Ziar+Z2aV+ ...)• (l + Zi-4-Z2^+ ...J 

aV 3 5 / 

Hence the definite integral of every term of the product 
forming the first member that is not independent of r, must 
equal zero ; 

... J (Z.Z,) = 0, and / (Z.)«--_— -. 

99' To find an expression for the sum of the reciprocals 
of the n*^ powers of the values of y, in the equation 

«-y + 0(y) = <>; 

where 0(y) denotes any function of y. 



• "*« 
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Let « - y + 0(y) = C(o, -y) (og - y).... (a. - y), 

«-y + <p{y) a^-y (h-y ^m-y 

Now developing the two members of this equation in 
powers of y^ the general term of the second member is 

and by Taylor's theorem the general term of the first member, 
putting l-dy0(y)=/(y), is 

But if /(y) = Jo + -4iy+ ... +-4„y"+ ... , 

smce = - +^ + ... +-7+1 + ••• > 

% — y % sr %^ 



^ f(y) ^o + ^«+ ... +^««" 

—as ... + ' 

x-y «f», 



= ... + :;^^^ 'y"+ ..., 



therefore the £:eneral term of is — — y", whatever be 

^ x-y «"+^ ^ 

the form of f(ss)^ if "^-^^ be restricted to such terms as in- 

volve negative powers of % only ; consequently the general 

term of ^.^lM is ^M^.y-, and therefore of 

hence equating coefiicients of y^ on both sides of equation (l), 



88 



Now putting for f(js) its value, the general term may be 
resolved into 



or 



» ^-,/£:>)^>)\ »+i^-./^'W 



(f 



/ 0— (af) 0- («) N _ « + l^_, /01(f)\ 






of which the first and last will be destroyed by the corres- 
ponding parts in the preceding and succeeding terms ; there- 
fore, changing n into w — 1, 



8 






But the second member of this equation is what Lagrange'*s 
Theorem gives for the development of y"% in the equation 
a? - y + (y) = ; therefore «.« is equal to the sum of those 
terms of the development of y" " which involve negative powers 
of %. 

100. By the help of the preceding proposition it may 
be shewn that Lagrange^s Theorem, applied to the solution 
of the equation 

« - y + (y) = Oj 

gives an approximation to the least root. By what precedes^ 
we have 

g-n jg* 1 js"-*-^ 1.2 V ;y*+^ / 

where each series is restricted to those terms which involve 
negative powers of ^; and the number of those terms in- 
creases with (w), and if n be supposed very great, then each 



^-n-r 
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ISeries may be taken ad infinitum. But whatever be the value 
of fij if both series go on ad infinitum, the value of the 
second member is 

as given by Lagrange^s Theorem, and which may be verified 
by actual multiplication; therefore when n is infinite, 

Umit of -i^ - 9% 

(f/^ denoting that function of » which is given for the develope- 
ment of ^ by Lagrange^s Theorem). 

But if Oi be the least of the roots Ou Og, Oj, &c., we 
have also (Theory of Algebraical Equations, Art. 156.) 

limit of — ^ «= a[y when n is infinite; 
.\ y^sOj, and y^aii the least root 



101. To find the sums of the series 



sin a 2 sin 2a S sin Sa ^ . 

4- — 5 — Ti + -15 — u + «^* 0° ^)) 



Aircosa A? cos 2a ApcosSa „ . v 

or the values of '^^?^, '^S-^^^^. 

First, we have 

—2^ — t-^ -:sc*'^+2(pcos^ + p*cos2ir+p'cos3a?+&c.)c"K 

l-2pcosa?+p« yr i- r / , 

12 
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and integrmtiog both ades from wmO to « ■ a^ and putting 
p « 1, we get 

* At -f- n* Ar + n* 

A; 



+ 2**S~ 



ife«+n« 



Again, integrating from ^mO to ^eSir, and putting 
p » 1, we find 

1 — «"*'"*' ^ k 

ir («-»** + 1) = ; + i^Sr -~ J - 2e-»»» »5" 



A» + n» k^ + n* 

Now subtract the foF^;oing result from this, and call 
the two sums we are in search of, i^i and St, then 

and changing the sign of k, 

___=_-_ 25-, -2*5,; 

therefore, adding and subtracting, 

^^ n&nna ^w c'^C'-") - e-*(«-«) 

* Af* + n* 2k ^^e-^ "^2A?' 



These formuke were first given by Poisson, and may be 
considered as embracing the chief results which have hitherto 
been obtained relating to the summation of series of the 
sines and cosines of multiple angles. 
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102. The definite integrals required in the preceding 
investigation may be found as follows* Suppose 



(1 -f) e-^'' 
1 — 2p COS a? +|>* 



y=. -:7Z „, „ . o > 



to be the equation to a curve ; and let it be proposed to 

find the limiting value of its area, from w = to w ^ a^ on 

the supposition that p approaches continually to unity. First, 

we observe that all the ordinates will be ultimately evanescent 

except those corresponding to a? = 0, Stt, 4^, &c., which will 

1 +© 
be very large, because for those values becomes a factor 

of the expression for y. We will begin by supposing that 
a<27r; and therefore the only ordinates that we are con- 
cerned with are those immediately succeeding that through 
the origin ; so that, making p = 1 - p, and then supposing 
J? very small, we get successively 

,»..4(l-,)sin«J P^' 

.-. Ty = 2tan--, 

Ts'i»=i = '^» making p^O. 

Next, suppose that the area is to be found from a? = 
to ^as27r; then, besides the area just found, there will be 
another portion immediately preceding the point for which 
,r = Sir ; to find this latter portion, put a? = 2 w - a/, and call 
the ordinate y \ 

1 — %p cos w ■\- p* 

therefore, the portion of the area immediately preceding the 
second limit = ^re"*'^'^; 
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103. As the following important theorems in Definite 
Integrals are based upon the preceding investigations^ we 
shall here give the proofs of them, though not strictly be- 
longing to the present subject. 

{^-P)A^) . /(-r)^.2i^«p«cosn(c-a?)^./(^)^ 



TT CI 

l-2pco8(c-^) — +p* 

a 

■•• "7/'y,=i - "7/ V(*) + a'-y '7/° cos n (c - ar) - ./(*). 

a 

Now, suppose c to be less than a, then for all values 
of iv between <r «= — a and ^ » + a, y is evanescent, except 
when 07 « c ; if therefore we write w^c + x, supposing z ex- 
ceedingly small, and then integrate with respect to %^ from 
x=^0 to AT = <y any small finite value, and double the result, 

and make p = 1, we shall obtain the value of "X^^y^ei* 
But making Jd « 1 - /^^ we have 

y ^ Pi^-f>)f(c^^) ^ M^ ultimately, 
p»+4(l-p)ffln«— p'+ — 

^afic) pair 

Arv 2a/(c) ^ , ^v 
TT pa 

.•.Typ-i-«/(c)5 

... 2a/(c)=-7/V(^)+2''S*'~7/"«>sn(c-*)-./(^r), (l.) 

a 

+ cos(c-a?) — . /(a?) + cos (c - 0?) — ./(^) +...}; 

a a 
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therefore, making a infinite, 



/(c)=:- limit of — ■"jT'^" Jcos(c-^) — +cos(c-.t) — + kc.]f(w), 
IT a a a 



or, /(c)«-7/"""X"*^"cos(c-a?)i8f./(a?), (Integ. Cal. Art. Il6), 

a theorem given by Fourier, and included, as we see, in the 
theorem (l), which was first given by Poisson. 



104. When in physical questions a definite integral arises 
whose value cannot be exhibited in finite terms, or in a form 
convenient for numerical calculation, the method of Quadra- 
tures is used as a substitute for Integration. This method 
consists in taking a series of values of the function to be in- 
tegrated, multiplying these by the differences of the cor- 
responding values of the independent variable, and adding 
together all the results. The sum of such results approxi. 
mates to the value of the definite integral, as the intervals of 
the independent variable are diminished. The method is 
equivalent to adopting, for the area of a curve which a 
definite integral gives, the sum of the areas of a system of 
inscribed or circumscribed polygons. 

Let the equal intervals into which the independent variable 
w is divided be taken equal to unity, and let f(at) be the 
function of of to be integrated between assigned limits a and, 6. 
Take a value of it /(a -f w - ^) at the middle of the n*^ 
of the intervals into which 6 - a is divided. If the intervals 
be very small, /(ooi) may be considered constant from 

w= a + n — I to ar = a4-n, or A^f(ai) 

may* be neglected ; and the value of the definite integral ia; 
2/(a + « - ^), where n is to receive in succession the values 

1, 2, 3, ... 6 — a. 
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For greater accuracy, suppose A'/(jr) to be constant bom 

w^a to w^by or ^f{a) 

to be neglected. Then 

/(a + n-i + sr)-/(a + n-i)+*A/(a-i-n.i) + -A«/(*), 

and the int^ral for the interval in question is 

^/(o + «-^ + «), from »=-^ to « = ^, 

The integral, consequently, from a to 6, is 
2|/(a + «_^)+lAY(«)}; 

n receiving the same successive values as in the former c^se. 
In calculating the value of a definite integral by this method, 
the intervals are to be taken smaller as the variation of the 
function is more rapid. 

Convergency and Divergency of Series. 

106. A series 

«*i + «2 + Ws + ... + «ii + ... ( to 00 ) 

is called convergent, if the sum «„ of any number n of its 
terms approaches continually to a finite quantity 8 as its 
limit, when n is indefinitely increased ; and divergent in the 
contrary case. 

106. When a series is convergent, the sum of any 
number of consecutive terms after the n^^ continually tends 
to zero as n increases. For 
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therefore, as n increases, the value of the first member con- 
tinually approaches to « — «, or zero, as its limit. This 
being true when ma 1, we see also that u^^i^ or the general 
term u^ continually tends to zero as n increases, or each term 
is greater than the following ; but this, although a necessary 
condition, is not sufficient to insure the convergency of a 
series. 

Thus in the series 

1111 

T + ;: + o + ••• + '"' 

12 3 n 

t^^ B ^ continually approaches to zero as n increases ; but 

1 I 1 

n+1 »+2 w+m 

is evidently greater than 

m 



n '\-m 



>\i if m = n; 



and therefore the sum of n consecutive terms after the n^, 
does not diminish indefinitely as n increases ; consequently, the 
series in question is divergent. 

107. In the geometrical progression 

a + aa? + tfa?*+ ... + aaf"* + ... (l.) 

1 -^p* 
1 -ar 

1 -af+* 
1 —m 



1 -a?" 



9 
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Henee ii m <\j when n is infinite «* « 0» 

a 



!-»' 



also «,«0, «.+,-«.«0; 

both wliidi results shew that the series is co&TergeDt. 

But if df > 1, then «. s aa^~^ increases indefinitelj with 
n, which alone shews that the series is divergent. 

Hence the geometrical progression (l) is convergent or 
divergent, according as j? is less or greater than unity ; and 
it may be used as the test of the convergency or divergency 
of other series. For if a proposed series can be shewn to 
have no term greater than the corresponding term of (l) 
when X <\^ then that series is convergent ; or if a proposed 
series can be shewn to have no term less than the correspond- 
hig term of (l) ;when a> \^ then that series is divergent. 

Thus in the series for .e, the base of the natural loga- 
rithms, 

111 1 

^ + 7 "•■ T~^ + , ^ ^ + -. + , 7 + ... 

1 1.2 1.2.3 [«- 1 

the terms which follow the ri^y viz. 

Ill 

+ &c. 



[» «+ 1 |nH-2 



are evidently less than the corresponding terms of the geome- 
trical progression 



the sum of which is 



1 

— + 
n 


1 
n 


1 
' n 


1 

+ 
n 


• -^ + &c., 


:h is 










1 


1 
1 - 


T "" 


1 


1 


n 


n - 


1 ' n-l' 



n 



97 

and continually tends to zero as n increases. Therefore the 
proposed series is convergent; and the error in taking the 
aggregate of its first n terms for its sum, is less than the 
quotient of the n^^ term divided by n - 1 . 

108. From the measure of convergency or divergency 
which a geometrical progression furnishes, we shall now pro- 
ceed to deduce one or two other tests as given by Cauchy, 
Caurs dCAnalyae AlgHmque. 

109. The series t«^ + %+... +!«.+ ••• is convergent, 

or will become so, if the superior limit of {u^*^ be less than 1, 
when n is infinite ; and divergent in the contrary case. 

1 
Let k denote the superior limit of {u^ when n is infinite ; 

and first suppose * < 1 ; also, let o be any magnitude between 

k and 1, so that k<a<\\ then when n is increased inde- 

finitely, {u^ cannot approach indefinitely near to k without 
finally becoming constantly less than a. Therefore it will be 
possible to take for n so large a value, that for that and all 
superior values, we may constantly have 

1 

Consequently, the proposed series will finish by always 
having its terms less than the corresponding terms of the 
geometrical progression 

a + o* + ... -f a" + o""*"^ + ... ; 

and as this series is convergent, a being < 1, it follows that 
the proposed series will end by being, i fortiori^ convergent. 

Secondly, suppose ft > 1 ; and take, as before, a between 

1 and ft, so that ft>a>l. Then, when n is indefinitely 

1 

increased, {u^)^ cannot approach indefinitely near to ft without 
finally becoming constantly >a; we shall therefore be able 
13 
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1 

to satisfy the condition (w^)»>a, or u^>a\ by taking n 
sufficiently large; and consequently we shall always find in 
the series Wj + w, + ... + w, + u»+i + ... , an indefinite number 
of tenns greater than the corresponding terms of the geo- 
metrical progression a + a' + ... a" + o*"*"* + ..., which is di- 
vergent, a being > 1 ; and therefore the proposed series will 
end by being divergent. 

110. It may be shewn that if, as n increases indefinitely, 
tt« remains positive, and the ratio --^ continually tends to 

become equal to a finite quantity k as its limit, the expression 

(u^ continually approaches to the same limit. Hence, the 
test of convergency or divergency in the last Article, may be 
changed into the following, which is more convenient in its 
application ; if the limit of the ratio of u^^i to u^ when n is 
infinite be less than 1, the series is convergent ; and divergent 
in the contrary case. 

If A « 1, this test gives no result in either form. 

111. Suppose that the series t«i + «^2 + ^ + •*• consists 
of both positive and negative terms; then if Vi, v%y Ds, &c. 
be the numerical values of these terms, so that u^ = Ji v^ , 
Ui^ ^ f>^i &c., it is evident that the sum of the proposed 
series can never surpass that of the series t^i + Vj^ + t9s+ ...; 
if therefore, the latter series be convergent, that is, if 






<i, 



the proposed series will be convergent ; or, if the latter series 
finish by having terms greater than any assignable magnitude, 
that is, if 



(— ) 



>i, 
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the same thing will happen to the proposed series, which will 
consequently be divergent Hence, the above test is ap- 
plicable to series consisting of both positive and negative terms, 
provided we use the numerical values of the terms without 
regard to signs; and it fails, as in the preceding case, when 

112. Let ao + aifl? + 02^+ ... +a«^ + ..• (!•) 

be a series arranged according to positive and ascending 
powers of the variable w; the coe£Scients being positive or 
negative; then, by what has been proved, this series will be 
convergent or divergent, according as kw, (where 



*-(«.) 






9 



when n is infinite,) is numerically less or greater than 1. 

Hence for all values of w between - - and + t the series 

k k 

will be convergent, and for all values of 4? beyond those 
limits it will be divergent. 

1 1.2 1.2.S 



Here -^ * = - !» when n is infinite; 

a« » + 1 

therefore the series is convergent or divergent, according as ^ 
lies between + 1 and -^ 1, or without those limits. 



Ex. - + - + T + &«•> 

12 3 



a" 
n 



• 
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e K a. when n is infinite ; 

«, n + 1 

therefore the series is convergent or divergent, acoording as 
a<or> I. 

Ill 1 



».+! 1 



0, when n is infinite. 



therefore the series converges. 

113. The series fii + f% +.•• + «, + ... is convergent, 
or will become so, if the inferior limit of be > 1, when 

n is infinite; and divergent in the contrary case. 

loff tf 
Let 4p denote the inferior limit of ^ * when n is in- 

finite, and first, suj^se A > 1 ; also let a be any number 
between k and 1, so that h> a> I. Then when n is in- 

- log — 

definitely increased, ' ■ * or its equal * cannot approach 

logl ***«" 

n 

indefinitely near to h without finally becoming oonslantly 
greater than a. Therefore it will be possible to take for n 
80 large a value, that for that and all superior values we 
may constantly have 

log — > a log n > log n , or — > »% or «, < -j • 
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Consequently the proposed series will finish by always having 
its terms less than the corresponding terms of the converging 
series (Ex. Art. 11^. 

Ill 1 1 „ 

I Sr fr n (» + 1)* 

and therefore will itself end by being, a fortwriy convergent. 

Similarly, if Ap < 1, it may be shewn that the proposed 
series will finish by being divergent. 

114. The series «o + ^i + ^s + tfs + ^4 -^ ^ + &<^- iX^y 

each term of which is less than the immediately preceding 
term, 

and the series Uo-^-ftu^-^ 4i% + 8«7 + l6«i» + &c. (2)^ 
are convergent or divergent at the same time. 

Suppose series (1) to be convergent and that its sum « a, 
then 

4W8<2«j| + SUB, 

8«7 < 2fl4 + 2f% + 2«e + 2«f7 



consequently series (2) is convergent. Next, suppose series 
(1) divergent, then 
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and is therefcve, diTergent. 
Ex. Let series (l) be 



1111 



- + - + - + jj+kc. (3), 



then series (2) is 

a geometric progression convergent when m > 1, and di- 
vergent in the contrary case; consequently series (3) will be 
convergent if m > I, and divergent if m « or < 1. 

115. The series Ui — f«, + Us — <^4 + &c., is convergent, 
if the numerical value of the terms decreases without limit. 

For, by writing it in the forms 

«i - Wj + (ms - ttj + &c. ; 

we see that it is >f«i — t^^ ^^^ <^i9 aQ<1 therefore is con- 
vergent. 

Thus, the sum of the series 

1111. 
+ - — + &c., 

12 3 4 ' 

lies between 1 and -; also the series 

2 

1 1 1 1 « 
+ + &c., 

1" 2" 3" 4* 

is convergent for all positive values of m* 
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116. It is also shewn in the work from which these tests 
of convergency are taken, that two converging series all whose 
term9 are positive, will by their addition or multiplication 
produce new converging series, whose sums result from the 
addition or multiplication of the sums of the former. 

117* An approximate value of 2i*, will be obtained by 
taking the aggregate of the converging terms only, in the 
series for 2u« involving ^t^^, t^«, and the differential coefficients 
of Ug^ and the error will be less than the last of the con- 
vergent, or the first of the divergent terms. 

We have by Art. 63^ (omitting the index of «, in the 
second member,) 

2tt, =5 C 4- Lu - 4 w + — - dgU - -r^ dlu + *.. 

1.2 \4t 

where 2 Jg, = ^'"^' di^^'u — ^^^i^ dJ«+3^ + &c. ; 



2» + 2 ' 2n + 4. 



or since in general, by Art. 67, 



B^n-l f 1 1 1 1 

(2irr+* (4»r)*"+« (6fl-)^+« ^ "• J "• " 



— &c. ; 



fH^ 



T^^. ^P5r>T^^^J J^^^ 



,•? ■W^." • 



' ^^F 



^^^ 



.1' 



4f.P^m^fimr*t.' 



'•('• 



» • 









I 






im^M.d^^m), 



mfi^^f ft hmtff M jffOfff/tTf iht tfm flHilnplKd by flBSoivx 
l# UftffH'pmtfff htAnff uhrtudy iotrodtieed in eqoatioD (l) ; 






+ iic 



\dT*'ui 



'''-^-'/{(54-«^(i^ 






rfj-fi, 



whtoh Unl q\i«iUlty U^ lH>twcon 



*!LLLj^^»«. 
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if these be the last of the convergent and the first of the. 
divergent terms respectively, of the series for 2w,. Conse- 
quently, the sum of all the diverging terms in the series for 
^u,f is less than the last of the convergent or the first of the 
divergent terms. 

Interpolation. 

118. When a series of values of a quantity is ob- 
tained either by observation, or laborious calculation, it 
is of great importance to be able to insert other values 
between them, such as would have resulted from a similar 
observation or calculation, without the labour of performing 
these. This is the object of Interpolation; and in this the 
Calculus of Finite Differences finds one of its chief uses. 
More strictly, Interpolation of Series is the inserting among 
the terms of a given series, new terms subject to the same law 
as the first. In doing this, the terms of the series are con- 
sidered as particular values of the function which expresses 
its general term, corresponding to a given regular succession 
of indices; and it is the business of Interpolation to dis- 
cover that general term ; or at least to assign such a function 
of the index as shall represent the given series of values, 
and, approximately, all intermediate values. The problem 
thus requiring us to assign the analytical expression of a 
function from a limited number of its numerical values, is 
plainly indeterminate; it is the same as to form the equa- 
tion to a curve which shall pass through a limited number 
of points, whose abscissae represent the values of the inde- 
pendent variable, and the ordinates those of the function, 
without giving the species of the curve ; which, as is evident, 
may be done in an infinite variety of ways. But if the 
given terms are numerous, and near to each other, the ex- 
pression for the general term, within the limits of the given 
quantities, may be found to a great degree of accuracy. 

119. There are two principal cases to be considered; 
first, when the given values of /(a?), namely, 

/(^i)» /(^i 4- A), /(a?i + 2A), ... f{wi + (« - 1)^1 , or 

Wl, U29 W3, •«., W„5 

14 
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as we shalL write them, correspond to values of the inde- 
pendent variable, a?„ Wi -f A, a?i + 2 A, &c., j?, + (« — 1) A, in 
arithmetical progression. And secondly, when the given 
values /(a?,), /(a?,), .../(ar,)* or «„ «,, ... ti„ correspond to 
values iVi, ^2> ... ^»» of the independent variable, not obeying 
any assigned law. 

120. Having given w,, «„ Wj, &c. «^, n values of a 
function /(^), corresponding to the n values of the indepen- 
dent variable a?!, a?i + A, ... a?j + (n - 1) A, to find an ex- 
pression for any intermediate value /(j?i + A). 

Since /(a? + nh) -/(a?) + n A/(4r) + ^^^^^^ A*/(a?) + &c.. 

(Art. SS) changing x into d^j, and then replacing nh by Ap, 
we get 

k k(k — h) 

/(ar, + *)-«, + - A«, + -Jjj^ A»tt, + &c. 

k(k-h)...{k-in-»)h\ 

^ |«-1.A'-' ^ "" 

in which, since 

(2). 



1 .z 

if we make A » 0, A, 2 A, &c., (n — l) A, the second member 
assumes the n values u^ t^g, ... t^„; and not only this, but if 
we assume for k any value whatever between and (n - 1) A, 
we shall obtain the value of the function corresponding to 
that value of the independent variable. 

121. In applying the above formula, the simplest mode 
is not to calculate Ai^i, A^t^i, &c. by equations (2), but by 
continued subtraction of the given terms; that is, we must 
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write down the series of given values t«|, t^, ... u^, and sub- 
tract each from the succeeding one ; next subtract each of 
these differences from the succeeding difference ; then perform 
the same operation upon the new differences ; and so on, till 
the process terminates ; the first terms of these series of dif- 
ferences are the values of Awi, A*Wi, ..* A*"*t/i. Unless 
the terms of the given series by continued subtraction lead to 
a constant difference, the expression for f(x^ + k) will have 
as many terms as the given series has. 

Ex. 1. Having given the values of sin SO®, sin 31®, sin 32% 
sin33®, to find sin (30® + ftO' ^ being between and 180'. Here 



Ui =! .5 


A 






u^ = .5150381 


150381 


A* 




^3 « .5299193 


148812 


-1569 


A" 


«4 = .5446390 


147197 


-1615 


-46; 



.-. Awi = 0-0150381, A* Wi «- 0*0001569, A' Wi «- 0-0000046, 

. , A , /v * . A? (A? - 60) , . 

.-. sm(30® + A')«.5 + -Awi + -^^-g^A*f*i 



A? (A? - 60) (A; - 120) 
2.3.6tf» 



A'w^. 



If ft a 20, it will be found that sin 30®. 20' = .5050299, 
which is too large only by a unit in the seventh place of 
decimals. 

Ex. 2. Having given 
log 3.14 =« 0.496929, log 3.15 =. 0.498310, log 3.16 = 0.499687, 
log 3.17 « 0.501059 ; shew that log 3.14159 - O.497149. 

Ex. S. Having given «^i9 1^9 ^£3, t^4, four right ascensions 
(declinations, longitudes, &c.) of the moon at intervals of 
12 hours, to find its value t hours after the time corresponding 
to the second value. Here A » 12, ft « 12 + ^^ and if the 
required right ascension -> t^, ^ ^, then 



lOB 



tm ^«'a *-A*^-A^ — -►-a* — : --A'f — ■ , 



« <5 ;I2 « li; 6 \\± 



I 



mhem 6er doped im poven of — . 

122. Betmtea ererj two coDsecutiTe terms of a 
MtJex, to interpfAzte aar Damber ci eqindisUmt 



Let f^i, tt^, tt», fcc^ «. be the fpTen aeiicsy and left m ^ 1 
be the number of equidistaDt terms to be inserted bel w eeia 
eftry two cfm%ecutite terms; tbeb the new series will be 

If therefore r^^^i denote the r + 1*^ term of this 
we hare 



^r^\ 



»^-/0 + ^)» 



or t^^^, « tt, + — Awi -»• T- A'«i + «c. 

i» 1.2 m 

Hence^ taking r from l to m- l, we get the terms in- 
serted between Ux and u^ ; next taking r from in + 1 to 
iim^lf we get the terms between U2 and «3; and so on. 
The differences At^M A^e^M Sec, are to be computed by con- 
tinued subtraction as in Art. (121) ; and the series for v^^^ 
will have as many terms as the proposed series has, unless 
thoifc terms by continued subtraction lead to a constant dif- 
ference. 

Ex. To insert three equidistant terms between every 
two consecutive ones of the series 1, 7, 15, 28, 49, &c. 
Here m -i 4, and 

A, 6, 8, 13, 21, 

A" 2, 5, 8, 
A' 3, 3, 



I. 
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6r f(r-4) r(r-4)(r-8) lg8 + 192r-4i^ + ^ 
'^' 4 16 128 128 



, . . . 317 504 695 ^ 1113 ^ 

and the series is 1, , • — , 7, • &c. 

128 128 128 128 



123. The formula of Art. 120 may be presented under a 
different form by changing k into «v — o^i, which gives 



n l,2.n 



{w - w^{w - a?i— A) ... {«!r - ^1 - (w - 2)A} ^ 



where /(^) is a function of of which, as m assumes the n values 
Wiy Wi + hf &c., ^j + (^ - l) A, successively assumes the cor* 
responding values u^^u^^ ... u^\ and for any other value of 
w within, or not far beyond, the limits x^ and w^ + {n^\)hj 
it gives the value of the corresponding interpolated term. If 
we put A s 1, the formula is adapted to the case where the 
increment of the principal variable is unity. 



124. In any series of consecutive equidistant values of 
a function, where one is deficient to insert that one. 

Let t^j, t^, t<s9 &c., u„ be the values of the function 
corresponding to the values /V|, j?} + A, ... a?j + (« — l) A for w. 
Then assuming that A"'"^Wi = 0, or that the (n-2)*^ dif- 
ferences are constant, which will almost always be the case 
in tabulated results, we have 

(w-l)(»-2) 
A-^Wi = w« - (»- l)w,-, + — W,»-8 - &C. ± M, « 0, 
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I 
I 

Ml cqiMtioD of the fint degree from whidi maj one at tbc 
▼does as Ur nuj be foundy if the rest be known. Hkving' 
thus completed the sjstem of Tallies, we may intcrpolsite mnj 
intennediate term /(dr, -^ k) bj the method of Art. ISO. If 
two values out of n mre deficient, then we most 



A*"*«, « 0, A'-'tf, - 0, or 

, ^ (n - 2) (ft - S) 
«..i-(«-«)«.-t+ r-T «._! - fee i 11, « O, 

(» - 2) (« - 3) 

which equations will suffice to determine any two of the values 
in terms of the rest. In the same way any number of 
deficient terms may be inserted^ 

Ex. 1. Given the cube roots of 121, 122, 124^ 125^ to 
find that of 12S. 

u,- 4.946088, tis » 4.959675, t^4 « 4.986631, fh^5, 

A*tti « «5 - 4^4 + 6ws — 4f^ 4- t^i B ; 

.-. «B - g {4 (i/j + u^) - ttj - w,} « 4.973190. 

Ex. 2. Given t^i, t^g, f^s, u^^ to find f^, t^^. 

A*|^, as t^5 - 4^4 + 6w, — 4f^2 + t«i B 0. 

A*Wji « Wj - 4w» + 6^4 - 4«8 + 1«2 " ; 
.-. W8=i5 (-3Wi + 10m, + 5wB-2We), 



y 



Ill 

K 125. We next come to the case where the given values 






correspond to values o^i, w^^.^.iv^^ not obeying any assigned 
law; and it is required to determine a rational integral 
function of n — 1 dimensions, /C^), which shall assume the 
n given values u^^ u^^.u^^ when for w the values w^ 
d^29 ^39 &C.9 ^„ are successively substituted. 

Since /(^) is of (n — 1) dimensions, we may assume 
f(w) C, . C, . <7, 

+ + •• • + 



Now make w^w^y w^y &c., a?^, successively; and observ* 
;. ing that the corresponding values of the first member are 

f^,9 t^S9***^ll9 ^6 get 

Wj =» Ci {wi - j?2) (<rj - ^s) ... (^1 - ^J, 

Ml 

u^ = C2 (a?a - a?i) (<r8 - .2?3) ... (d?8 - a?.) 



(«r - ^2) C^ - ^3) ••- C^ - w^) 



• /(^)«««i 



(a?i - a?2) (^1 - ^3) •.. (^1 - ^n) 



{W — a7i) (O? — a?3) ... (a? — W^ 

(j7 — a?j) (a? — ^2) ••• (^ — ^%-\) 
which is Lagrange^s Theorem for interpolation. 
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Ek. To find a function of w which, when a? s i, 3, 6, 12, 
shall assume the values 1, 7, 10, - 8. 

•^^ ^*"" 2.5.11 ' 2. 3.9 

(a? - 1) (a? - 3) (/I? - 12) (jv - 1) (a? - 3) (^ - 6) 

— ' — 4 . — • 

S.S 11. 9.3 

126. To determine the maximum or minimum value of 
a function, from three of its values near its maximum or 
minimum, and the three corresponding values of the indepen* 
dent variable. 

If Uxy t^29 ^3 ^ ^^ given values of Uj and w^^ w^, w^ 
those of ^, we have 

(/p-a?2) (a?— ti^s) (^-»Vi) (^—^3) (w-Wi) (d?— d?2) 

w=«^i7 77 rrz +^^27 77 r + ^7 r? c- 

(^1-^2) (^1-^3) (^«-^i) w-^s) («»3-^i) (^3-^2) 

Hence, putting d^u = 0, we find 

Ui (a?2 - ^3) (2^ — ^2 "" ^3) + % (^8 — ^1) (2^ — a?i — a?,) 
+ % (j^i — a?2) (2^ — ^1 — a?2) = ; 



^ = 



2^1 (^2 - ^3) - 2% (^1 - ^3) + 2W3 (a?i - ^72) ' 



the value of w at which «^ is a maximum or minimum. This 
formula is useful in various Astronomical problems, as for 
instance, to determine the meridian altitude of a heavenly body, 
when an observation exactly on the meridian cannot be ob- 
tained. 
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